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ABSTRACT 


Fractal  Geometry  is  a  recent  synthesis  of  old  mathematical 
constructs.  It  was  first  popularized  by  complex  renderings  of  terrain 
on  a  computer  graphics  medium.  Fractal  geometry  has  since 
spawned  research  in  many  diverse  scientific  disciplines.  Its  rapid 
acceptance  has  been  achieved  due  to  its  ability  to  model  phenomena 
that  defy  discrete  computation  due  to  roughneaa  and  diaeontinuities. 
With  its  quick  acceptance  has  come  problems.  Fractal  geometry  is 
a  misunderstood  idea  that  is  quickly  becoming  buried  under  grandi¬ 
ose  terminology  that  serves  no  purpose.  Its  essence  is  induction 
using  simple  geometric  constructs  to  transform  initiating  objects. 
The  fractal  objects  that  we  create  with  this  process  often  resemble 
natural  phenomenon.  The  purpose  of  this  work  is  to  present  fractal 
geometry  to  the  graphics  programmer  as  a  simple  workable  tech¬ 
nique.  We  hope  to  demystify  the  concepts  of  fractal  geometry  and 
m^e  it  available  to  all  who  are  interested. 
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I.  AN  INTRODUCTION  TO  FRACTAL  GEOMETRY 


A.  MATHEMATICS  AS  A  MODEL  FOR  OUR  UNIVERSE 

The  various  branches  of  mathematics  have  through  time  developed  as  a 
response  to  the  need  for  more  detailed  models  to  describe  new  developments, 
both  technological  and  philosophical.  This  was  true  when  Newton  developed 
calculus  and  also  true  during  the  late  1800*8  through  the  1020’s  when  a  schism 
developed  between  the  classical  mathematicians  and  some  brilliant  mnovaUve 
thinkers. 

1.  The  Mathematical  Crises  of  the  Early  10th  Century 

One  of  man's  greatest  strengths  is  his  ability  to  question  his  surroundings 
and  beliefs  and  through  this  questioning  develop  new  insight  and  innovation. 
Most  mathematical  systems  are  developed  for  use  in  applications.  Man's  natural 
Inquisitiveness  (dten  leads  him  to  develt^  his  systems  beyond  the  application  and 
into  abstract  theory.  This  theory  drives  him  to  investigate  the  applications  and 
often  yields  direction  for  new  discoveries  that  were  not  previously  foreseen  or  that 
defy  intuition. 

Georg  Cantor  (184S-1018)  was  the  most  notable  of  a  number  of 
mathematicians  who  questioned  the  basic  precepts  of  mathematics  and  developed 
the  modem  set  theory.  Some  of  Cantor's  discoveries  seemed  to  invalidate  many 
of  the  long  held  beliefs  of  mathematics.  Cantor  and  his  peers  became  deeply 
involved  in  controversy  over  their  finding.  Their  discovery  of  functions  which 
seemed  to  violate  the  basic  rules  of  geometry  and  calculus  were  deemed  as 
mon$ter9  and  unworthy  of  consideration  by  reasonable  men  because  they  lacked 
usefulness  to  any  application  then  known  (Ref.  l:pp.  6).  These  new  concepts 
would  remain  in  the  arena  of  pure  theoretical  mathematics  until  science 
developed  to  a  point  where  the  old  modeb  could  no  longer  adequately  describe  its 
processes  and  would  look  to  the  new  mathematics  for  a  new  perspective. 

It  was  from  these  discoveries  that  Ptmetai  Gtomtiry  was  bom  (Ref. 
i:Chap.  2).  It  will  be  seen  in  the  following  chapters  that  fractal  geometry 
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is  a  synthesis  many  of  the  concepts  which  developed  from  the  mathematical 
schism  of  the  10th  century,  most  notably  set  theory  and  topology. 

2.  What  is  a  Mathematical  Model 

Reference  2  deitnes  a  mathematical  model  in  the  following  fashion  [Ref. 
2:pp  1>3]: 

A  mathematical  model  is  a  .mathematical  characterisation  of  a  phenomenon  or 
procen.  It  has  three  esMntial  parts:  a  process  or  phenomenon  which  'S  to  be 
modeled,  a  inathematical  structure  capable  of  expressing  the  important  proper* 
ties  of  the  object  to  be  modeled,  and  an  expliat  correspondence  between  the 
two. 

Although  the  phenomenon  of  interest  need  pot  be  taken  from  the  real  world, 
they  usually  are.  The  real  world  component  is  described  quantitatively  by  su^ 
things  as  parameter  values  and  at  which  time  things  occur. 

The  second  component  of  a  model  is  an  abstract  mathematical  structure.  In  it* 
self,  the  structure  is  abstract  and  has  no  intrinsic  relation  to  the  real  world. 
However  because  of  its  abstractness  it  can  be  used  to  model  many  different 
phenomena.  Every  mathematical  structure  hiu  an  associated  language  for  mak¬ 
ing  assertions.  If  the  mathematical  model  is  successful,  the  language  of  its 
mathematical  structure  can  be  used  to  make  assertions  about  the  object  being 
modeled. 

The  third  component  of  a  model  is  a  speciilcation  of  the  way  in  which  the  real 
world  is  represented  by  the  mathematical  structure,  that  is,  a  correspondence 
between  the  elements  of  the  first  component  and  those  of  the  second. 

3.  The  Euclidean  Model 

When  using  mathematics  to  describe  man^soade  objects,  the  Euclidean 
model  (standard  Euclidean  geometry)  b  usually  satlafactory.  Its  structure  b 
simple  and  pure,  which  appeals  to  an  engineer’s  nature.  But  as  technology 
expands  and  we  need  to  describe  processes  that  are  not  well  hcAoued,  we  need  to 
develop  a  geometry  that  can  adequately  model  our  process  within  a  cerUin 
closeness  of  scale. 

No  model  can  completely  describe  a  natural  object  because  nature  docs 
ncA  follow  the  man-made  rules  that  vre  impose  on  our  model.  But  at  a  given 
scale,  the  model  (if  it  is  accurate)  can  describe  the  object  with  enough  precision 
to  be  of  help  in  constructing  it.  Engineem  use  the  geometry  of  a  straight  line  to 
describe  a  wall  but  this  wall,  when  viewed  closely  enough,  is  not  straight  at  all. 
Thb  is  of  no  matter  to  the  engineer  because  his  model  is  accurate  for  his  scale  of 
reference. 


8 


B.  FRACTAL  GEOMETRY 

One  man  who  saw  a  need  for  a  new  geometry  was  Benoit  B.  Mandlebrot.  He 
felt  that  Euclidean  geometry  was  not  satisfactory  as  a  model  for  natural  objects. 
To  anyone  who  has  tried  to  draw  a  picture  of  a  nonregular  object  (such  as  a  tree) 
on  a  computer  graphics  screen,  using  the  Euclidean  drawing  primitives  usually 
provided,  this  is  an  obvious  statement.  The  strength  of  Maidlebrot’s  finding  was 
his  research  into  the  findings  of  the  earlier  mathematicians  and  the  development 
of  a  practical  application  their  theory.  Mandlebrot  coined  the  term  Frmetat  to 
describe  a  class  oS  functions  first  discovered  by  Cantor  (Cantor’s  dust),  Koch  (the 
Koch  curve)  and  Peano.  He  showed  how  these  functions  yield  valuable  insight 
into  the  creation  of  modeb  for  natural  objects  such  as  coastlines  and  mountains. 
Mandlebrot  popularized  the  notion  of  a  fractal  geometry  for  these  types  of 
objects.  Although  he  did  not  invent  the  ideas  he  presents,  Mandlebrot  must  be 
considered  important  because  of  his  synthesis  of  the  theory  at  a  time  when 
science  was  reaching  out  for  new  more  accurate  models  to  describe  its  processes. 

C.  GOALS  OF  THIS  RESEARCH 

There  are  two  approaches  that  can  be  taken  in  the  investigation  of  fractal 
geometry  and  computer  graphics. 

>  To  view  the  computer  as  a  tool  to  enhance  the  investigation  of  fractal 
geometry. 

or 

•  To  view  fractal  geometry  as  a  tool  to  enhance  the  realism  of  computer 
graphics. 

This  research  will  take  the  later  ^iproach'.  It  is  designed  to  investigate  the 
mathematics  of  fractal  geometry  and  to  show  its  application  to  computer 
graphics.  I  hope  to  be  able  to  tame  the  subject  of  fractal  geometry  by  makini  its 
mathematics  and  technittuc  accessible  to  the  average  computer  scientist. 


'  When  MAttdeibrot  Uioie  th«  ibrtew. 


U.  THE  MATHEMATICAL  BASIS  OF  FRACTAL  GEOMETRY 


This  chapter  is  a  brief  introduction  to  the  mathematical  foundations  that 
underlie  the  theory  of  fractals.  Little  technique  currently  exists  for  the  practical 
application  to  attain  complete  mathematical  rigor  when  using  fractiU  functions 
(t.e.  it  is  very  difficult  to  prove  that  a  set  is  fractal).  This  causes  the  non^ 
mathematician  to  accept  much  of  what  he  does  with  fractals  on  faith.  It  is 
instead  important  to  understand  the  theory  intuitively.  This  can  be  gained  by  a 
cursory  look  at  the  mathematical  foundations  for  fractals. 

A.  PRELIMINARIES 

A  complete  definition  of  fractab  is  given  later  in  this  chapter  but  before  we 
can  understand  that  definition,  we  must  establish  a  foundation  in  set  theory. 
Fractals  were  discovered  in  set  theixry  and  topology.  They  can  be  considered  as 
an  outgrowth  of  investigations  into  these  related  fields. 

1.  What  is  a  Set 

A  set  is  defined  in  (Ref.  3;pp.  U]  in  the  following  fashion: 

A  set  is  formed  bv  the  grouping  together  of  single  objects  into  a  whole.  A  jset  is 
a  plurality  thought  of  as  «  unit.  We  can  consider  these  statements  as  espository, 
as  references  to  a  primitive  t  mcept,  familiar  to  us  ail,  whose  resolution  into 
more  fundamental  concepts  would  perhaps  be  neither  compeienl  nor  necessary. 

We  will  content  ourselves  with  this  consiructian  and  will  assume  as  an  axiom 
that  an  object  A#  inherently  determines  certain  other  objects  a,  k,  c,  ...  in  some 
undefined  way,  and  vice  versa.  We  express  thU  rdatlon  by  the  words:  The  set  .V 
consists  of  the  objects  a,  t,  t,  ... 

This  defiuition  is  iuteotioually  vague  to  allow  the  set  to  become  the  basic 
building  block  for  all  matheoaatical  constructs. 

2.  Some  Set  Thcorcijc  Concepts 

This  section  presents  some  background  definitions  for  concepts  used  in 
the  body  of  this  chapter.  The  reader  is  directed  to  the  references  for  a  detailed 
explanation  or  proof  (Ref.  2:Cbap  Sj. 

DefimitismM: 

Cardinality 

Two  sets  5  and  T  are  said  to  have  the  sanv  number  of  rirrtiento,  or  to  have  the 
same  cardimafitg,  if  there  is  a  ooc»oae  fuMtton  /from  5  to  T. 
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A  Set  5  is  eaid  to  be  finite  if  S  hat  the  same  cardinality  as  or  if  there  is  a  po¬ 
sitive  ipteger  n  sucii  ^at  S  haa  the  same  cardinality  as  ).  Otnu* 

wise  S  ts  smd  to  be  mfimte. 

Countability 

A  set  5  is  9i!iid  to  be  enmnUMc  if  S  has  the  same  cardinality  as  a  subset  of  N,  the 
set  of  positive  integers.  Otherwise,  5is  said  to  be  wummMe. 

PnfpMftietu: 

a) .  Any  subset  of  a  finite  set  is  finite. 

b) .  Any  subset  of  any  countable  set  S  is  countable. 

c) .  The  set  of  natural  numbers  N  is  countable. 

d) .  The  set  of  rational  numbers  Q  is  countable. 

d) .  The  set  of  irrational  numbers  is  countable. 

e) .  The  union  of  a  countable  coltection  of  countable  sets  b  cou  .liable. 

f) .  The  set  of  real  numbers  R  is  uncountable. 

3.  Some  Topological  Concepts 

This  section  presents  some  topological  background  concepts  used  in  the 
body  of  this  chapter.  The  reader  is  directed  to  the  reference  for  a  detailed 
explanation  (Ref.  3|. 


CnmtepU: 

Metric  Space 

muttie  is  a  set  tn  which  we  have  a  measure  of  the  closeness  or  proximity 
of  two  elements  of  the  set,  that  it«  we  have  a  dliunce  aefUied  on  the  seV.  For  ex* 
ample,  a  metric  on  S*  would  be  the  pythugorcan  meUtc: 


Covering  a  Sat 

Let  A’  be  a  topological  space  and  5  a  subset  of  A.  A  csser  cl  the  set  5  Is  exsetSy 
what  Its  narne  impUes,  a  coUcclion  of  subeete  of  X  which  cover  S,  that  U,  whose 
union  cootams 

4,  What  is  a  Funcllon 


A  function  is  defined  in  (Ref,  2:pp.  103*194]  as 

A  function  from  aseiAtoaictRisa  rule  which  speciftci  an  eiemcnl  of  R  for 
each  ekment  of  4. 

Let  4  and  R  be  set*.  A  Unetion  (or  map,  or  ttentfotmetio**)  /  from  4  to  R, 
denoted  /  :  A  B  ,  is  a  relatioa  from  .4  to  R  such  that  for  every  a  whkh  Is  an 
element  cf  4.  there  exists  a  uoiquc  h  in  R  such  that  <a,4i-  is  an  eletmmt  of  /. 
We  write /fa/ - 


If  /  is  function  from  A  to  B,  then  A  is  called  the  domam  of  the  function  /  and  B 
is  called  the  eodomtdn  of  /. 

To  completely  define  a  function  we  must  specify  the  domain,  the  codomain,  and 
the  value  f(x)  for  each  possible  argument  x. 

Functions  can  be  viewed  as  a  specification  of  a  method  to  describe  the 
creation  of  a  set  from  other  sets  using  some  agreed  upon  mathematical 
symbolism.  The  functions  can  yield  powerful  results  when  the  target  set  (co¬ 
domain)  is  complex  and  not  easily  described  by  set  theoretical  constructs.  This  is 
especially  true  in  fractal  functions  when  the  domain  is  and  the  object  created 
(set,  co-domain)  is  a  nonregular  shape.  This  is  one  reason  why  the  computer 
graphics  system  is  useful  in  the  investigation  of  fractal  functions.  The  computer 
can  model  the  infinite  function  and  display  a  finlce  approximation  of  the  created 
fractal  set. 

5.  Useful  Functional  Concepts 

It  is  often  helpful  to  clearly  understand  the  universe  of  discourse  within 
which  a  function  exists.  The  function  can  be  rigorously  defined  within  the  above 
constructs  but  lack  intuitive  appeal  due  to  its  complexity.  Mathematicians  have 
defined  many  useful  concepts  to  describe  functions.  The  concepts  applicable  to 
this  study  are  described  below. 

a.  Partial  Functions 

Most  of  the  fractal  functions  in  this  study  have  as  their  domain  some 
undefined  subset  of  R^.  It  is  useful  then  to  consider  them  as  partial  functions 
and  not  concern  ourselves  with  a  rigorous  description  of  the  domain  of  the 
fractal.  We  take  our  definition  from  [Ref.  2;pp.  201-202]. 

It  is  often  convenient  to  consider  a  function  from  a  subset  d*of  4  to  a  set  B 
without  exactly  specifying  the  domain  A  ’of  the  function.  Alternatively,  we  can 
view  such  a  situation  as  one  where  a  function  has  a  domain  A  and  codomain  R, 
but  the  value  of  the  function  does  not  exist  for  some  arguments  of  A.  This  is 
called  a  partiml  function. 

Definition: 

Let  A  and  B  be  sets.  A  partial  function  f  with  domain  A  and  codomain  B  is  any 
function  from  A’  to  B,  where  A*  is  a  subset  of  A.  For  any  s  which  is  an  element 
of  A  A’,  the  value  of  /  (z)  is  said  to  be  undefined. 

b.  Bijectivity 

It  is  often  useful  to  know  to  what  extent  a  function  maps  from  the 
uomain  to  the  codomain.  If  a  function  is  not  a  partial  function  and  every  point  of 
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the  domain  maps  to  a  point  in  the  codomain  then  we  want  to  know  if  all  points 
of  the  domain  A  in  the  mapping  /  (A)  map  to  distinct  points  in  the  codomain  B. 
We  may  also  want  to  know  to  what  extent  the  mapping  /(A)  covers  the  set  B. 
The  deHnition  of  b^jective,  surjective  and  mjective  functions  is  from  (Ref  2:pp. 
204]. 

Defimtiim: 

Let  /  be  a  function  /:  A-»B. 

(a)  /  is  surjective  (onto)  if  /(A)  =  B, 

(b)  y  is  injective  (one>to>one)  if  a  #  a'  implies  /(a)  #  /(a'), 

(c)  /  is  hijective  (one-to-one  and  onto)  if  /  b  both  surjective  ari  injective. 

6.  Functions  From  -♦ 

A  point  in  R^  space  is  specified  by  an  n-iuple  of  the  form 

(xi,22,X3, . ).  To  completely  specify  a  function  from  -»  R^  each  point  in 

th<  d-''':ain  must  map  to  a  point  in  the  codomain.  An  example  b: 

/  ;R*-»R2 

/  ((*1.  *2))  »  (*1*.»2*) 

This  function  is  well  defmed.  For  each  point  in  the  domain  of  the  function  we 
have  specified  a  unique  point  in  the  codomain. 

Most  of  the  functions  that  are  covered  in  thb  study  are  mappings  within 

or  Fractal  sets  exbt  in  all  finite  dimensions  but  it  b  impractical  at  this 
point  to  use  fractal  functions  beyond  the  fourth  dimension  in  view  of  the  graphics 
dbplay  medium’s  limitation  to  two  dimensions.  The  use  of  fractal  functions 
whose  dimension  b  between  3  and  4  b  currently  being  investigated  by  allowing 
the  function  to  roam  the  fourth  dimension  and  then  taking  time  slices  which 
yield  three  dimension  approximations  of  the  set  (Ref.  4). 

7.  Inductive  Definitions  of  Sets 

Functional  construcU  do  not  always  provide  a  convenient  means  of 
charactering  an  infinite  set.  It  b  sometimes  more  eloquent  and  powerful  to  use 
the  inductive  method  to  characterize  a  set. 
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Our  definition  of  inductive  definitions  of  sets  is  from  (Ref.  2:pp.  109-201]. 

An  inductive^  definition  of  a  set  always  consists  of  three  distinct  components. 

1.  The  6a«M.  or  basic  clause,  of  the,  definition  establishes  that  certain  objiticti. 
are  in  the  set.  The  basic  clause  establishes  that  a  set  is  not  empty  and  charac¬ 
terizes  the  "building  blocks"  (the  seeds  of  the  induction)  which  are  used  to  con¬ 
struct  the  set  from  the  inductive  clause. 

2.  The  mduetton.  or  ind»elit<e  clause,  of  an  inductive  definition  establishes 
the  ways  in  which  elements  of  the  set  can  be  combined  to  obtmn  new  elements. 

The  inductive  clause  always  asserts  that  if  objects  z,y,...,z  are  elements  of  the 
set,  then  they  cim  be  combined  in  certun  specific  ways  to  create  othe,r  elements 
of  the  set  (thus  from  the  basic  clause  (or  seeds)  of  the  induction  we  induce  the 
remaining  elements  of  the  set). 

3.  The  extremal  c/ause,  asserts  that  unleu  an  object  can  be  shown  to  be  a 
member  of  a  set  by  applying  the  basis  and  inductive  clauses  a  finite  number  of 
times,  then  the  object  is  not  a  member  of  the  set. 

An  evemplc  of  an  inductively  defined  sei  is: 

fBasisJ 

Ot  A 

(Induction) 

If  n  <  A,  then  (n  -1-2)  t  A 
(Extremal) 

No  integer  is  an  element  of  A  unless  it  can  be  shown  to  be  so  in  a  finite  cumber 
of  applications  of  clauses  1  and  2  above. 

The  set  that  we  defined  is  the  set  of  all  even  nonnegative  integers. 

S.  The  Path  To  Fractals 

The  path  to  fractals  by  the  non-mathematician  is  not  through  theory  but 
through  the  investigation  of  their  functions  and  methods  of  construction.  This 
investigation  (and  experimentation)  yields  considerable  insight  into  the  nature  of 
fractals. 

The  choice  of  which  set-descriptive  methodology  to  use  (functional  or 
inductive)  in  describing  a  set  is  often  a  matter  of  style  but  can  be  dictated  by 
necessity  if  one  method  is  inordinately  tedious. 

Most  of  the  frac.al  functions  that  are  introduced  in  this  study  use  the 
inductive  method  as  the  primary  functional  tool.  In  fact,  these  functions  arc  a 
hybrid  of  the  functional  and  inductive  constructs  described  above. 


*  Oflea  called  a  recurrcRce  defiaitba. 


14 


B.  DIMENSION 

The  classification  of  fractal  sets  from  non>fractal  sets  is  based  on  the 
dimensional  qualities  of  the  set.  To  understand  fractals  you  must  have  an 
appreciation  for  these  differences. 

The  concept  of  dimension  is  one  rife  with  difficulties.  Many  of  the  great 
mathematicians  have  attempted  to  define  dimension  as  a  rigorous  concept 
consistent  with  the  known  mathematical  systems.  For  each  of  their  attempts 
however,  the  concept  becomes  more  prone  to  contradiction  and  paradoxes. 

There  currently  exist  five  definitiona  of  dimension  that  date  back  to  the  late 
1800’s^.  The  classiHcation  of  Fractal  sets  into  a  class  of  sets  is  the  result  of  the 
discovery  of  functions  that  created  sets  which  did  not  fit  comfortably  into  the 
topological  definition  of  dimension  (which  was  the  accepted  definition  at  the 
ti^^e).  Fractal  sets  are  rigorously  classified  as  those  sets  that  demonstrate  a 
difference  between  the  Hausdorff-Besicovitch  dimension  and  the  standard 
topological  dimension. 

1.  An  In.uitive  Approach  to  Dimension 

Dimension  is  a  concept  that  seems  intuitive  when  it  is  first  introduced  in 
Euclidean  geometrv  a^  the  standard  three  dimensions.  Long  after  Euclid  made 
the  first  attempts  at  defining  dimension  and  concurrent  with  the  discovery  of 
atomic  particle  pbysi:s  the  concept  of  dimension  was  rethought  by  the  prominent 
mathematicians  of  the  wime  Ihis  was  necessary  to  realign  the  mathematical 
model  for  the  geometry  of  objects  with  the  new  view  of  what  those  objects  were 
made  of.  Our  increasing  ability  to  focus  on  tne  nature  of  matter  inevitably 
causes  the  models  we  use  to  change. 

The  dilemma  that  arose  from  the  new  concepts  of  dimension  quickly 
developed  into  a  theoretu  d  debate  that  left  intuition  behind.  When  human 
Intuition  fails,  wc  rausi  rely  upon  well  founded  modeb  that  are  based  on  axioim 
of  basic  mathematical  truth.  It  is  only  through  the  rigorous  investigation  of  our 
mathematical  models  that  allows  us  to  go  beyond  ir*uition  and  investigate  the 

*  1).  CseUH’  «nd  Minkowaki;  2).  Bouligud  and  Minkowiki;  3).  Pontrajgin,  Schnirelman  and 
Kotomogorov,  Tihomirov;  4).  HausdorlT>H«aicovUcb  and  S).  the  topological  dimension  (there  are 
cHhera).  Moat  of  these  definitiotu  etc  concerned  with  the  most  efTicient  method  of  covering  a  set 
(i.e.  are  topological  coneeima|. 
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true  physical  nature  of  the  objects  we  model.  The  debate  still  rages  today  and 
borders  on  the  philosophical.  Two  examples  should  suffice  to  demonstrate  the 
complexity  and  possible  paradoxes  that  can  arise  from  dimension  theory. 

The  first  is  from  [Ref.  5:pp.  323-344]. 


Consider  the  ww  in  which  we  deflne  the  density  of  air  at  a  given  point  and  at  a 

given  moment.  We  picture  a  sphere  of  volume  V  centered, at  that  point  and  in- 
uding  the  mass  M.  The  quotient  M/V  is  the  mean  density  within  the  sphere, 
and  by  the  true  density  we  denote  some  limiting  value  of  this  quotient.  This  no¬ 
tion,  however,  implies  that  at  the, given  moment  the  mean  density  is  practically 
constant  for  spheres  below  a  certain  volume.  This  mean  density  may  be  notably 
different  for  spheres  containing  1,000  cubic  meters  and  1  cubic  centimeter 
respectively. 


Supimse  the  volume  becomes  continually  smaller.  Instead  of  becoming  less  and 
less  important,  these  fluctuations  come  to  increase.  For  scales  at  which  the 
brownian  motion  shows  great  activity,  fluctuations  may  attain  1  part  in  1,C90, 
and  they  become  of  the  order  of  1  part  in  5  when  the  radius  of  the  hypothetical 
spherule  becomes  of  the  order  of  a  Hundredth  of  a  micron. 


One  step  further  and  our  spherule  becomes  of  the  order  of  a  molecule  radius.  In 
a  gas  it  will  generally  lie  in  intcrmolecular  spue,  where  its  mean  density  will 
henceforth  vaniak.  At  our  point  the  true  density  will  also  vanish.  But  about 
once  in  a  thousand  times  that  point  will  lie  within  a  molecule,  and  the  mean 
density  will  be  a  thousand  times  higher  than  the  value  we  usually  take  to  be 
the  true  density  of  the  gas. 


Let  our  spherule  grow  steadily  smaller.  Soon,  except  under  ,exceptionid  cir¬ 
cumstances,  it  will  necome  empty  and  remain  so  henceforth  owing  to  the  intra- 
atomic  space;  the  true  density  vanishes  almost  everywhere,  except  at  an  infinite 
number  of  isolated  points,  where  it  reaches  an  infinite  value. 


The  second  is  from  [Ref.  l:pp.  17-18]. 

Consider  a  ball  of  10  cm  diameter  made,  of  a  thick  thread  of  1  mm  diameter  that 
(in  latent  fashion)  possesses  several  distinct  effective  dimensions. 

To  an  observer  placed  far  away,  the  ball  appears  as  a  zero-dimensional  figure:  a 
point.  As  seen  from  a  distance  of  10  cm  resolution,  the  ball  of  thread  is  a  three- 
dimensional  figure.  At  to  mm.  it  is  a  mess  of  one-dimensional  threads.  At  0.1 
mm,  each  thread  becomes  a  column  and  the  whole  becomes  a  thre^dtmensional 
figure  again.  At  0.01  mm,  each  column  dissolves  into  fibers,  and  the  ball  again 
b^omes  one-dimensional,  and  so  on,  with  the  dimension  crossing  over  repeated¬ 
ly  from  one  value  to  another.  When  the  ball  is  represented  by  a  huite  number  of 
atomlike  pinpoints,  it  becomes  zerodimensionaJ  again. 

It  b  interesting  to  note  that  each  of  these  examples  demonstrate 
dimension  as  a  reflection  of  physical  properties  dependent  on  the  observers  point 
of  reference.  Each  ends  with  reference  to  the  paradox  of  atomic  particles.  That 
paradox  is,  for  any  collection  of  finite  (or  countably  infinite)  points,  the 
dimension  b  zero  [Ref.  6:pp.  1-8].  Since  the  earth  and  sun  each  have  a  finite 
collection  of  atoms  then  accordingly  their  dimension  b  zero.  Dimension  cxbts 
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only  for  a  mathematical  continuum  and  as  such  lacks  application  to  the  physical 
universe  as  we  currently  know  it^. 

The  two  properties  (continuity  and  dimension)  cannot  be  separated. 
Before  the  advent  of  atomic  theory,  matter  was  viewed  as  continuous  and 
composed  of  basic  elements  that  were  indivisible.  While  the  debate  raged  over  the 
practical  and  philosophical  aspects  of  the  nature  of  matter  it  became  apparent 
that  the  mathematical  models  which  represent  matter  would  have  to  change.  It  is 
not  practical  to  represent  objects  by  representing  each  atom  and  its  position 
relative  to  the  entire  set.  The  power  of  modeling  would  thus  be  lost;  that  is,  the 
ability  to  model  complex  objects  and  their  interaction  by  relatively  simple 
constructs.  Thus  the  fact  is  reinforced  that  models  can  only  represent  objects 
through  gross  approximations  and  that  the  model  is  only  effective  for  a  restricted 
frame  of  reference.  Without  thb  realization,  dimension  would  have  very  little 
application. 

2.  Topological  Dimension 
a.  An  Intuitive  Approach 

The  concept  of  dimension  is  very  old.  It  is  based  on  the  algebraic 
concepts  of  Euclidean  n  space  and  the  notion  that  a  set  has  dimension  n  if  the 
least  number  of  real  parameters  needed  to  describe  its  points  was  n.  This  fuzzy 
definition  was  accepted  for  a  very  long  time  until  the  advent  of  Cantor  and  set 
theory.  Cantor  showed  that  dimension  can  be  changed  by  a  1^1  transformation 
from  an  interval  to  a  planar  object.  The  fuzzy  notion  of  dimension,  as  defined, 
was  challenged  and  required  rethinking. 

The  mathematicians  who  did  not  accept  many  of  the  findings  of  set 
theory  at  the  time  (but  who  could  not  disregard  Cantors  findings)  began  to 
consider  ways  of  explicitly  defining  dimension.  The  new  definition  would  have  to 
be  applicable  to  the  6i2orre  functions  of  Cantor,  Koch  and  Peano  as  well  as  the 

*  The  Ml  ihcorelieel  concepl*  of  linileneM,  countably  infinite  and  uncounlably  infinite  (cor* 
tinuous)  eeti  tarty  with  them  very  profound  implicaliont.  It  it  premature  to  view  matter  aj  mere¬ 
ly  collections  of  finite  atoms.  Science  may  yet  find  true  continuity  (in  the  mathematical  sense)  in 
atomic  matter  and  the  universe.  For  now,  matter  is  what  it  U  and  our  proauncUtixMit  upon  it  will 
not  ckaa|e  its  true  texture. 
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relatively  simple  objects  that  had  previously  fit  into  the  old  definition  without 
contradiction. 

There  was  one  crucial  problem  that  Cantor’s  findings  raised 
(Ref.  6:pp.  4-6]: 

An  extremely  important  question  was  left  open:  Is  it  possible  to  establish  a 
correspondmce  Mtween  Euclidean  n  space  ancf  Euclidean  m  space  combiniUe  the 
feature  of  both  Cantor’s  and  IVano’s  constructions,  j.e.  a  correspondence  v^ich 
IS  both  1:1  ^  and  continuous?  The  question  is  crucial  since  the  existence  of  a 
transformation  of  the  stated  type  between  Euclidean  n  space  and  Euclidean  m 
space  would  signify  that  dimension  (in  the  natural  sense  that  Euclidean  n  space 
has  dimension  n  )  has  no  topological  meaning  whatsoever. 

This  fundamental  problem  was  answered  in  1011  by  Brouwer.  He 
proved  that  Euclidean  n  space  and  Euclidean  m  space  were  not  homeomorphic 
unless  n  equals  m.  To  say  that  two  spaces  A  and  B  are  homeomorphic  means 
that  a  mapping  f  :A  ~*B  exists,  such  that  /  is  continuous  over  A  and  b(jective. 
Additionally,  the  inverse  of  this  mapping  /  ,  is  continuous  over  B  and 

bijective.  If  two  spaces  are  homeomorphic  then  it  is  analogous  to  saying  that  they 
are  topologically  eguivalerU. 

Further  research  was  done  and  a  precise  definition  of  topological 
dimension  of  a  set  was  developed.  This  definition  assigned  an  integer  value  as  the 
dimension  of  any  set  based  on  its  topological  properties, 
b.  Definition  of  Topological  Dimension  n 

The  rigorous  investigation  of  dimension  is  beyond  the  scope  of  this 
study  but  the  following  definition  is  included  for  completeness  [Ref.  6:pp.  24): 

Rouchty  speaking,  we  may  say  that  a  space  has  dimension  $  n  if  an  arbitrarily 
small  piece  of  the  spwe  surrounding  each  .poiijt  may  be  delimited  by  subsets  of 
dimension  £  n  -  1.  This  method  of  definition  is  inductive,  and  an  elegant  start- 
trig  point  for  the  induction  b  given  by  prescribing  the  null  set  as  the  (-1)  dimen¬ 
sional  space. 

Dtfimtiou: 


The  empty  set  and  only  the  empty  set  has  dimenaioit  -1. 

A  space  X  has  dimension  «  n  (  n  ^  0)  at  o  point  p  if  p  has  arbilrarily  small 
neighborhoods  whose  boundaries  have  dimension  ^  n  -  1. 

X  has  dimenaton  4  n,  dim  A  <  n,  if  A  has  dimension  ^  n  at  each  of  its  points. 

A  has  dimension  n  at  point  p  if  it  is  true  that  A  lias  dimension  <  n  at  p  and  it  b 

false  that  A  has  dimension  4  n  -  1  at  p. 

A'  hn$  dimension  a  if  dim  A  ^  n  b  true  and  dim  A*  ^  n  -  1  b  false. 

A  has  dimenaibii  oe  if  dim  A  ^  n  b  false  for  each  a. 
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The  topological  dimension  is  rigorous  and  consistent  for  all  sets  that 
eirist  within  a  mdrte  space.  The  problem  that  arises  with  fractal  sets  and  its 
topological  dimension  is  not  that  the  topological  dimension  is  wrong.  Fractal  sets 
like  all  sets  in  a  metric  space  exhibit  a  topological  dimension.  The  question  is 
then,  is  the  topological  dimension  an  accurate  description  of  the  dimension  of  the 
set  or  can  we  find  a  better  way  to  characterize  the  dimensional  qualities  of  the 
set?  This  question  can  be  extended;  is  the  topological  definition  of  dimension 
useful  and  consistent  with  the  of  dimension  and  space?  Can  we  devise  a 

better  concept  which  can  further  refine  dimension  and  make  it  more  useful^. 

3.  The  Hausdofff-Besicovitch  Dimension 

This  section  is  intentionally  brief  due  to  the  subject’s  complexity  and  to 
the  lack  of  practical  technique  that  it  yields.  The  Hausdorff  measure  of  a  set  is  a 
complex  characterization  of  a  method  for  covering  a  set.  HausdorfTs  theorem  is 
proved  using  the  existetUial  qualities  of  infinite  sets  in  a  metric  space  0.  While 
the  theorem  may  be  important  to  mathematical  theory,  it  proves  unfortunate 
that  there  is  no  straightforward  practical  method  for  determining  the  Hausdorff 
measure  of  a  set. 

a.  An  Intuitive  Approach  to  the  Hausdorff  Dimension 

The  acceptance  of  the  Hausdorff  method  for  covering  a  set  as  a 
measure  of  dimension  is  not  universal  (Ref.  6:pp.  102]  and  (Ref.  l:pp.  363*365]. 
The  debate  is  between  the  dueipiinca  of  topology  and  metrics  and  is  not  wholly 
germane  to  this  study.  It  is  beneficial  to  divorce  ourselves  from  the  debate  and 
consider  both  the  topological  dimension  Df>  and  the  Hausdorif-Besicovitch 
dimension  (HB)  as  merely  measures  of  iptaUtiat  of  a  set’s  structure.  Certainly 
sets  exist  that  have  a  topological  dimension  equal  to  1  but  in  no  way  resemble  a 
simple  Euclidean  curve  If  the  Hausdorff  dimension  yields  a  better  measure  of  a 
set’s  structure  that  provides  a  mathematical  and  intuitive  difference  that  is  useful 
to  us,  it  would  be  beneficial  for  us  to  investigate  it. 


*  Try  not  to  aKrtbr  § randtOM  implicalionr  to  a  set's  dimension  (the  fourth  dimension  as  time 
or  some  such)  as  this  is  premature  at  best.  Rather,  view  a  set's  dimension  as  merely  descriptive 
teritiinolofy  much  tike  the  terminology  of  biitcliviiy  describes  a  function'a  characleritlics.  The 
problem  most  people  have  with  this  mental  abstraction  is  the  visual  reinforcement  that  they  re¬ 
ceived  from  the  notion  of  the  atandaid  three  dimensions. 
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The  HausdorfT  measure  of  a  set  was  developed  during  the  same 
period  that  the  new  topological  dimension  was  invented  to  solve  the  paradoxes  of 
Cantor.  The  topological  dimension  was  based  on  the  idea  of  a  neighborhood  of  a 
point  within  a  Euclidean  space  of  R^.  The  connection  to  metric  spaces  and  the 
idea  of  measure  is  obvious  when  you  consider  that  the  HausdorfT  measure  of  a  set 
is  also  based  on  this  notion  of  a  spherical  neighborhood  and  what  HausdorfT  calls 
the  test  fitneium  of  a  set.  The  test  function  of  a  set  denoted  h(p )  is  a  function 
that  characterizes  the  **iest**  method  of  covering  a  point  with  the  spherical  ball 
of  radius  p  that  covers  points  of  the  set,  which  in  their  union,  cover  the  entire 
set. 

Consider  for  example  the  test  function  for  a  surface  within  R^.  A 
surface  can  be  covered  by  discs  (circles).  The  formula  for  the  area  of  a  circle 
becomes  the  test  function  for  the  surface.  The  formula  for  the  area  of  a  circle 
always  contains  the  constant  factor  ir  multiplied  by  the  square  of  the  radius  r. 
This  radius  is  the  measure  p  as  above.  This  leaves  us  with  a  test  function  for  a 
planar  shape  in  R*  of  h(p )  =  irp*. 

You  might  expect  that  the  test  function  for  a  spherical  neighborhood 
in  a  Euclidean  space  above  R^  would  be  very  difficult  to  imagine,  and  indeed  it 
is.  HausdorfT  further  complicated  the  idea  of  test  functions  (even  within  the  lower 
dimensions)  by  allowing  a  test  function  to  assume  a  non-integer  parameter  d  so 
that  the  test  function  h(p )  =  1&(d)p  ^  could  have  a  real-valued  parameter  d. 
This  further  refinement  of  the  test  function  allowed  HausdorfT  to  mahe  assertions 
about  how  this  test  function  h(p )  behaved  when  the  parameters  p  and  d  were 
allowed  to  vary. 

HausdorfT  imagined  the  parameter  p  reducing  in  size  until  it 
approached  zero.  The  effect  of  this  on  our  disc  example  is  increasingly  smaller 
and  smaller  discs  around  points  of  the  planar  set.  As  the  disc  size  is  decreased, 
fewer  points  of  the  set  are  contained  in  each  disc  neighborhood.  This  requires 
more  discs  to  cover  the  set.  As  the  parameter  p  becomes  infinitely  small  the 
number  of  discs  required  to  cover  the  planar  set  approaches  oo.  \Ve  allow  the 
parameter  to  grow  arbitrarily  small.  It  is  interesting  to  study  the  test  funetbn 


and  8«e  what  happens  to  the  total  area  when  the  areas  of  the  collection  of  discs 
which  cover  the  planar  set  are  summed. 

Let’s  reflect  upon  the  mathematical  process  that  we  are  developing. 
When  we  attempt  to  approximate  the  area  of  the  planar  mi  by  the  union  of  the 
discs  which  make  up  its  cover,  we  are  essentially  observing  small  patehee  of  the 
surface  and  approximating  the  area  of  the  net  by  making  assertions  about  the 
intrinsic  qualities  of  the  patch.  The  notion  that  this  measure  is  merely  an 
approximation  is  important.  As  the  size  of  our  patch  grows  increasingly  small,  we 
can  expect  that  we  will  get  a  better  fit  with  our  patchy  and  hence  a  better 
approximation  of  the  area.  The  notions  of  approximation  and  fit  become 
increasingly  helpful  when  you  consider  functions  which  describe  sets  of  infinitely 
rough  texture  as  we  And  in  fractal  sets. 

The  importance  of  Hausdorffs  discovery  lies  in  the  fact  that  for  a 
test  function  of  a  set,  the  parameter  d  is  epeeiat.  As  p-*0  he  discovered  that 
there  existed  a  unique  real  number  d  such  that  for  d'<d,  the  inflnum^  defmed 
by  the  test  function  using  d^  approaches  oq  (for  any  countably  infinite  set).  And 
for  a  d'>d  the  corresponding  inflnum  approaches  zero.  This  means  that 
set  has  a  parameter  which  can  be  associated  with  it  that  is  closely  relatet^ 
amount  of  epaee  that  it  occupies.  This  number  is  the  Hausdorff-Besicovi>'x:h 
dimension  (measure)  of  the  set. 

These  results  only  tell  us  that  a  number  and  function  exist.  They  do 
not  tell  us  how  to  compute  them  in  the  general  case.  This  gives  us  the  quandary 
of  dealing  with  the  HB  dimension  as  a  known  concept  that  we  can  make 
allusions  to,  but  can  rarely  compute  (at  least  at  the  present  time), 
b.  Definition  of  the  HausdorflT  Dimension 

The  formal  defmition  of  the  HausdorfT  dimension  requires  that  we 
formalize  Uie  intuitive  discussion  above.  We  first  define  what  a  p>measure  of  a  set 
is  (analogous  to  the  disc  above)  (Ref.  6:pp.  102' 103]. 

fip•(iimen|loniJ  measure  for  each  non-negative  real  number  p  was  deftned  by 
ausdorfi  for  arbitrary  metric  snacea.  This  measure  is  a  melrtcal  concept,  white 
dinwnslon  is  purely  (opdf(WC«f.  Nevertheless  there  is  a  strong  connection 
between  the  two  coneepi;^.  for  it  turns  out  that  a  space  of  (tppotofieaJn  dimen¬ 
sion  R  must  have  positive  R-ditneasiooaJ  n.easure  (lac  Hau*dorffme*»ur^, 

*  For  our  example  this  would  be  the  aura  of  the  area*  of  the  discs. 


Defbutiom 

Let  X  be  a  space  and  p  an  arbitrary  real  number,  0  ^  p  <oo.  Given  c  >0  let 

m%  =  infjj 
<•1 

where  X  =  is  any  decomuosition  of  X  into  a  countable 

number  of  subsets  of  diameter  less  than  e,  and  the  superscript  p  denotes  ex¬ 
ponentiation.  Let 


nif  (X)  =  luj^  m'p  (X). 

(X)  is  coiled  the  p>measure  or  (p^imenaional)  measure  of  X. 

Propoaitwn: 

If  p  <f  then  (X)  ^  m,  (X);  in  fact  p<q  and  (X)  <oo  imply  m,  (X)  =  0. 

Conversely,  if  p  ^  then  (X)  ^  (X);  and  if  p  >q  and 

(X)  <00  inpliet  iq,  (X!)  =  oo. 


c.  Mandelbrot’s  Misgivings  about  the  HB  Measure 

It  is  clear  from  the  previous  definition  that  for  practical  applications 
the  Hausdorff  dimension  is  difficult  to  compute  directly.  In  (Ref.  l:pp.  14-10]  and 
(Ref  7|,  Mandelbrot  expressed  a  distaste  for  the  focusing  of  attention  on  the  HD 
dimension.  He  states: 

developed  the  derinition  of  fractals  using  the  tomlogic^  and  Hausdorff  di¬ 
mensions  in  response  to  colleagues  who  urged  me  to  do  so.  They  felt  that  it  was 
necessary  to  rigorously  define  the  concept  within  firm  mathematical  criterja.  f 
have  come  to  believe  that  an  empirical  definition  would  be  more  beneficial  at 
this  time  because  the  present  definition  denies  the  inclusion  of  some  shapes  that 
could  best  be  describe  as  fractals.** 

Mandelbrot  believes  that  the  definition  of  the  Hausdorff  dimension  is 
too  diffieutt  to  deal  with  and  is  perhaps  too  restrictive.  He  prefers  to  focus  on  the 
befuivioral  aspects  of  the  recurrence  relationship  involved  in  fractal  and  the 
empirical  results  from  these  equations. 

A  practical  application  of  fractal  functions  in  computer  graphics 
does,  by  necessity,  bend  to  this  same  paradigm.  This  realixation  should  not  blind 
us  to  the  fundamental  nature  of  a  fractal  equation's  uniqueness,  however.  It  is 
important  to  understand  the  dimensional  aspect  of  the  fractal  discourse  to 
appreciate  the  poientiat  importance  of  fractals. 

It  is  not  preordained  that  fractal  equations  model  nature  with  a 
greater  degree  of  accuracy  then  does  the  Euclidean  model.  The  future  may  prove 
the  fractal  model  the  superior  method,  however.  The  dimensiotud  qualities  of 
fractal  functions  may  be  the  aspect  that  proves  this  to  be  so. 
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4.  Why  Ccrjider  Dimension 

For  the  purpose  of  this  study,  it  is  not  important  that  a  rigorous  feel  for 
the  mathematical  properties  of  dimension  theory  be  grasped.  In  fact  one  needs  no 
knowledge  of  dimension  to  use  fractal  tecl‘':’que8  in  the  generation  ot  computer 
graphics  terrain.  The  literature  is  rife  \dth  articles  about  fractal  objects  in 
computer  graphics  and  it  seems  de  rigue’jr  to  include  an  approximate  fractal 
dimension  as  pert  of  its  description.  The  techniques  used  to  approximate 
dimension  as  presented  in  [Ref.  l:pp.  5&-ii7]  are  mathematically  unproven^.  More 
importantly,  the  dimension  yields  littla  intuitive  insight;  one  is  hard  pressed  to 
describe  the  differences  between  an  object  with  an  approximate  dimension  of  1^.37 
and  another  with  a  dimension  of  2.45.  The  pictures  are  much  more  descriptive. 

One  use  of  approximate  fractal  dimension  is  to  describe  an  object’s 
relative  roughness.  It  is  beneficial  to  view  a  fractal  dimension  as  degrees  of 
roughness  between  the  standard  three  dimensions.  If  the  dimension  is  between  1 
and  2  then  the  object  should  be  a  vesy  irregular  curve.  If  the  dimension  of  that 
curve  approaches  1  then  the  curve  is  probably  not  very  rough  and  would  lack  any 
interesting  diversion  from  an  ordinary  plane  curve.  If  the  dimeniuon  of  the  curve 
approaches  2  then  the  curve  becomes  like  plane  or  filled  polygon  and  again 
lacks  appeal.  The  most  interesting  fractal  curves  are  those  which  demonstrate 
dimensions  nearer  the  center  of  the  scale  between  the  standrrd  Euclidean 
dimensions.  A  similar  argument  can  be  made  for  soiid  objects  with  dimensions 
between  2  and  3. 

This  is  not  to  say  that  fractal  geometry  is  not  a  powerful  tool  for  the 
graphics  programmer.  The  evidence  of  the  power  of  fractals  to  model  objects  of 
considerable  complexity  is  clearly  demonstrated.  To  date,  this  power  has  not 
been  matched  by  other  standard  methods. 

The  graphics  programmer  should  not  concern  himself  greatly  with  the 
dimension  tliai  is  demonstrated  at  difrcieni  leveb  of  object  construction.  He  roust 
concern  himself  with  the  techniques  of  construction  and  the  realism  that  is 
achieved. 

^  Msiuielbrot's  meUiod  for  the  Hiiutdi^ff-lhiKovkek  dimctistoa  for  Ron^-msdoni 

•ett  buitl  Uuovfti  self  simiUr  will  be  laUoditced  ia  secUoa  C  *fUr  Um  Hoek  twvc  is 

described. 
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C.  FRACTAL  CURVES  AND  SETS 

1.  Derinition  of  Fractal  Seti 

We  t«ke  our  definition  of  fractal  seta  from  (Ref  l:Chap  3  and  pp.  361]. 

Definitiom: 

A  frmdal  set  is  a  set  for  which  the  Bmudar^Buieornttk  dimension  strictly 

excecos  the  tofolopeal  dimension. 

As  we  have  established  and  is  emphasized  by  Mandelbrot,  this  definition 
is  not  very  useful.  The  definitions  of  topological  and  HausdorfT  dimensions  are 
very  involved.  It  is  a  gargantuan  effort  to  prove  that  a  set  has  a  topological  or 
Hausdoiff  dimension  (if  one  desires  complete  rigor,  typically  the  topological 
dimension  is  derived  by  the  least  parameter  approach  (section  2)).  When  using 
fractal  functions  then,  it  is  practical  to  asasMie  that  because  the  functional 
method  you  use  is  Jraetal'tike  that  the  dimension  is  fractal. 

The  functional  techniques  to  be  introduced  have  a  certain  methodology 
that  creates  fractal  sets  with  a  behsvior  that  is  disciplined  and  predictable.  The 
assumption  b,  since  these  methodologies  are  well  behaved,  that  any  set  created 
by  these  methods  (with  some  careful  restrictions)  will  itself  be  a  fractal  set. 

We  are  left  with  a  practical  methodology  whereby  we  discover  fractal 
functional  methods,  prove  that  the  set  created  b  fractal,  characterize  the  fractal 
part  of  the  functional  method  (carefully)  and  then  ewitriiia  that  method  as  a 
fractal  method.  If  <mc’§  purpose  b  a  practical  application  of  fractal  techniques* 
and  not  a  rigcurous  mathematical  investigation  then  thb  b  a  reasonable  and 
practical  approach.  Thb  approach  b  taleen  in  the  remainder  of  thb  study. 

2.  Constructing  Fractal  Sets 

In  order  to  describe  the  construction  of  the  Koch  curve,  it  b  necessary  to 
present  terminology  introduced  by  Mandelbrot  |Bef.  l:pp.  34*35). 

We  use  a  geometrical  shape  (at  first  a  straight  tine)  and  call  Uib  shape 
an  iNiTiATOSL  We  create  another  shape  that  b  constructed  with  shapes 
similar  to  the  initiator  and  call  thb  a  GENEBA  TOB  We  d^ne  a  sequence  of 

*  A  woriciai  model  or  rqnsltoa  wber*  you  Ut  oaly  eoMeriMd  abottl  behavior  aad  sol  the 
exact  matltematkal  properik*. 
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traiuforinations  upon  all  current  initiators  (suppose  there  are  m  such  initiators) 
in  the  construction  by  applying  the  generator  to  all  initiators.  This  creates  a  new 
construction  that  consists  of  m  x  r  (where  r  is  the  number  of  distinct  parts  of  the 
generator)  sides  where  each  side  is  a  shape  that  is  similar  to  the  initiator.  The 
next  step  is  to  again  apply  the  generator  to  all  initiators. 

This  recursive  defuiition  has  no  terminating  event  but  is  continued  ad 
infim'tum.  Thb  functional  process  is  well  suited  for  recursion  because  the 
application  of  the  generator  to  the  initiator  is  constant  with  respect  to  method 
and  varies  only  to  scale.  It  is  also  well  suited  for  parallel  processing  (in  the 
computer  science  sensc']  because  each  application  of  a  generator  on  all  current 
initiators  is  independent. 

These  concepts  are  probably  confusing  at  this  point  and  were  especially 
difficult  to  visualize  when  they  were  first  envisioned  because  the  authors  had  few 
toob  beyond  mental  imagery  to  convey  their  point.  This  b  probably  why  they 
were  largely  ignored  for  70  years.  It  b  much  easier  to  vbualize  these  functions 
when  they  are  shown  on  a  computer  graphics  dbplay. 

3.  The  Koch  curve. 

The  mathematicians  Koch  and  Cantor  developed  functions  which 
attempted  to  challenge  the  mathematical  modeb  of  continuity  and 
differentiability.  These  equations  were  developed  during  the  great  debate  on  set 
theory  and  were  used  by  Cantor  in  arguing  for  bb  theory.  These  funetioos  were 
Uhe  none  before,  using  constructions  which  played  upon  natural  geometric 
constructs  but  when  combined  with  the  power  of  infinite  recursion  became  sets 
which  defied  intuition.  It  was  not  until  much  later  that  mathematicians  were  able 
to  reconcile  these  functions  with  algebra  and  set  theory.  The  function  to  be 
introduced  has  been  proven  to  have  a  Hausdorff-BesieoiMtek  dimension  which 
creeds  its  topological  dimension^. 

The  Koch  cur.'c  is  a  very  beautiful  curve  that  at  first  gives  the  observer 
the  impression  of  a  snowflaiie  or  a  coastline  (Fig  2.1).  Mandelbrot  uses  thb  type 
of  construction  (with  variation  — >  i.e.  randomness  or  less  behaved  generators)  to 
draw  coastlines  that  look  very  realistic.  To  understand  thb  construction  b  to 

*  Tkttt  hrseuit  dd  clutt  sad  il  »  potibb  lo  ri(ore«tIr  ptove  u  iS  be  so. 
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understand  one  method  of  obtaining  a  fractal  set  from  a  well  defined  non-fractal 
set  (R^)  using  non-random  techniques.  This  method  of  construction  (in  the 
general  sense)  is  very  powerful  and  is  used  throughout  this  thesis. 

To  construct  the  Koch  curve,  we  use  three  initiators  (line  segments)  of 
equal  length  and  join  them  to  form  an  equilateral  triangle^*^.  To  construct  the 
generator  we  use  four  line  segments  that  are  each  the  length  of  the  initiators 

and  apply  these  to  each  initiator  (Fig  2.1).  This  yields  a  new  geometric  figure 
with  12  sides  versus  the  original  3  and  a  total  perimeter  length  of  4  units  of 
length  versus  the  original  3  units. 

Figure  2.1  demonstrates  the  first  and  second  recursive  iterations  of 
building  the  Koch  curve.  Observe  how  the  progression  develops  to  yield  the  final 
figure  in  Figure  2.1.  Imagine  this  progression  occurring  indefinitely. 

-  With  each  iteration  of  applying  the  generators  the  total  perimeter  length 

1 

increases  by  —  over  the  previous  perimeter. 

3 

-  The  lentfih  of  the  curve  begins  to  increase  without  bound  even  though  the 
length  of  the  initiator  decreases  to  an  infinitely  small  length.  Hence  the 
curve  s  length  is  unbounded  with  no  point  intersecting  but  yet  is  contained 
in  a  small  bounded  two  dimensional  area. 

-  The  points  of  the  curve  are  by  construction  only  the  end-points  of  each 
initiator  and  each  point  is  clearly  distinct  from  the  other  (  no  two  points  are 
connected) . 

-  Although  each  point  is  distinct  at  any  one  level  of  the  curve  construction,  it 
can  be  proven  that  the  curve  when  viewed  in  the  limit  is  continuous  at  every 
point. 

-  That  due  to  the  above  qualities  the  curve  is  not  differentiable  at  ANY  point. 

It  is  important  to  realize  that  the  endpoints  of  the  lines  (initiators)  are  the  only 
points  of  the  curve.  The  line  only  serves  as  a  vehicle  by  which  the  points  may  be 
easily  determined.  The  exact  same  set  could  be  built  using  ISO**  arcs  as  initiators. 

An  algorithm  and  computer  graphics  program  for  the  construction  of  the 
Koch  curve  is  presented  in  Chapter  4  and  Appendix  A. 


The  choice  of  an  equilateral  triangle  was  arbitrary.  We  could  have  chosen  any  shape  as 
long  as  it  wu  made  up  of  Initiators  a;(d  avoided  intersecting  lines  during  recursion. 


THE  KOCH  CURVE 
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Figure  2.1  The  Kcch  Curve. 
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4.  Mandelbrot’s  Dimension  Approximation  Function 

In  the  above  section,  one  functional  construction  technique  for  building  a 
fractal  set  has  been  introduced.  It  is  possible  to  approximate  the  dimension  of  a 
fractal  curve  that  is  built  using  these  constructs.  In  [Ref.  l:pp  56-57],  Mandelbrot 
introduces  a  function  that  is  based  on  the  similarity  properties  of  the  above 
technique.  This  function  has  a  real  exponent  D  that  is  t^^  approximate 
dimension  of  the  fractal  set^^. 

Consider  a  method  of  paving  (covering?)  a  Euclidean  shape.  Divide  a  line 

segment  into  N  segments  with  each  segment  a  part  of  the  original  segment  such 

that  the  sum  of  the  lengths  of  the  N  segments  equals  the  length  of  the  original 

segment.  It  follows  then,  that  the  sum  of  the  ratios  of  the  divided  segments 

,  ^  .  •  I  *  1  Divided  segment  a 

lengths  to  the  origmal  segment  length  e.g.  r.  =  ■. — ■ — -  must 

Qrtgtnal  aegment 

,  i 

equal  1. 

N 

We  know  that  the  dimenmon  of  a  line  is  equal  to  1.  If  we  raise  each  of  the  above 
ratios  to  the  power  D  (where  D  »  1)  the  equivalence  still  holds. 


Let's  allow  the  Koch  function  to  assume  a  similar  dimensional  relationship  but 
treat  D  as  a  real  valued  unknown.  Refer  once  again  to  Figure  2.1.  Notice  that 
the  length  of  the  four  line  segments  that  make  up  the  generator  have  a  length 

ratio  of  ~  to  the  initiator.  Call  this  ratio  !*„.  Notice  that  the  generator  is  made 
3 

up  of  four  initiator  shapes.  Call  thb  number  N. 


Substituting  and  Solving  for  D: 


\D 

N  X  |rn,j  =  1 


D  =  ^251  »  1.2618 

logs 


Which  is  equal  to  the  Hausdorff-Besicovitch  dimension  of  the  Koch  curve^^. 

This  is  not  a  proof  of  a  general  equation  for  the  fractal  (Hausdorff) 
dimension  of  a  self  similar  fractal  set  but  implies  that  a  general  dimension 
generating  function  is  possible: 

N  /  \D 

G(I>)  =  E  (r«  =  1 

Whtre  N  =  number  of  sides  of  the  nenerator. 

Where  Rq  =  ratio  of  tide  m  to  the  mitiotor. 

Mandelbrot  claims  that  experimental  evidence  suggests  that  this  equation  holds 
whenever  this  functional  method  is  used. 

When  each  segment  of  the  generator  is  a  fixed  ratio  to  the  initiator  (as  is 
the  case  with  the  Koch  curve)  then  the  solution  to  this  equation  is  trivial: 

D  ^  .  J2?IgL 

log(rr) 

Pia 

S*  Functional  Characterttation  of  the  Koch  Method 

A  complete  and  rigorous  inductive  defmition  of  most  fractal  functions 
can  stretch  the  notational  capabilities  of  the  symbolic  aspects  of  the  inductive 
and  functional  methods.  It  b  thus  generally  impractical  to  use  these  methods 


'*The  Koch  curve  wu  proven  to  have  a  Hausdorff-Beatcoviicb  dimension  equal 
»  1.2018  and  a  topolo§ical  dimension  equal  to  1,  jHausdorfl:  Dimension  und  ausseres  Mass}. 

logs 


except  in  a  verbose  and  non-rigorotw  manner.  It  can  be  insightful  however,  to 
dissect  the  beast  (once!)  and  hopefully  gain  further  intuitive  insight. 

To  simplify  the  process,  we  define  a  Koch  half-line  as  in  the  above  fractal 
but  with  a  single  line  segment  as  the  Initiator  (versus  the  equilateral  tciangk). 
We  consider  the  line  interval  of  [(0,0), (1,0)]  within  the  set  This  restricts  the 
fractal  shape  that  is  drawn  to  a  partial  function  on:  [0,1]  X  [0,1]  -*  [0,1]  X  [0,1]. 

Using  the  inductive  process  to  deHne  the  essence  of  the  fractal 
sequencing,  we  have  the  following  definition: 

(Basis) 

Step  0 

Oo  =  {(0.0, 0.0), (1.0, 0.0)} 

(Induction) 

Step  k 

Label  the  ordered  set  of  2-tuple  points  from  Step  k-1  as: 

n,-.  -  ) 

where 

n  =  4.  1 

Determine  the  new  set  of  ordered  points  for  Step  k  as: 

n,  =  > 

where 

m  4*  +  1 
where 


where 


3,^*4)  =  P‘-', 

(p‘„p*„p‘,)  =  p*-‘,  [jRjJp'-'s 

(P‘iO,P‘iiP‘h)  =  P*-‘3[3R4]p*-‘4 

(P*m-3,P‘m-3.P‘m-l)  =  P*'*.-,  [.-.R.Jp*''. 

and  where  the  relation  [i-iRi]  geometrical  relationship  between  the  two 
end  points  of  the  initiator  line  segment  from  step  h-i  and  the  five  points  of  the 
generator  that  compose  the  ordered  subset  for  step  k  as  above.  For  purposes  of 
brevity,  the  full  functional  definition  for  this  geometrical  relationship  is  explained 
in  detail  in  Chapter  4. 

(Extremal) 

No  point  is  an  element  of  0  unless  it  can  be  shown  to  be  so  in  a  finite 
number  of  applications  of  clauses  1  and  2  above. 

It  is  thus  possible  (but  tedious)  to  rigorously  characterize  the  Koch 
fractal  set  within  the  well  defmed  constructs  of  induction  and  the  functional 
technique. 

6.  Another  Fractal  Set,  Cantor  *b  Dust 

Cantor  developed  a  function’^  which  used  the  same  functional  technique 
as  the  Koch  curve  but  with  a  reverse  twist.  Cantor's  function  takes  an  initiator 
(the  unit  interval  [0,lj)  and  diteolvee  it  into  a  discontinuous  set  which  b  as  rich 
in  points  as  the  interval  (0,lj  but  contains  no  interval  itself  {Ref,  6:pp.  22-23], 
The  points  of  the  set  are  all  dbtinct  but  the  set  has  the  same  cardinality  as  the 
unit  interval.  The  best  description  for  this  set  b  that  it  resembles  a  du$L 

'*  Abo  referred  to  u  Castor's  discoatisuuin  or  Castor'e  triadic  set. 
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Cantor's  Dust  is  difficult  to  demonstrate  on  a  graphics  display  because  it 
quickly  dissolves  below  the  resolution  of  the  display.  Refer  to  Figure  2.2  to 
visualize  the  initiator  and  the  generator.  The  initiator  is  a  line  interval  [Qi^] 


1 


where  0  ^  a  <  6  ^  1  and  the  generator  is  two  intervals  each  —  the  length  of 

3 


the  initiator  such  that  interval  1  is  a ,  o  + 


(6-a)  X  1 


and  interval  2  is 


[b- 


(i-a)x^ 


,6].  After  the  initial  application  of  the  generator  to  the  unit 


1  2 

interval  there  will  be  two  intervals  in  the  current  construction  (0,—)  and  {—,1). 

3  3 


The  second  iteration  of  the  recursive  routine  will  yield  four  intervals 


Every  initiator  that  is  created  by  the  generator  has  an  infinite  sequence 
that  is  begun  at  the  next  application  of  the  generator.  This  causes  a  series  of 
convergent  sequences  to  each  end  point  of  each  initiator.  Thus  each  initiator 


■INITIATOR 


,GgNBfATOE„ 


CANTOR’S  DUST 


Initiating  Structure: 
The  Interval  [0,1] 


Dt  =  0 


HB  =  .6309 


1st  Iteration 
2nd  Iteration 


1 


Initiator  Length  ^ 

1 


.....  Initiator  Length 


0 


3rd  Iteration  ..  .. 


Initiator  Length  2^;^ 


4th  Iteration 


.  Initiator  Length 


Figure  2.2  Cantor’s  Dust. 
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tpawfu  a  convergent  sequence  toward  its  endpoints.  But  for  each  initiator  there 
are  also  two  spawned  sequences  toward  the  center.  It  is  possible  to  imagine  this 
as  an  infmitely  dividing  organism  which  leaves  behind  four  points  (eggs)  each 
time  it  divides  and  then  each  egg  itself  replicates  ‘  ad  infinitum. 

This  functional  method  of  dtssolnn#  a  line  is  very  powerful  as  a  tool  in 
computer  graphics  because  it  can  be  used  to  cause  many  special  effects  from 
ordinary  objects.  Mandelbrot  has  used  variations  of  this  method  to  create  images 
of  star  systems  for  example. 

The  topological  dimension  of  Cantor’s  dust  is  0  and  by  Mandelbrot’s 
dimension  generating  function  we  have  the  fractal  (Hausdori!)  dimension  equal 
to: 

N  /  \D 

g(D)=  eU  =1 

Where  N  =  2 
Where  \ 

Subetituting  and  aoiving : 

\D 

2  X  ~  =  I 

3 

D  =  »  .6309 

log  3 

Cantor’s  dust  was  proven  to  have  a  Hausdorff-Besicovitch  dimension  of 

iSil  .6309  (Ref.  l:pp.  77.78). 
logs 

7.  A  Note  on  the  Concept  of  Similarity  and  Fractals  Sets 

The  use  of  the  mathematical  concept  of  similarity  continually  shows  up 
in  the  investigation  of  fractal  functions.  This  can  be  expected  because  of  the  type 
of  functional  building  tools  that  are  used  for  most  of  these  sets.  The  relationship 
the  generator  to  the  initiator  has  similarity  built  in.  Thus  Mandelbrot  is  able 
to  use  similarity  properties  in  formulating  a  dimension  generating  function  and  in 
making  claims  about  the  set’s  inherent  structure. 
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The  functional  method  invented  by  Koch  and  Cantor  is  but  one  method 
of  determining  t,  fractal  set.  Similarity  fractals  may  eventually  be  grouped  into  a 
class  of  fractals  (important  but  restricted).  This  type  of  functional  method 
provides  a  vehicle  for  the  creation  of  diseipiined  fractal  sets  but  makes  you 
wonder  about  the  rest  of  the  spaet  that  (ills  the  gaps  between  the  standard 
Euclidean  shapes*^,  the  self'similar  fractals  and  the  infinitude  of  R^. 


'*  Much  like  (be  (reiueeiMUiitel  numben  (ir  end  *  for  example  )  in  ftU  ibeir  uncoualftbly  rich 
expuue  fill  the  geiM  between  the  fuailinr  (ud  well  beheved)  seU  in  R. 
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III.  THE  IMPLEMENTATION  OF  FRACTALS  IN  COMPUTER  GRAPHICS 


This  chapter  introduces  the  reader  to  the  practical  aspects  of  implementing 
Fractals  on  a  two  dimensional  computer  graphics  display.  We  are  forced  to 
confront  the  issues  of  economy  of  scale  between  the  inHnite  fractal  function  and 
the  finite  computational  environment  of  the  computer.  Understanding  these 
compromise  is  a  necessary  bridge  to  successful  fractal  programming. 

A.  THE  IMPLEMENTATION  PROBLEM 
1.  Infinite  Recursion,  Stacks  and  Data  sets 

It  is  naive  to  view  the  computer  as  a  truly  infinite  abstract  machine 
which  is  capable  of  may  binary  computation  of  any  length.  Infinity  is  a  concept 
that  when  applied  to  physical  objects  quickly  breaks  down  as  soon  as  that  object 
is  bounded  in  any  way. 

It  is  possible  to  model  infinite  behavior  in  computers  though  mathematics 
(automata  theory)  and  gain  useful  insight  into  possible  capabilities  of  the 
computer  (the  use  of  push  down  automata  in  compilers  for  instance).  But  in 
order  for  automata  theory  to  make  assertions  it  is  dten  necessary  to  make  the 
assumption  that  the  automata  (computer  or  an  abstract  machine)  is  in  fact 
infinite.  When  the  assunq>tion  of  infinity  is  made,  there  are  many  powerful 
mathematical  tools  which  can  be  brought  to  bear  upon  non-intuitive  abstract 
problems  that  would  otherwise  be  functionally*^  intractable  if  the  automata  was 
considered  to  have  an  arbitrarily  large  but  bounded  space.  The  question  arises, 
can  we  consider  the  implementation  of  these  constructs  (insights)  as  valid?  The 
answer  is  yes,  but  only  in  the  context  of  some  bounded  space  (for  instance,  the 
maximum  stack  space  for  a  push  down  compiler). 

The  question  of  how  many  valid  programs  can  be  recognized  by  such  a 
compiler  is  usually  too  difficult  to  determine.  The  compiler's  sdution  to  the 
problem  is  to  use  a  passive  sensor  to  detect  stack  overflow  and  notify  the  user 

Fuaetioiialtr  in  (he  lense  (htl  (he  problem  may  be  decidable  but  the  toluiioo  apace  ia  ao 
Urge  aad  uBdcTuted  that  it  could  not  be  deiermiaed  ia  «  reasonable  amouat  of  time. 


that  his  problem  is  to  large  for  the  current  stack  space.  Although  the  compiler  is 
theoretically  a  recognizer  for  an  infinite  set  of  programs,  the  finiteuess  of  the 
computer  is  the  grounding  factor. 

A  similar  paradigm  exists  for  implementing  mfmite  functions  like  fractals 
in  computers.  Fractal  functions  use  recursion,  randomness,  massive  floating  point 
computations  and  large  amounts  of  primary  memory.  In  order  to  use  the  fractal 
function  productively,  we  must  manage  the  methods  we  use  and  produce  a  finite 
approximation  of  the  fractal  set  we  create.  This  can  be  done  by  a  passive  or 
active  means. 

2.  The  Bounded  Stack  Limits  Recursion 

The  Koch'like  fractal  functional  method  is  by  definition  a  recurrence 
equation.  All  fractal  methods  introduced  in  this  thesis  have  a  similar  recurrence 
definition.  Thus  it  is  impossible  to  avoid  the  use  of  recursion  in  the  production  of 
fractal  graphics^^ 

Recursion  on  most  Von  Neuman  type  computers  is  implemented  on  the 
system  stack  (which  may  be  hardware  (fixed)  or  software  (in  primary  memory 
and  therefore  expandable)).  Such  a  stack  is  always  bounded  by  a  fixed  upper 
limit  of  allocatable  memory  space.  The  formal  definition  of  the  Koch>like  fractal 
method  has  no  recurnve  terminating  event.  In  order  to  use  these  methods,  we 
must  determine  the  precision  that  we  require  and  develop  a  termination  event  to 
signal  the  beginning  of  the  recursive  ascent  when  this  precision  is  met. 

Failure  U)  manage  the  recursive  descent  inevitably  causes  the  exhaustion 
of  the  system  stack  which  stores  the  program  state  on  each  recursive  catl'^  The 
programmer  should  be  careful  to  keep  his  local  variables  at  a  minimum  in  the 
recursive  subroutine.  By  doing  so,  the  total  data  stored  during  each  recursive 
call  is  reduced. 


'*  Why  ihouid  yoo  try  ?  The  r«eursi‘->e  p«rt  of  fuaelioos  U  the  metheoieticeUy  bc«u- 

tiful  upecl  which  makes  them  so  eloquent  and  poweeful. 

This  could  be  the  recursive  tertniaalioii  eveat  if  you  sie  not  csjreful.  It  might  be  a  good 
idea  to  try  this  experiaieBtaUy  oa  your  computer  to  determiac  this  maximum  recursive  doKeat 
disuace. 


3.  The  Computer  Graphics  Set  Paradigm 

The  paradigm  that  is  used  for  displaying  objects  on  the  raster  graphics 
screen  is  inherently  discrete  and  two-dimensional.  A  typical  system  consists  of 
prhmtive§  that  allow  the  user  a  view  of  his  modeling  world  as  a  largely 
unrestricted  three  dimensional  space.  These  primitives  are  limited  in  their  power 
and  are  usually  based  on  the  Euclidean  geometry  (lines,  points,  polygons  etc.). 
The  user  is  required  to  supply  oieiwfag  fcrameUrt  that  deflne  a  limited  three 
dimensional  viewing  space.  These  commands  are  processed  by  the  graphics 
system  which  projt  cU-  the  objects  contained  in  the  viewing  space  onto  a  two- 
dimensional  space  thut  is  the  display  screen.  The  display  screen  can  be  divided 
into  discrete  entities  called  pixtUt  each  of  which  represents  one  point  on  the 
screen.  The  number  of  pixels  defines  the  resolution  of  the  screen. 

Given  the  graphics  paradigm,  it  is  natural  to  model  objects  that  are 
defmed  on  the  display  as  a  collection  of  discrete  points.  When  viewed  on  the 
screen,  these  pictures  can  accurately  model  objects  of  considerable  complexity 
[Ref.  ij.  This  view  is  a  departure  from  the  normal  view  of  the  graphics 
application  programmer.  The  application  programmer  views  the  objects  he 
creates  through  a  collection  of  Euclidean  geometric  primitives  (lines  polygons 
etc.)  that  are  abstracted  from  the  actual  display**.  Both  views  are  useful  models 
in  describing  objects  but  for  the  fractal  graphics  pr<%rammer  the  former  b  the 
more  powerful  because  it  neatly  nuqjs  to  the  fractal  method  of  object 
construction. 

For  thb  study,  we  view  the  world  and  the  objects  it  contains  as 
collections  (sets)  of  discrete  three  dimcnsbnal  pmnts.  When  an  object  b  built 
through  fractal  techniques,  the  length  of  the  pixel  b  a  natural  termination  point 
for  the  fractal  recursive  process.  Thb  yields  an  attractive  byective  mapping  from 
our  object  to  the  dbplay  screen. 


Ail  computer  iraphic*  u  ea  abtlractioa;  reaiUm  i*  achieved  throttfh  dKeplton  dt  the  eye 
by  a  very  »mail  coltccltoa  of  cobred  powu.  Our  goal  Utco,  u  to  ftad  the  moat  cQkieat  method  of 
defiaiag  thoae  poiota. 
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4.  Summary 

From  our  previous  discussion,  it  should  be  obvious  that  fractal  sets  do 
not  exist  in  our  computer.  The  sets  that  we  create  are  finite  approximations  of 
the  actual  set.  The  same  is  true  of  the  fractal  pictures  that  we  create  and  to 
which  we  ascribe  fractal  dimensioru.  It  is  an  illusion  of  the  eye  that  we  create  by 
taking  the  fractal  computations  below  the  resolution  of  the  screen. 

B.  MAPPING  FRACTALS  TO  A  BOUNDED  SPACE 

In  view  of  the  graphics  set  paradigm,  it  behooves  us  to  develop  a  completely 
bounded  set  space  which  has  a  one  to  one  correspondence  between  the  points 
computed  by  the  fractal  equations  and  the  entities  of  the  screen.  This  abstraction 
is  appropriate  because  these  entities  called  pixeb  (which  are  nothing  more  than 
colored  points)  are  the  only  components  of  the  picture*^  The  only  reason  for  not 
using,  thb  methodology  before  is  that  it  was  functionally  difficult  to  compute 
(without  some  simplifying  abstractions)  the  large  number  of  pixels  in  the  display 
set 

Fractal  Recursive  Termination  Event 

Most  graphics  software  packages  provide  a  means  by  which  the  user  can 
define  the  space  on  the  display  screen  onto  which  he  wishes  to  map.  This  is  done 
by  providing  to  the  graphics  system,  viewing  parameters,  which  define  a  bounded 
three  dimensional  space  which  b  then  projected  onto  the  two  dimensional  screen. 

Many  mappings  are  possible  within  the  above  constructs  but  for  i:he 
fractal  pri^rammer  U  b  convenient  to  establbh  a  mapping  that  provides  for  a 
one  to  one  (or  byective)  relationship  between  the  3^  pixel  set  and  the  real>valued 
coordinate  space  in  which  we  compute  a  fractal  object.  We  develop  thb  byective 
relaltonship  by  defining  a  fractal  part  to  be  determined  when  the  dbunce 
between  generating  points  (in  real  coordiuotee)  b  less  than  the  dbtance  of  one 
pixel  in  the  mapping  of  real  to  screen  coordinates  (SCR),  Figure  3.1.  Thb 
mapping  b  explicit  when  the  windov,  vieaipori  and  s^dipping  are  defined.  Thb 
mapping  limits  the  number  of  pixeb  that  can  be  associated  with  a  given  fractal 
initiolor  to  a  one  to  one  mapping  of  fractal  parts  to  the  pixel  set.  The  sise  of  the 

'*  la  fact,  tbu  is  Uie  way  Iks'  fAiciutesa  fraf^ks  absintels  such  u  a  tuw  an  napped  to  tge 
sema.  At  soanc  poki  ia  tit*  ditpUy  pcoecaa  bis  must  talu>  pUc«. 
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pixel  is  a  cubic  space  which  clumps  together  the  inPinite  fractal  space  into  a 
discrete  number  of  cubic  spaces  equal  to  the  number  of  pixels  in  the  pixel  set. 

The  fractal  programmer  should  choose  his  viewport  and  window  careftilly 
so  as  to  elici*  the  most  attractive  mapping  possible  between  his  object  and  the 
screen.  The  viewport  (SCR)  rectangle  should  be  defined  such  that  the  ratio  of 
the  X  distance  (horizontal)  to  the  Y  distance  (vertical)  is  equal  to  the  ratio  of  the 
X-Y  distances  of  his  (real-valued)  window.  The  Z-distance  (front  and  back 
clipping  plane  distance)  should  be  equal  to  the  distance  in  the  window 
coordinates  that  defines  the  length  of  the  pixel  multiplied  by  the  desired  Z-depth, 
Figure  3.2.  Formally: 

_  Xn’to 

Ysce  ^  wiD 


Recursive  Terepination  Event 

Recursive  Xteretion  i 


Thm  LmagtSi  Lmmw 

Titma  rhm  FSXSL_SXSS,  So  Mmp 
T&*  Poiat-m  To  Ttsa  FixmJ  Smt. 


Recursive  Iteration  i  At 
The  POXliT  PLOT  Routine 


TS*  FiattJm  lf#r#  Umppxei  To 
Thtir  Corratfpottiiiog  Foxicioa 
40  T&*  FSaimi  Sat. 


Figure  3.2.  The  Fractal  Recursive  Termination  Event. 


^  ■*•■»*<  ».  «  »  *■  y  •  ►  *  7.  • 
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If  the  above  relationship  holds  then  the  size  of  the  pixel  in  world  coordinates  is 
equal  to  the  X  (or  Y)  world  distance  divided  by  the  number  of  pixels,  Figure  3.2. 


PIX_SIZE  = 

X5CJ1 

The  front  and  back  Z-clipping  planes  should  be  established  such  that  the  distance 
from  front  to  back  is  equal  to  the  desired  Z  depth  (expressed  in  pixeb)  multiplied 
by  the  pixel  size  (in  world  coordinates),  Figure  3.2. 

If  the  above  viewing  relationships  hold,  then  the  inclusion  of  a  simple 

check: 


IF  (Initiator  J)lstance  <  PK  JIZE)  THEN 
return; 

ELSE 

continue; 


terminates  the  current  recursive  descent  and  begins  backtracking. 

2.  Memory  Requirements 

The  amount  of  main  memory  required  by  the  fractal  programmer  is 
directly  proportional  to  the  size  of  the  pixel  space  that  is  being  used  to  display  a 
fractal  figure  and  the  sophistication  of  the  desired  display, 
a.  Data  Locality  During  Recursive  Descent 

If  the  fractal  in  question  is  a  1-2  dimensional  display  that  is 
displayed  as  a  collection  of  points  on  a  two  dimensional  plane  then  the 
requirement  to  store  the  3-tuple  points  is  eliminated  {(z  ,v  ,2 )  coordinates  which 
represent  a  point  in  R^).  The  points  can  be  computed  by  the  function  and 
displayed  on  the  screen  in  the  immediate  mode  and  then  destroyed.  The  only 
requirement  for  memory  is  the  data  that  is  germane  to  the  program  (“globals** 
and  thus  fixed  at  run  time)  and  the  amount  of  data  pushed  onto  the  run-time 
stack  for  the  recursive  calls  up  to  the  point  of  the  deepest  recursive  level.  This 
requirement  b  minimal  and  does  not  pr^nt  any  significant  limitations. 

An  example  cf  such  a  fractal  b  the  Koch  curve  The  algorithm 
presented  in  chapter  4  uses  the  data  that  b  local  to  the  subroutine  that  b 
computing  the  current  generator  with  the  input  coordinates  of  the  initiator.  The 
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inductive  nature  of  the  fractal  method  provides  data  independence  so  that  the 
programmer  does  not  have  to  have  available  to  his  subroutine  all  points 
computed  thus  far. 

b.  Memory  Requirements  for  the  Fractal  Set  Paradigm 

If  the  programmer  requires  a  more  sophisticated  display  that  includes 
a  requirement  for  hidden  surface  removal  and  (or)  lighting  enhancement  then  the 
fractal  method  becomes  put  a  part  of  the  overall  display  process.  The 
programmer  has  the  fractal  process  compute  the  points  of  the  picture  and  these 
points  are  stored  in  a  memory  structure.  This  structure  is  processed  by  a  hidden 
surface  algorithm  and  a  lighting  algorithm  and  is  then  projected  onto  the  screen. 

The  ideal  Fractal  computer  would  have  an  exhaustive  memory  and 
unlimited  <*omputing  power  to  be  able  to  allow  us  to  store  the  entire  pixel  set  in 
memory.  This  however,  can  stretch  the  capabilities  of  most  present  day 
computers.  For  example;  if  the  pixel  set  is  defined  such  that  its  dimensions  are 
1000x1000x1000  (which  is  certainly  reasonable),  the  amount  of  storage  required 
would  be  an  array  of  10*  x  24  bits*®  or  24  billion  bits. 

Current  techniques  exist  where  this  storage  can  be  minimized.  A  Z- 
buffer  array  can  be  used  to  store  the  current  Z  coordinate  of  the  forward  most 
displayed  pbcel.  By  using  this  method  the  fractal  computation  can  be  made  and 
then  a  hidden  surface  computation  can  be  immediately  invoked  to  check  the 
viaibility  of  the  point  by  checking  it  against  any  other  point  in  the  same  position 
on  the  X'Y  plane.  This  technique  reduces  the  above  space  requirements  to 
10®  X  24  bits  +  (Z- buffer  =  10®  x  32  bits)**  or  56  million  bits. 

A  fractal  programmer  must  manage  his  memory  resources  carefully 
in  order  to  maximize  the  computational  resources  available  to  him.  This  usually 
requires  a  tradeoff  between  efficiency,  realism  and  the  memory  space  utilized. 
The  algorithms  for  hidden  surface  and  lighting  are  well  covered  :n  the  literature 
and  although  they  are  integral  aspects  of  the  fractal  realism  issue  they  arc  not 
germane  to  this  study.  Throughout  the  remainder  of  this  study,  wc  are  not 

”  34  biu  for  «  machine  uMuminf  the  RGB  color  tyslem  with  34  bit  plages.  Each  color;  red. 
|r«eg  and  blue  would  then  have  6  bits  of  precision. 

Assuming  the  Z  coordinate  u  stored  as  a  33  bit  floalins  point  number;  tbit  number  can  be 
further  reduced  if  we  restrict  the  Z  precision. 


concerned  with  the  methods  of  minimizing  the  the  pixel  set.  We  assume  that  we 
have  the  full  three  dimensional  set. 

3.  Concurrency 

A  question  arises,  is  it  p<»sible  to  generate  fractals  in  Rtai  Time  We 
have  enough  knowledge  at  this  juncture  to  discuss  the  possibility  for  concurrent 
operations  during  the  fractal  recursive  descent.  There  are  two  basic  strategies 
that  can  used  if  concurrent  operations  are  available  on  the  computer.  We  could 
assign  a  process  to  each  of  the  initial  Initiators  and  process  each  to  its  desired 
precision.  We  could  also  allow  the  imagined  processor  to  have  a  means  by  which 
at  any  level  of  the  recursive  descent,  we  could  spawn  a  process  to  apply  the 
generator  to  the  local  initiator^'^. 

The  local  data  independence  of  the  fractal  recursive  function  allows  such 
a  spawning  because  the  application  cf  the  generator  to  an  initiator  is  local  to  the 
initiator.  This  application  is  a  result  of  the  inductive  process  of  the  K-1  steps 
that  led  up  to  the  step.  The  mathematical  process  of  applying  the  generator 
at  the  level  needs  only  the  initiator  data  from  the  K-1**  step  and  requires  no 
knowledge  of  the  entire  fractal  set.  The  power  of  induction  is  then 
computationally  realized.  The  number  of  calculations  before  the  entire  pixel  set  is 
determined  is  not  greater  than  the  number  of  computations  required  to  apply  a 
generator  time.s  the  maximum  recursive  descent  distance  from  any  initiator  to  the 
recursive  termination  event.  A  short  algorithm  describing  the  processing  aspects 
of  the  second  method  is  shown  in  Figure  3.3. 


**  To  b««t  the  bumu  eye  and  achieve  the  ultimate  graphics  illusion,  typically  by  completing 
all  computattans  and  screen  I/O  prior  to  the  refresh  rate  (for  a  SOHZ  display)  of  most 

raster  graphics  displays. 

**  This  of  course  would  require  an  enormous  computational  power  that  does  not  currently  ex¬ 
ist,  so  in  reality  the  spawning  would  have  to  be  bounded  in  tome  way. 


mainO 

LOAD  INITIATOR  AND  GENERATOR  DATA 

concurrent; 

generate(INITIATOR  1); 
generatejlNITlATOR  2); 
generatejiNITUTOR  3); 


generate(INITIATOR  n); 
end  concurrent; 
end  main; 

generate(INITIATOR  I) 

BUILD  GENERATOR  FROM  INITIATOR 

if  (Distance  <  PIX_8ize)  then 
PLOT  POINT 
return; 
endif 

concurrent; 

generate(INITIATOR  1); 
generate(INITIATOR  2); 
generatejlNITIATOR  3); 


generate(INlTIATOR  m); 
end  concurrent; 

end  generate; 


Figure  3.3.  An  Algorithm  Which  Utilizes  Concurrent  Operations. 


IV.  FRACTAL  COMPUTATION  IN  R* 


The  algorithm  introduced  in  the  following  section  is  capable  of  computing  a 
very  broad  class  of  Kochrlikt  fractal  curves  within  R^.  It  provides  the  fractal 
graphics  programmer  with  the  basic  tools  for  fractal  computation  and  an 
algorithmic  template  that  can  be  used  for  many  applications.  The  graphics 
programmer  needs  to  fully  understand  fractal  programming  within  R^  before  he 
attempts  the  more  complicated  concepts  of  fractal  terrain  modeling  in  R^. 


A.  THE  GEOMETRY  OF  INITIATOR  GENERATOR 

The  geometric  relationship  between  the  Koch  curve  initiator  and  generator 
can  be  described  through  a  very  simple  set  of  data  that  captures  the  essence  of 
the  Koch  curve.  The  method  introduced  also  allows  us  to  vary  the  data  which 
defmes  the  generator  and  compute  many  different  fractal  shapes^^. 

The  general  strategy  for  computing  generator  points  from  a  set  of  initiator 
points  is  to  determine  two  lines  that  intersect  at  an  unknown  generator  point, 
Figure  4.1.a.  The  unknown  generator  point  is  defined  by  constant  relationships 
between  the  initiator  and  generator.  The  first  line  is  the  perpendicular  from  the 
generator  point  to  the  initiator  line.  The  second  line  is  formed  between  the  first 
point  of  the  initiator  and  the  unknown  generator  point.  All  data  to  compute  the 
line  equations  are  derivable  from  the  two  endpoints  of  the  initiator  or  from 
constant  initiator/generator  ratios.  For  simplicity's  sake,  we  ignore  the  divide-by- 
zero  problems  that  are  encountered  when  any  of  the  lines  are  parallel  to  the  X  or 
Y  axis.  These  situations  only  simplify  the  computations  and  their  solution  is 
demonstrated  in  Appendix  A. 

If  we  label  the  endpoints  of  the  initiator  as  P|»(X),Y|}  and  P2=(X2,Y3) 
(Figure  4.1. a)  then  the  slope  of  the  initiator  line  is: 


Slope.iiiit 


Xj- A, 


We  must  be  cvefut  in  our  lettninolosy  because  the  relationship  described  can  draw  shapes 
that  do  not  avoid  self  iaierseclion  and  thus  roust  be  coaeidered  quasi*fractal. 
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The  slope  of  the  perpendicular  intercept  line  is  the  negative  inverse  of  the 
slope  of  the  initiator,  hence  the  slope  of  this  lins;  is: 


Slope.perp  » 


-1 

F’.ope.iait 


The  intercept  point  on  the  initiator  can  be  determined  by  using  a  constant 
dbtance  ratio  as  shown  by  the  following  equations: 


X. perp  = 

Y. perp  = 


X]  +  {Generator. ratio.eonstant  x  X2) 
1  +  Gef%erator.ratio.con$tant 

Yi  +  {Generator. ratio.eonetant  x  Tj) 
1  +  Generator.ratio.eonatant 


The  value  of  ''Generator,ratio.eoneiwti"  is  a  fixed  constant  (although  it  might  be 
interesting  to  randomize  it)  where  you  determine  at  what  point  the  generator 
intercept  point  intersects  the  initiator  and  express  it  as  demonstrated  in  Figure 
4.1.b.  This  ratio  can  be  determined  graphically,  through  hand  calculation  or  via 
an  interactive  automated  means. 

With  the  slope  and  a  point  of  the  line  (X.perp,Y.perp)  we  can  determine  the 
line  equation  for  the  perpendicular  line  by  determining  the  Y  intercept: 

Y.interct^t.peip  Y.perp  -  (Slope.perp  x  X.perp) 

Thb  yields  the  line  equation  for  the  perpendicular  line: 

Y  =  (Slope.perp  x  X)  -f  Y.intercept.perp 

To  deteii  dne  the  generator  line  equation  for  the  line  segment  which  connects 
the  first  point  of  the  initiator  and  the  unknown  generator  point,  we  need 
constant  information  about  the  angle  between  the  initiator  and  this  line.  This 
angle  is  always  constant  with  respect  to  the  initiator  and  like  the  ratio 
information  above,  it  is  recorded  as  a  constant  at  run  time  (see  Figure  4.1.c), 
note  that  the  angle  may  be  positive  or  negative. 


I 

Generator  Line  •'  j 


Unknown 

» 

Generator  Point 


^(x^,ya) 

—Initiator 


^■W  WL 


'Perpendicular  Intercept 


Fig  4.1. a  Intersecting  Lines  Determine  a  Generator  Point. 
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Fig  4.1.b.  The  Generator 

Ratio  Constant 

Point  with  positive  angle  Point  with  negative  angle 


B 


Fig  4.I.C.  The  Angled  '''qP“ 


Figure  4.1  Building  the  Generator  with  Xnterse\:ting  Lines. 


W  *n. V  i,r 


We  record  the  data  about  the  angle  $  as  the  tangent  of  the  angle.  With  this 
information,  the  slope  of  the  generator  line  can  be  determined  with  the  following 
equation: 


Slope.gen  = 


Tan  6  -f-  Slope.init 
(1  -  Tan  )  X  Slope.init 


The  Y  intercept  for  the  generator  line  can  now  be  determined: 


Y.intercept.gen  =  Yj  -  (Slope.gen  x  Xj) 


This  yields  the  line  equation  for  the  generator  line: 

Y  =  (Slope.gen  x  X)  +  Y.intercept.gen 


The  Cartesian  points  of  the  unknown  generator  point  can  now  be  determined 
by  intersecting  the  two  line  equations  and  solving  for  (Figure  4.1.a)  then 
substituting  into  one  of  the  line  equations  and  solving  for  Yg^g.  The 

equations  follow: 


Y.intercept.perp  -• 


Y.intercept.gen 
Slope.gen  -  Slope.perp 


=  (Slope.gen  x  Xg,g)  +  Y.intercept.gen 


gen/ 


The  constant  data  for  the  Koch  curve  that  corresponds  to  this  geometric 
method  is  illustrated  in  Figure  4.2. 

There  are  many  different  ways  to  build  a  generator  given  an  initiator  using 
standard  geometric  constructs  but  the  method  introduced  allows  experimentation 
with  the  Koch  function  to  discover  new  shapes.  By  varying  the  tangent  of  6  and 
the  fixed  ratio  of  Figure  4.1  we  can  describe  new  generator  constructions.  These 
new  constructions  can  be  used  in  the  same  algorithms  that  compute  the  Koch 
curve.  We  initiate  the  recursion  on  a  generating  structure  built  of  line  segment 
initiators  and  allow  the  recursion  to  progress  until  a  terminating  event  creating 
many  diverse  shapes.  Figures  4.6  through  4.0  demonstrate  some  of  these  shapes. 


**  Tlie  leaglh  of  •  pixel  or  u  vbilrary  line  length. 
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Figure  4.2  Ratio  Constants  for  the  Koch  Curve  Generator. 

B.  THE  MID-POINT  DISPLACEMENT  TECHNIQUE 

Mandelbrot  [Ref.  l;pp.  43  and  pp.  233-234]  uses  an  alternate  technique  to 
draw  the  Koch-like  curves  that  is  equally  valid.  He  calls  his  technique  mid-point 
displacement  because  he  determines  the  generator  via  fixed  relationships  that 
displace  a  point  from  the  mid-point  of  the  initiator.  This  method  uses  many  of 
the  same  geometric  relationships  that  are  used  above  but  provides  a  different 
progression  to  the  method  of  building  the  generator.  This  new  view  allows  us  to 
look  at  the  relationship  between  the  initiator  and  generator  in  a  slightly  different 
light.  By  so  doing  we  are  provided  new  insight  as  to  how  we  might  alter  the 
relationship  to  create  new  images  and  sets. 

The  mid-point  displacement  technique  can  be  useful  for  two  other  reasons.  It 
is  the  best  known  method  for  fractal  set  building  and  because  of  this,  it  facilitates 
communication  between  fractal  programmers.  Most  the  terrain  models  that 
have  been  developed  use  the  mid-point  technique.  It  also  provides  an  easier 
method  to  avoid  line  intersection. 
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The  Midpoint  Displacement 

Technique 


Figure  4. 3. a  First  Midpoint  Displacement. 


Figure  4.3.b  Second  Midpoint  Displacement. 


Figure  4.3.  The  Koch  Generator  Using  Midpoint  Displacement. 


The  mid'poiiU  dtspSacement  method  is  demonstrated  in  Figure  4.3.  The 
method  progresses  by  taking  the  initial  initiator  and  applying  the  first  midpoint 
displacement.  This  yields  the  figure  demonstrated  in  Figure  4.3.a.  The  next  step 
performs  a  mid-point  displacement  on  the  left  initiator  created  by  step  1.  This 
yields  the  figure  demonstrated  in  Figure  4.3.b.  The  third  and  final  step  is  to 
replace  the  right  initiator  created  by  step  1  as  demonstrated  in  Figure  4.3.c. 


The  inversion  of  the  direction  of  the  mid-point  application  (in  Figure  4.3.a 
the  displacement  is  above  the  initial  initiator  where  in  Figure  4.3.b  it  is  below) 
can  be  accomplished  with  a  single  computational  procedure.  We  only  need  to 
invert  the  position  of  point  1  and  point  2  in  the  parameters  of  that  procedure. 
The  parameter  mversion  changes  the  orientation  of  the  initiator  in  space  with 
respect  to  the  computation  cf  the  midpoint  displacement.  The  procedure  blindly 
computes  the  mid-point  displacement  relative  to  a  fixed  relationship  to  the 
initiator  input  points.  This  operation  implicitly  defmes  the  orientation  of  the 
generator  in  space. 

The  geometry  for  computing  the  midpoint  given  two  initiator  points  can  be 
computed  in  precisely  the  same  manner  as  the  intersecting  line  algorithm  (The 
Koch  midpoint  ratios  are  1.0  and  the  angle  $  is  .437  radians).  An  alternative 
method  is  to  use  the  equations  of  a  line  (or  a  plane)  normal.  The  second  method 
(utilizing  the  normal)  provides  a  geometric  relationship  which  is  intuitively 
appealing.  Its  appeal  comes  from  the  desire  to  modify  the  length  of  the 
displacement  relative  to  the  initiator  with  a  random  scaling  factor. 

The  mid-point  displacement  technique  has  some  advantages  over ‘the  line 
internetion  algorithm.  Random  modification  of  the  length  of  the  displacement 
along  the  normal  (from  the  computed  generator  point  to  the  initiator  using  the 
line  intersection  algorithm)  is  not  intuitively  appealing.  It  requires  a  translation 
of  the  desired  displacement  into  an  angle  (the  angle  d  (Figure  4.1.c)  between  the 
initiator  line  and  the  unknown  generator  intercept  line).  The  control  of  that 
angle  is  less  intuitive  than  the  control  of  a  displacement  length.  The  geometry 
for  mid-point  displacement  using  the  initiaUvr  normal  is  introduced  in  Chapter  5. 

C.  A  KOCH-UKE  FRACTAL  ALGORITHM 

Implementing  the  above  function  U  a  relatively  easy  process  that  k 
demonstrated  in  thb  section  and  Appendix  A  via  a  gradual  unwroppinp  of  the 
layers  df  complexity  that  are  required  to  successfully  implement  the  algorithm. 
The  algorithm  roughly  follows  the  template  used  in  Chapter  3  to  demonstrate 
concurrent  processing.  A  C-like  language  is  used  for  the  algorithms. 

The  first  algorithm  (Figure  4.4)  is  a  template  that  delineates  the  basic 
processing  steps.  This  recursive  process  b  typical  of  fractal  functions  and  can  be 


used  as  a  template  for  many  fractal  programs.  The  second  algorithm  (Figure  4.5) 
is  an  expansion  of  the  first  and  demonstrates  the  replacement  of  a  given  initiator 
using  the  line  intersection  method. 

Appendix  A  is  a  complete  Fractal  program.  Thu  program  was  used  to 
produce  the  data  for  Figure  2.1  and  Figures  4.6  through  4.0.  This  program 
demonstrates  the  precautions  that  must  be  taken  to  avoid  dividt~b}f‘Zero  when 
lines  are  parallel  to  the  X  or  Y  axis. 


mainO 

i 


Load  Initiator  Coordinates; 

Load  Generator  Relationship  Values; 

for  1=1;  1<=  Number  of  Initiators;  1=1+1; 

{ 

generate{X(I)„Y(I),JC(I),.Y{I)j) 

} 


generate(Xi,Y  iXjiYj) 

{ 


Determine  Distance  Between  the  Endpoints  of  the  Initiator 

if  (DIST  <  Pixel.length)  return; 

Replace  Initiator  with  the  Computed  Generator; 

Load  the  New  Initiator  Data  into  Local  Generator  Array 


for  J=l;  J<=  Number  of  Generator  Segments;  J=J+r, 

{ 

gcnerate(Xgen(J),,  Ygen(J)|,  Xgen(J)j,  Yg€a(J)2); 

> 


Figure  4.4.  High-level  View  of  the  Koch  Algorithm. 
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mainO 

{ 


Load  Initiator  Coordinates; 

Load  Generator  Relationship  Values; 

for  1=1;  I<=  Number  of  Initiators;  I=l-f  1; 

{ 

generate(X(I)„Y(l),^(I)j.y(I)j) 

) 


generale(X|  ;Y  |  .Xj.Y  2) 


j*  Dettrmine  Diatonee  Between  the  Endpcini*  oj  the  /nitiaior  */ 
DIST  =  s<pt{(X2-X,)*^2  4-  (Y2-Yi)**2); 

/•  //  Diatanee  i$  teas  than.  Piiei  lenfth ;  Piot  and  Return  */ 
if  (DIST  <  PixelJength) 

{ 

plot.pointO:  /*  Y<Htr  Grophie#  Point  Plotting  Routine  */ 

return;  /*  Point  I  and  2  plot  the  same  pixel  */ 

) 

/*  Load  The  Endpeinta  of  the  imtiotor  into  the  Generator  Array  */ 
Generator.XfOj  =  Xj*. 

Generator. vjoj  ^  Yf  i 

Generator.XjNumbcr.of.generator.pojnts  +  1|  =  X-; 
Getterator,YjNumber.of4eacraior.poiDts  +  ij  «  Y3; 

/*  Determine  Slope  of  ike  Initiator  *{ 

SIopc,wit  =  (Y3- Y|)  /  (Xa-X,); 


Figure  4.5.  Detailed  View  qH  the  Koch  Algorithm. 


j*  Calculate  the  Unknown  Generator  Point  via  Intereeeting  Linee  */ 
for  J=l;  J<=  Number.of.gC3xerator. points;  J=J+1; 

{ 

/*  Determine  the  Generator  Point  intercept  on  the  Initiator  */ 

X. perp  =  {Xii- (Generator. ratio.eonstant  [3]  *  X2))  / 

(1-^Generator.ratio.eonstant  [J]); 

Y. perp  =  (Y 1+ (Generator. ratio.eoiutant  [3]  *  Yj))  / 

(l+Gerurator.ratio.eonstant  [J]); 

/*  Determine  the  Slope  of  the  Perpendicular  Line  */ 

Slope.perp  =  (-1  /  Slope.init); 

/*  Determine  the  Y -intercept  of  the  Perpendicular  Line  */ 
Y.intercept.perp  =  Y.perp  -  (Slope.perp  *  X.perp); 

/*  Determine  the  Slope  of  the  Generator  Line  */ 

Slope.gen  =  Generator.tan.theda  [J]  +  Slope.init)  / 
(l-Gcnerator.tan.theda  [J]  *  Slope.init); 

!*  Determine  the  Y -intercept  of  the  Generator  Line  */ 
Y.intercept.gen  =  Yj  -  (Slope.gen  *  Xj); 

/*  Determine  the  Unknown  Generator  Point  */ 

Generator.X(J]  =  (Y.intercept.perp  -  Y.intercept.gen)  / 

(Slope.gen  -  Slope.perp); 

Generator. Y(Jj  ==  Slope.gen  *  Generator.X[J]  +  Y.intercept.gen: 


} 

for  K=0;  K<=  Number  .of .generator, points  +  1;  K=K+1; 

{ 

generate(Generator.X(K],Generator.Y[K), 

Generator.X(K+ 1], Generator.  Y(K+ 1)) ; 

} 

}  /*  End  Generate  */ 


Figure  4.5.  Detailed  View  of  the  Koch  Algorithm  (continued). 


D.  IMPLEMENTATION  STRATEGIES 

There  are  numerous  ways  to  display  the  fractal  shapes  that  the  above 
algorithm  is  capable  of  computing.  The  graphics  primitives  required  are  limited 
to  the  standard  initiation  and  termination  conunands  coupled  with  the  ability  to 
plot  a  point  (or  alternatively  a  line).  Any  raster  graphics  system,  plotter  or 
similar  technology  sufHces. 

The  algorithm  can  be  extended  to  include: 

-  Online  generator  draw'ing  to  compose  a  generator  relationship  visually. 

-  Rotation  in  3  dimensions  (if  your  system  has  this  capability). 

-  Variation  of  the  inductive  application  of  the  generator  by  the  inclusion  of 
randomness  with  respect  to  the  generator  constants. 

Figures  4.6  through  4.9  represent  a  few  of  the  shapes  that  this  algorithm  can 
compute.  Each  figure  has  the  generator  data  used  to  compute  the  shape  and  a 
progression  that  shows  the  first  two  recursive  iterations. 

E.  SUMMARY 

The  line  intersection  algorithm  as  it  stands  is  not  very  useful  for  the 
production  of  graphics  images  of  realistically  textured  terrain.  Its  importance 
results  from  its  encapsulation  of  the  essence  of  the  non-random  inductive  fractal 
method.  This  algorithm  demonstrates  the  idea  and  intent  of  fractal  functions  and 
their  implementation  within  computer  graphics.  The  potential  fractal 
programmer  must  throughly  understand  the  salient  parts  of  this  chapter  before 
successfully  attempting  fractal  images  in  three  dimensions  (i.e.  before  climbing  the 
mountain  Chapter  5). 


Boxes  Ad  Infinitum 


f  ¥ 


i 


An  Exaggerated  Koch  Curve 


GENERATOR  DATA 

Number  of  poiatm  =  3 
Peiatl t 

Aaglo  a  0.0000  rmdm 
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Figure  4.8.  An  Exaggerated  Koch  Curve. 
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Figure  43.  A  Plane  Filling  C«rve. 


V.  FRACTAL  GEOMETRY  FOR  GRAPHICS  TERRAIN 


One  of  the  most  widely  recognized  fractal  images  found  in  the  literature  is  of 
the  mountain  scene.  This  type  of  terrain  modeling  is  perfectly  attuned  to  the 
fractal  technique.  The  reason  for  this  is  that  mountains  are  highly  irregular 
shapes,  with  a  rough  but  consistent  texture  when  viewed  from  a  distant  vantage 
point.  It  is  appropriate  then,  to  introduce  graphics  terrain  simulation  techniques 
through  this  model. 

This  chapter  describes  the  theory  and  techniques  of  simulating  mountainous 
terrain  with  computer  graphics.  It  provides  the  blueprint  for  fractal  graphics 
programming  within  by  providing  general  tools  and  a  methodology  that  is 
easily  adapted  to  many  other  modeling  needs. 

A.  MODELING  MOUNTAINOUS  TERRAIN 

For  the  programmer  who  fully  understands  the  essence  of  the  method  of 
fractal  programming  introduced  in  Chapter  4,  the  movement  into  programming 
in  is  not  difficult.  The  primary  differences  lie  in  the  quantum  jump  in 
computing  resources  that  are  required  and  the  requirement  to  perform  the 
generator  geometry  in  R*  versus  R*.  The  theory  and  technique  of  fractals  does 
not  change  substantially. 

Chapter  3  provided  a  rough  framework  to  begin  the  coalescence  of  fractal 
progranuning  into  a  workable  technique.  We  need  to  develop  a  number  of  toots 
from  that  chapter  and  use  standard  computer  graphics  techniques  to  manage 
those  tools.  To  this  framework,  we  add  new  fractal  functions  which  provide  the 
texture  of  realism  for  our  simulated  mountain. 

L  The  Artist's  Model 

One  way  for  an  artist  to  build  a  physical  relief  model  of  a  mountain  is  to 
use  a  framework  to  provide  structure  to  the  model  and  a  texturizing  elay  to 
provide  realism.  The  artist  might  use  chicken  wire  on  top  of  small  boxes  as  the 
frame  with  modeling  clay  as  the  texturizing  dement.  His  choice  of  clay  is 
predicated  by  the  type  of  look  that  he  wants  to  achieve.  The  chicken  wire 


provides  an  inexpensive  and  disguised  method  to  quickly  build  the  mountainous 
shape  and  structure.  This  method  minimizes  the  cost  and  time  to  build  up  the 
clay. 

The  artist  continues  the  modeling  process  after  the  development  of  the 
basic  mountain  shape  to  achieve  hues  and  contrast  in  the  coloration.  He  might 
achieve  this  by  the  use  of  natural  lighting  to  cast  shadows  or  by  a  careful 
painting  of  prominent  features. 

2.  The  Fractal  Programmer’s  Model 

There  is  very  little  in  science  that  is  truly  new  or  innovative.  We  borrow 
the  essence  of  the  above  idea  to  guide  us  in  developing  a  model  for  the  discrete 
computation  of  our  two  dimensional  picture  of  the  mountain.  This  section 
describes  the  process  intuitively  and  leaves  the  implementation  details  to  later 
sections. 


a.  The  Lattice  Control  Structure 

The  pixel  space  that  we  developed  in  previous  chapters  can  be 
divided  into  discrete  cubic  units  by  use  of  a  concept  from  mathematics  called  a 
lattiee  (in  our  case  we  can  view  it  as  three  dimensional  graph  paper).  This  lattice 
serves  as  our  controlling  structure,  the  equivalent  of  the  chicken  wire  structure 
above.  It  is  beneficial  to  build  the  lattice  as  a  structure  with  well-formed 
relationships,  where  the  number  of  lines  evenly  divides  the  boundaries  of  the  pixel 
space  and  each  line  is  a  constant  distance  from  its  neighbors.  By  this  method,  we 
do  not  have  to  store  the  lattice  but  can  express  it  as  a  mathematical  function  of 
the  pixel  space. 

The  lattice  can  be  very  useful  in  developing  a  rough  approximation 
of  the  mountain  that  we  wish  to  model.  This  can  be  done  in  many  diHerent  ways 
but  should  result  in  a  itiek  frame  model  of  the  mountain  (a  connected  polygon 
mesh  like  that  of  Figure  5.1). 

The  frame  can  be  deveio])ed  through  an  online  graphics  interface 
that  allows  the  progranuner  to  select  a  gtound  letfd  plane  of  the  lattice  and 
provide  a  means  to  visually  seU*ct  |K>ints  for  the  rough  outline  (essentially  draw 
the  framework).  This  approach  is  useful  when  a  particular  shape  is  desired. 
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A  frame  can  also  be  developed  using  fractal  functions  to  pervert  the 
lattice  into  a  controlled  random  shape  from  a  given  plane  of  the  lattice.  This  is  a 
powerful  method  that  can  be  controlled  via  bounds  on  the  random  tools, 
heuristics  or  discrete  functional  bounding  of  the  fractal  function.  This  approach  is 
most  useful  when  a  class  of  mountain  shapes  are  required  but  no  particular 
mountain  needs  to  be  modeled,  i.e.  when  random  landscapes  sufHce. 
Alternatively,  the  stick  frame  of  the  mountain  can  be  determined  via  manual 
(hand  computation)  means.  This  approach  is  tedious  and  limited  in  its  flexibility, 
and  is  not  recommended. 

b.  Surface  Texture  via  Fractal  Functions 

The  next  step  in  the  creation  of  a  mountain  is  to  provide  the 
frapfcte*  day  to  cover  our  stick  frame  model.  This  clay  is  a  fractal  function  which 
closes  the  polygons  of  the  frame  model  with  an  inductive  process  that  provides  a 
continuous  pixel  surface  for  the  entire  structure  of  the  mountain. 

The  initiator /generator  paradigm  is  used.  The  initial  set  of  initiators 
is  the  frame  described  above  and  the  generator  is  a  similar  geometrical  shape  that 
reduces  in  size  continually  until  it  becomes  the  size  of  a  pixel  and  is  mapped. 

After  the  stick  frame  of  the  mountain  is  developed,  thb  texturizing  of 
the  surface  becomes  an  automatic  process  that  terminates  when  each  geometrical 
shape  that  makes  up  the  framework  is  reduced  to  a  continuous  set  of  pixels  in  the 
pixel  set.  At  this  point,  the  mountain  exists  in  the  pixel  set  (memory)  but  roust 
be  provided  color  and  tight  to  bring  it  to  life. 

c.  Hidden  Surface  Blimination 

The  pixel  set  has  the  entire  structure  the  mountain  in  memory, 
but  we  can  project  only  a  two  dimensional  image  of  one  plane  onto  the  screen. 
The  fractal  function  which  texturizes  the  surface  docs  not  concern  itself  with 
local  computations  so  many  overlapping  pixels  are  mapped  to  the  pixel  set. 
There  are  two  reasons  then,  why  wc  need  hidden  surface  removal  (in  this  case 
better  referred  to  ss  hidden  pixel  eUmination).  First  we  have  to  eliminate  the 
back  or  hidden  sides  of  the  mouiitaiii  by  projecting  only  those  pixels  which  are 
visible  along  the  axis  of  sight  to  the  per{>endiruiar  planar  surface  of  the  display 
screen.  The  second  kind  of  hidden  pixel  removal  is  caused  by  mapped  pixels 


which  were  covered  up  by  other  recursive  fractal  descents  either  before  or  after 
the  pixel  was  mapped. 

The  removal  of  hidden  pixels  is  greatly  facilitated  by  the  use  of  the 
concept  of  the  pixel  set.  In  standard  computer  graphics  hidden  surface 
elimination,  the  programmer  is  confronted  with  graphics  primitives  which  are 
fitntiionaUjf  continuous  Euclidean  shapes.  To  effectively  remove  the  hidden  parts 
of  these  shapes  is  in  general  very  tedious  and  mathematically  complicated.  Since 
the  graphics  programmer  is  shielded  from  the  primitive  -*  pixel  mapping,  he  is 
functionally  denied  access  past  the  simplifying  abstraction^^  of  graphics  Euclidean 
primitives.  The  fractal  programmer  must  have  access  to  this  level  of  the  graphics 
mapping  and  thus  can  use  simple  techniques  to  determine  if  a  pixel  is  hidden  or 
visible. 

The  simplest  and  most  economical  means  available  to  provide  hidden 
pixel  elimination  is  through  the  use  of  Z-hnffer  algorithms.  With  this  nsethod,  the 
determination  of  whether  a  pixel  is  hidden  can  be  appended  to  the  pixel  set 
mapping  process.  The  Z  coordinate  of  a  pixel  that  is  to  be  mapped  is  checked 
against  the  Z  coordinate  of  the  pixel  currently  in  the  pixel  set  at  the  same  row 
and  column  of  the  three  dimensional  array  used  to  store  the  pixel  set.  If  the  pixel 
is  ctoeer  to  the  planar  surface  of  the  display  screen,  then  the  Z  coordinate  is 
changed  to  reflect  the  position  of  the  newly  mapped  pixel. 

The  Z-bufTer  approach,  while  powerful,  docs  limit  the  fractal 
programmer's  flexibility.  The  axis  of  sight  toward  the  mountain  must  be 
determined  prior  to  the  fractal  recursive  process  so  Uiat  the  determination  of  the 
line  through  the  (now  two  dimensional)  pixel  set  is  known.  The  Z-buffer  becomes 
an  adjacency  matrix  to  the  pixel  set  and  can  retain  information  about  forwardly 
displayed  pixels  only.  All  information  is  lost  about  other  pixels  that  were 
computed  in  the  fractal  process.  If  another  view  of  the  tnountain  is  required  then 
the  entire  pixel  set  has  to  be  recomputed  with  a  new  axis  of  sight.  If  the  fractal 
function  uses  (non>tabular)  random  techniques  then  the  mountain  varies  with 
each  view. 

**  The  kbsireclton  {)r«vid<d  by  Kueltdren  primitive*  is  m  prvwerful  one  when  the  klieranttve  of 
pixel  rnxppinf  U  eoRStdered.  Wiihout  some  powerful  mapping  tool  (such  «s  fractal  fuacitoo*),  the 
pixel  level  modeling  proccu  is  in  general  very  difTicuil. 


Most  fractal  pictures  consume  such  vast  computing  resources  that 
only  one  view  is  computed  for  a  given  picture.  As  more  requirements  for  graphics 
terrain  are  determined,  a  more  powerful  method  has  to  be  used  to  retain  all  of 
the  computed  pixels  in  the  three  dimensional  pixel  set.  This  method  requires 
that  all  pixels  be  stored  in  the  three  dimensional  array  previously  described.  The 
hidden  surface  calculation  can  then  be  performed  during  the  pixel  mapping 
operation  or  as  a  separate  calculation  that  is  performed  after  all  fractal  recursion 
has  terminated. 

As  specified  in  Chapter  3,  the  full  array  approach  requires  large 
amounts  of  memory.  This  method,  however,  allows  the  computed  fractal 
mountain  to  become  an  entity  that  can  be  manipulated  versus  an  instance  of  the 
fractal  mountain  as  above. 

Both  methods  are  viable  but  the  latter  approach  provides  more 
flexibility  for  the  programmer  whereas  the  first  approach  is  a  response  to  the 
economies  of  scale  of  data  processing.  As  new  architectures  are  developed^^  with 
capacities  geared  toward  fractal  image  computatic*'  the  first  method  can  be 
eliminated. 

d.  Illuminating  the  Mountain 

If  you  stood  on  the  dark  side  of  the  moon  without  illumination,  the 
mountains  and  craters  of  the  moon  would  not  be  visible.  They  still  exbt  however, 
just  as  our  imaginary  mountaiii  exists  in  memory.  In  order  to  vbualize  them,  we 
must  illuminate  them. 

Illumination  in  computer  graphics  is  achieved  by  varying  the  light 
intensities  of  pixels  displayed  on  the  screen.  The  color  mixture  of  these  discrete 
points  determines  the  lighting  effect  that  a  viewer  perceives.  Thb  perception  b 
not  reality  but  another  deception  caused  by  scale  and  composition.  A  lighting 
model  then,  is  one  which  b  able  to  abstract  the  essence  of  color  from  a  real  world 


An  f4<tt  nnhilffiuiv  is  on»  with  n  lufr  tnnin  mfmoty  and  puraiki  (vrocmtng  capahilUMt 
wah  K  procesaioc  (Wnwau  (wti«fv  H  is  grralrr  ikaa  l(i»  maxtatum  munis*  dtsccitl  dtsiaac*). 


object  and  transform  that  essence  into  a  set  of  color  values  (intensities)  that 
accurately  deceive  the  human  eye  via  the  graphics  medium. 

The  literature  on  computer  graphics  contains  many  lighting  models 
with  diverse  approaches  to  the  same  problem.  Many  of  these  models  (like  those  of 
hidden  surface)  concern  themselves  with  illumination  of  continuous  Euclidean 
surfaces  and  as  such,  are  not  directly  germane  to  our  study 

An  object  in  space  is  a  composition  of  basic  elements.  These  elements 
interact  with  the  physics  of  light  reflection  to  create  the  spectrum  of  light  that 
our  eyes  decode.  In  a  graphics  image,  this  process  has  to  be  simulated  with 
discrete  lighting  intensity  values  for  each  pixel.  Thus,  the  illumination  of  the 
mountain  is  a  two  step  process;  the  fractal  entities  that  are  mapped  to  the  pixel 
set  must  be  provided  with  a  basic  color,  and  these  colors  have  to  be  highlighted 
and  dimmed  by  the  lighting  algorithm. 

The  basic  color  can  be  determined  during  the  pi.xel  mapping  event  of 
the  fractal  recursion  process  or  as  a  separate  process  prior  to  or  in  conjunction 
with  the  lighting  algorithm.  This  color  can  add  realism  to  the  picture  through 
heuristics  which  the  programmer  defines.  Mc»t  mountains  are  composed  of 
different  types  of  rocks  and  flora  and  these  elements  change  at  different  altitudes. 
This  type  of  heuristic  combined  with  some  random  control  structure  (i.e.  to  vary 
the  snow  peak)  can  provide  for  improved  realism  (versus  making  the  whole 
mountain  brown).  The  process  of  determining  the  basic  color  must  be 
accomplished  prior  to  applying  the  lighting  algorithm  since  the  lighting  algorithm 
can  only  vary  the  intensities  of  an  existing  color^^.  Developing  the  process  of 
basic  color  detcmiination  is  best  accomplished  Uirough  trial  and  error.  It  is  the 
artistic  aspect  of  developing  fractal  mountains. 

The  general  process  of  computer  graphics  iiiumination  concerns  itself 
with  casting  shadows  from  one  object  to  another  given  a  direction  from  an 
imaginary  light  source  and  with  highlighting  surfaces  which  are  directly  exjKised 
to  the  source.  A  surface  is  highlighted  relative  to  the  angle  at  which  the  light 

^  Thf  Gouraud  model  (inttnulr  intefpoUKon  shading)  inuance 

^  !iince  so  many  4ivm«  cotw  mo«l«b  eaut.  the  detatU  of  color  reprcscaUUoa  are  not 
covered  here. 
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source's  rays  strike  the  surface.  This  poses  special  problems  for  fractal  surfaces 
due  to  their  discontinuity  at  every  point. 

The  process  of  illuminating  a  fractal  surface  is  best  aided  by 
divorcing  the  lighting  process  from  the  fractal  computation  process  (except  as 
noted  above).  It  is  beneficial  to  view  the  pixel  set  as  a  collection  of  pe6Mes  which 
have  size  and  position.  This  abstraction  allows  us  to  view  the  pixel  as  a 
continuous  space  that  can  block  light  (cast  shadows)  and  for  which  an  angle  of 
illumination  can  be  determined  (usually  in  conjunction  with  neighboring  pixels). 

A  well  formed  fractal  mountain  surface  is  completely  connected  (no 
space  between  adjacent  pixeb  in  the  pixel  set).  Thus  the  surface  can  also  be 
viewed  as  a  continuous  (while  very  rough)  surface  where  reflected  light  can  be 
cast  from  or  to  adjacent  pixels. 

One  lighting  model  which  fits  the  fractal  process  is  the  Tonune^- 
Spamm>  model  (Ref.  8:pp.  S78-S79]. 

This  model  views  an  obiect  as  a  collection  of  facets  which  is  each  a  perfect  reflec- 
(i.e.  does  not  absorb  light).  The  orientation  of  each  facet  is  given  by  the 
Gaussian  probability  distribution  function  (i.e.  the  smooth  surface  of  ihe  Eu* 
didean  object  is  roughed  by  the  Gaussian  relationship).  The  geometry  of  the 
face^  and  the  ajrettion  of  light  (assumi^  to  be  from  an  infinitely  distant  source^ 
so  ail  rays  are  parallel)  determines  the  intensity  and  direction  of  specular  reflec¬ 
tion  as  a  function  of  the  light  source  intensity,  the  normal  to  the  average  sur- 
face,  the  direction  to  the  light  source  and  the  direction  to  the  viewpoint. 

This  model  has  to  be  modined  to  adapt  to  the  fractal  &t?>'  method.  In  the  fractal 
method,  there  is  no  need  to  rough  the  surface  to  provide  reflection  because  the 
surface  is  by  design  roughly  textured.  A  method  of  assigning  planar  fronts  to 
each  pixel  space  has  to  be  deteruiined  and  the  geometry  of  connecting  these 
fronts  identified.  Witli  these  modificatiot^  to  the  lighting  model,  each  mdivldua} 
pixel’s  color  intensity  can  be  tnodified  for  the  increase  in  intensity  associated  with 
the  light  which  fatb  u{)on  it. 

The  model  also  allows  diffuse  reflection  (light  feHceted  from  one 
object  60  another)  which  is  criti*'nl  to  bring  out  ehtiig  of  the  fractal  image.  For 
further  iuforfimtiou  on  the  trtodd  tiic  reader  b  referred  to  llie  referesce. 


e.  Summary  of  the  Fractal  Mountain  Paradigm 

To  summarize  the  methodology  we  can  view  the  process  as  a  five 
step  process; 

-  Build  the  initiator  framettfork  or  stick  frame  model  of  the  mountain. 

-  Give  the  frame’s  surface  texture  with  fractal  functions. 

-  Remove  hidden  surfaces  (pixels)  from  the  display. 

-  Bluminate  the  surface  with  lighting  algorithms. 

-  Project  the  surface  to  the  screen. 


B.  FRACTAL  TOOLS  FOR  TERRAIN  MODELING 


The  tools  presented  in  this  section  can  be  used  in  the  creation  of  fractal 
images  within  R^.  The  list  provides  a  basic  set  of  programming  tools  to  guide  the 
creation  process. 

1.  Equations  of  the  Lattice 

The  lattice  (or  controlling  structure)  can  be  very  useful  to  the  graphics 
programmer  to  implement  heuristics  or  bounding  functions  on  the  essentially 
random  progression  of  the  fractal  figure.  The  graphics  programmer  may  wish  to 


limit  the  growth  of  the  mountain  by  implementing  a  ground  level  plane  of  the 
lattice  and  a  maximum  height  that  the  mountain  can  obtain.  He  accomplishes 
this  by  arbitrarily  assigning  another  plane  of  the  lattice  as  the  upper  bounding 
plane.  The  height  of  the  mountain  can  then  be  checked  during  any  level  of  the 
fractal  recursive  descent  against  this  fixed  plane.  The  programmer  can  then  clip 
the  height  by  adjusting  the  random  equation  that  controls  the  upward  trend  to 
tend  towards  the  ground  again.  This  is  an  example  of  a  heuristic  applied  to  the 
fractal  recursion  that  controls  the  external  qualities  of  the  function. 

A  fractal  programmer  can  use  the  lattice  to  assign  the  initial  colors  to 
the  mountain  via  a  user  designed  $e<.  of  rules.  The  lattice  aids  the  user  in  the 
implementation  of  the  rules  by  giving  reference  points  for  inclusion  of  branching 
conditions  (tree  line  to  snow  line  etc.)  and  can  be  used  in  conjunction  with 
decision  weights  to  add  a  varied  transition  between  textures,  The  determination 
of  the  initial  color  of  a  mapped  pixel  i.s  usually  a  controlled  randojn  pfOces.s,  one 
of  the  primary  m,..hods  of  control  Wing  the  lattice  or  some  derivative  thcre<>f.  As 


the  height  of  the  mountain  increases  (lattice  level),  it  becomes  increasingly  more 
likely  that  it  will  transition  to  another  texture.  This  can  be  controlled  by  adding 
the  lattice  level  as  a  factor  to  the  rule  that  decides  color. 

An  example  of  a  potential  lattice  equation  and  how  it  might  relate  to  the 
pixel  space  is  demonstrated  in  Figure  5.1.  The  actual  lattice  has  been  extended 
from  the  pixel  space  in  order  to  visually  demonstrate  how  it  relates  to  the  pixei 
set.  In  actuality,  this  is  not  the  case.  The  lattice  coincides  with  the  boundaries  of 
the  pixel  space.  Although  the  lattice  can  have  a  one-to-one  relationship  with  the 
pixei  space,  this  defeats  the  purpose  of  the  lattice  (macro  control).  By  grouping 
cubic  sets  of  pixels  into  a  well-formed  relationship,  we  can  better  implement 
heuristics  and  bounding  functions. 

As  a  lattice  example,  consider  a  pixel  space  that  is  created  by  abstracting 
the  real  world  coordinate  space  for  our  mountain  as  described  below.  We  desire  a 
real  world  space  to  be  a  cubic  area  established  by  the  box  20,000  ft. 
(z  coordinate )  by  15,000  ft.  (y )  by  20,000  ft.  (z )  This  can  be  sectioned  into  a 
lattice  by  establishing  the  increment  of  distance  between  adjacent  lattice  points 
to  be  1000  ft.  and  establishing  the  corner  lattice  point  as  (0,0,0)“. 

The  mapping  function  between  the  lattice  and  pixel  space  is  then 
straightforward.  The  size  of  the  pixel  (recall  equation  from  Chapter  3)  is 
20  000  ft 

— - - — 20  ft.  and  the  lattice  cubic  sections  contain  10*  cubic  feet  or 

1,000  ft. 

equivalently  75,000  cubic  pixels. 

The  ground  level  can  then  be  identified  as  the  2000  foot  level  and  the 
bounding  height  can  l)e  assigned  a  level  of  105W)  feet.  If  wo  wish,  we  can  make 
the  bounding  heuristic  more  realistic  by  sectioning  the  lattice  into  mountainous 
areas,  each  having  different  bounding  Icveb. 


^  A  compUuly  arbiirary  s«t  of  dimmsions,  incrvinraui  tioinU. 


The  Lattice  Control  Structure 


2.  A  Fractal  Function  for  Contouring  Mountains 

The  usual  method  for  contouring  mountains  uses  a  randomized  variation 
of  the  mid-point  displacement  method  introduced  in  Chapter  4.  The  planar 
structure  is  typically  the  triangle®*  imbedded  in  R®.  The  basic  methodology  is 
demonstrated  in  Figure  5.2.  Figure  5.2.a  shows  a  triangle  with  its  first  iteration  of 
mid-point  displacement.  This  process  continues  until  all  triangles  have  reached 
the  desired  level  of  precision.  One  completed  structure  is  demonstrated  in  Figure 

5.2. b.  Th?  precision  is  typically  lower  (pixel  level)  than  that  demonstrated  in 
Figure  5.2.b  but  it  was  terminated  at  a  higher  level  to  better  demonstrate  the 
idea.  Random  techniques  (described  below)  are  used  to  produce  the  relatively 
accurate  picture  of  a  mountain  frame  as  depicted  in  Figure  5.2.c. 

In  practice,  the  random  techniques  are  implemented  with  the  mid-point 
displacement  function  during  the  fractal  recursive  descent.  The  random 
techniques  provide  local  disorder  to  the  fractal  function  which  provides  the 
computational  structure.  Results  have  shown  that  very  little  randonmess  needs  to 
be  applied  to  the  regular  structure  of  Figure  5.2.b  to  achieve  satisfactory  results. 
The  mountains  created  for  the  film  Star  Trek:  The  Search  For  Spoek  used  a 
limited  random  number  look-up  table  consisting  of  fewer  than  300  entries  (Ref. 

f|- 

3.  The  Geom  try  for  Mid-Point  Displacement 

The  general  approach  to  building  a  fractal  shape  as  illustrated  in  Figure 

5.2. C  is  to  use  the  algorithm  of  midpoint  triangle  displacement  combined  with  a 
randomized  displacement  along  the  normal  to  the  X-Z  plane  of  a  cartesian  tliircc 
space  coordinate  system.  A  recursive  procedure  which  computes  this  relationship 
requires  as  inputs  the  points  of  the  triangle.  It  computes  the  midpoints  of  each 
line  of  the  triangle  and  inscrilies  a  triangle  inside  of  the  initiating  triangle  by 
connecting  each  midpoint.  Figure  5. 3. a.  This  process  yields  four  triangles 
coincident  with  the  plane  of  the  initiating  triangle.  When  we  fix  the  X-Z  normal 
at  any  of  the  midpoints,  we  can  displace  the  tnid[>oint  by  a  discrete  distance 
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Fig  5. 2. a  The  1st  Iteration  of  Midpoint  Displacement. 
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along  the  normal  and  determine  a  point,  Figure  5.3.b.  Since  the  normal  is  to  the 
X-Z  plane,  it  is  sufficient  to  simply  modify  the  Y  coordinate  according  to  a 
'positive  or  negative  value.  This  is  equivalent  to  displacing  the  midpoint  along  the 
X'Z  normal  up  or  down.  We  perform  this  displacement  to  each  midpoint  normal 
and  replace  the  midpoint  with  these  new  points.  This  yields  a  new  structure  that 
still  consists  of  four  triangles  but  with  each  coincident  with  a  different  plane. 
Figure  5.3.c. 

a.  Midpoint  of  a  Line  in  R 

The  determination  of  the  midpoints  of  the  lines  of  the  initiating 
triangle  is  a  simple  process  that  uses  the  equation  of  Chapter  4,  and  fixes  the 

generator  ratio  constant  at  1.  This  simplifies  the  general  equation  of: 

Xj  -h  ( Generator .ratio.constant  x  Xj) 

- ^ - 

1  +  Generator.ratio.constant 

to  the  well-known  midpoint  relationships  of: 

X,  t 

-  -T- 

Y,  + Y, 

Y».d  -  — — 

2 

Km  ^  j 

The  above  equations  completely  determine  the  midpoints  of  the  lines  formed  by 
each  endpoint  of  the  initiating  triangle. 

b.  Displacement  along  the  X-  Z  Nonnal 

The  process  of  displacing  the  midpoint  along  the  X-Z  normal  is  a 
simple  one.  We  nce<l  a  factor  such  that  the  displacement  can  obtain  a  varied 
magnitude.  This  is  best  aided  by  the  inclusion  of  a  random  variable  as  a  multiple 
of  some  scaling  factor  that  is  added  to  the  Y  coordinate  of  each  computed 
midpoint.  This  process  is  demonstrated  in  the  following  code  segment: 

Randvar  =  get  rand (HchmJ); 

PoinH(yj  «  Point l(y)  +  (Scale  *  Randvar); 
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Figure  5. 3. a.  Triangular  Midpoint  Displacement. 


Figure  3.3.C  Completed  Random  Fractal  Triangle. 
Figure  5.3.  The  Random  Midpoint  Displacement  Technique. 


A  valid  question  is,  why  the  normal  to  the  X-Z  plane?  There  are 
three  good  answers  to  this  question.  Using  the  normal  to  a  fixed  plane  simpliHes 
the  computation  (eliminates  the  need  to  perform  planar  computations  at  each 
recursive  division).  It  also  is  generally  the  direction  that  we  want  the  mountain 
to  grow.  The  most  important  reason  however,  is  related  to  the  gapping  problem 
(described  below).  With  a  fixed  direction  for  displacement,  there  is  no  need  to 
eommunieate  the  direction  of  displacement  along  the  normal  between  adjacent 
side  computations.  The  recursive  levels  that  compute  adjacent  sides  are 
functionally  discordant.  It  is  demonstrated  below  that  the  solution  to  the  gapping 
problem  (inconsistent  random  numbers)  which  creates  the  need  to  communicate 
along  discordant  recursive  levels  is  algorithmically  diflicult  to  solve  and  thus 
should  be  avoided. 

c.  The  Gapping  Problem 

One  problem  exists  for  the  midpoint  displacement  procedure  which 
utilizes  a  random  displacement  along  the  X>Z  normal  line.  It  is  indirectly  caused 
by  the  data  locality  aspect  of  the  inductive  process  of  the  recursive  fractal 
descent.  The  problem  exists  when  two  adjacent  sides  of  two  adjacent  triangles  are 
not  displaced  with  the  same  value.  Each  side  is  computed  during  independent 
levels  of  the  recursive  descent  so  there  is  no  practical  method  to  communicate  the 
random  numbers  for  the  displacement. 

The  gapping  problem  is  illustrated  in  Figure  5.4.  For  the  two 
triangles  that  are  extrapolated  from  the  structure,  there  is  an  unknown 
relationship  that  is  the  random  variable  used  to  displace  the  common  midpoint. 
If  triangle  A  uses  Rond  ^  0.3  and  triangle  B  us<^  Rond  ^  '1.07,  then  the 
displacement  for  each  adjacent  midpoint  (which  are  at  the  start  coincident)  is 
skewed  in  the  opposite  direction.  This  creates  a  §ap  in  the  fractal  landscape  that 
will  (in  ail  likelyhood)  not  be  filled  by  other  fractal  shapes  from  neighboring 
triangles.  We  need  an  algorithm  which  can  insure  that  each  midpoint  (which  is 
always  shared  by  two  triangles)  has  the  same  displacement  along  the  normal  to 
the  plane. 
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d.  Solving  the  Gapping  Problem  via  Random  Tables 

The  solution  to  the  gapping  problem  is  straightforward  if  the 
programmer  adopts  the  random  number  table  as  his  random  function 
impl^entation.  The  goal  is  to  match  adjacent  triangles  with  a  seed  or 
displacement  within  the  random  table  so  that  the  random  number  returned  is 
equivalent  for  each  coincident  midpoint. 

There  exists  a  symmetry  within  the  triangle  of  Figure  5.5.a  that 
allows  such  an  approach.  Ideally  we  want  the  point  to  be  displaced  by  the 
same  magnitude  when  triangles  Tj  and  T4  (highlighted  by  textures)  compute 
their  random  numbers  for  M,.  This  can  be  facilitated  by  the  inclusion  of  a  table 
seed  for  each  recursive  call  to  the  midpoint  displacement  routine  and  by  rotating 
the  orientation  of  the  midpoint  triangle  (the  triangle  created  by  the  three 
computed  midpoints)  labeled  T4  in  Figure  S.S.a.  This  rotation  b  performed  in 
relation  to  the  random  table  and  not  in  relation  to  the  Cartesian  space,  it  b 
accomplbhed  by  adjusting  the  order  of  the  points  in  the  recursive  call. 

The  order  of  the  points  for  triangles  T|  Tj  and  T3  are  as  described 
in  Figure  5.5.b  and  for  T4  as  described  in  Figure  5.5, c.  All  four  triangles  generate 
a  recursive  sequence  and  use  the  same  seed  to  the  random  number  table.  The 
random  numbers  retrieved  from  the  table  must  observe  the  order  of  assignment 
that  b  demonstrated  in  Figure  5.5.d.  For  example,  the  line  segment  formed  by 
the  first  two  {mints  (P|  and  P^)  input  to  the  midpoint  dbplacement  routine 
determine  the  midpoint  R,.  This  midpoint  b  assigned  the  random  dbplacement 
from  the  table  corresponding  to  the  entry  ,  aeed.  The  next  midpoint  retrieves  the 
table  entry  corre${K}nding  to  #eed  4  i  and  so  on. 

If  thb  technique  b  followed  the  sequence  of  random  numbers  will 
mateh~up  as  demonstrated  in  Figure  5,5.e-  The  recursive  calb  correspond  to  the 
code  segment  in  Figure  5.6. 
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Figure  5.4  The  G&pping  Problem. 


Solving  the  Gapping  Problem 
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Figure  5.5.  A  Solution  to  the  Gapping  Probleo. 


MainO 

{ 

LOAD  THE  INITIATING  TRUNGLE 
Seed  —  1; 

frac_triangle{Pj,P2,P3,Seed) 

} 

lrac_triangle(P|,P24’3,Seed) 

{ 

DETERMINE  DISTANCE  BETWEEN  ENDPOINTS  OF  AN  INITUTOR 

If  (DIST  <  PbcelJength) 

PlotjiixelO; 

return; 

} 

COMPUTE  THE  MIDPOINTS  (M,»M2.M3) 

ADJUST  THE  Y  COORDINATE  FOR  M,  Using  RandtabIcfSecd) 
ADJUST  THE  Y  COORDINATE  FOR  Using  Randtable{Sc€d+l) 
ADJUST  THE  Y  COORDINATE  FOR  M3  Using  RaadUble(S€€d+2) 

Seed  “  Seed  3; 

/*  Tfiangtc  T,  */ 
frac  jtriang!r(Mj  .Pj.Mj.Seed) 

/*  Tfiajigle  Tj  */ 
fracJtriangl<*(M3,M5,p3.Secd) 

Triangle  T,  */ 

frar  ^riangteiP  I  .M|  Ms^^ieed) 

/•  Triangle  T4  */ 
frac_tnangle(M3.M3,Mt«^^) 


Figure  5.G.  An  Aigoriihm  for  the  Mid{x»int  Duptaceineni  Techni<iae. 


4.  Random  (Stocastic)  Fractals 

One  common  complaint  about  computer  graphics  images  and  animations 
is  the  artificial  perfection  of  the  displayed  shapes.  Our  mind  subconsciously  rebels 
against  the  order  that  b  dbplayed,  our  expectations  about  the  rough  reality  of 
nature  are  not  satbfied.  The  use  of  randomness  in  generating  fractal  images  b 
necessary  to  approximate  the  observed  disorder  of  nature.  An  example  b  the 
Koch  curve.  Although  it  resembles  a  snowflake,  it  lacks  the  realbtic  look  that 
experience  trains  our  eyes  to  see.  In  a  mathematical  sense,  the  Koch  curve  b 
beautiful;  as  an  approximate  to  nature  it  lacks  appeal. 

To  approximate  the  rough  texture  of  nature,  we  are  ferced  to  modify  the 
well-behaved  mathematical  relationship  of  the  initiator-^ generator  in  a  controlled 
manner  to  add  variety  to  our  computed  image.  Thb  modificatbn  b  usually  by 
the  inclusion  of  a  random  variable  into  the  control  structure  within  the  fractal 
equation.  The  random  variable  must  exhibit  restraint.  It  cannot  be  allowed  to 
vary  wildly  without  structure. 

One  of  the  most  appealing  random  functions  which  provides  very 
satbfactory  results  in  fractal  images  b  the  normal  dbtribution^*.  The  normal 
dbtribution  (as  opposed  to  a  uniform  dbtribution)  approximates  the  expected 
local  disorder  in  nature  (at  least  experimentally). 

a.  The  Normal  (Gaussian)  Dbtribution 

The  normal  dbtribution  b  used  tliroughout  the  natural  sekne^  for 
many  applicatitms.  It  was  first  derivctl  as  an  empirical  result  of  the  observed 
error  about  a  true  value  that  normally  occurs  when  measurements  are  taken  a 
natural  event.  The  symmetry  that  was  observed  from  error  measurement  and 
sampling  suggested  that  there  was  a  natural  order  to  such  observations.  These 
empirical  results  spurred  natural  scientbts  and  mathematicians  to  trj’  to  Ht  a 
curve  to  the  otuiervcd  graph  that  behaves  as  probability  requires  (i.e.  the  sum 
the  area  under  the  curve  equals  unity).  Many  of  the  early  scientbts 
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referred  to  the  norma!  distribution  as  the  law  of  error  in  deference  to  its  root?  in 
experimental  natural  science. 

Many  functional  character!  ‘-'ns  of  the  normal  distribution  were 
developed^,  but  credit  is  usually  attributed  to  Carl  Frederic  Gauss  [Ref.  9:pp.  1> 
llj  who  formulated  a  least  squar<^  approach,  published  in  1809  in  Theoria  Motiu 
Corporum  Coehatium.  The  form  of  the  normal  distribution  was  not  finalized 
until  the  early  20‘^  century. 

We  take  our  definition  of  the  normal  distribution  from  [Ref.  9:pp  18), 
refer  to  Figure  5.7. 

Definilion: 

The  probability  density  function  of  a  normal  random  variable  X  is  given  by: 

M  is  the  mean,  9  is  the  standard  deviation 
and  0*  is  the  variance. 

«  *.  I  I*-#*)* 

vkhere  co<x<oq  -oo<>*<oq  and  a>0 
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me  Simos  de  Lai4are  and  Abraham  dr  Motm  Krf  9pp  tl 
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This  definition  is  the  general  case  of  the  normal  distribution.  We  are 
interested  in  the  behavior  of  the  function  and  need  a  practical  way  to  determine 
a  random  number  that  we  car  use  in  the  parameter  of  the  normal  to  the  plane  in 
the  geometrical  relationship  described  above.  To  facilitate  this,  we  simplify  the 
general  normal  distribution  to  the  well  known  standard  normal  distribution, 
illustrated  in  Figure  5.8.  The  standard  normal  distribution  function  is  the  special 
case  where  fi  =  0  and  cr  =  1.  This  reduces  the  general  equation  to  the  simplified 
equatif'n: 

1  1  xO 

The  above  functions  describe  the  behavior  of  a  normal  random 
variable.  We  need  a  function  that  returns  values  from  that  function  which  will 
observe  the  period  of  the  normal  distribution.  This  means  we  need  a  string  of 
real  numbers  t  an  assigned  range  about  a  mean  that  will  observe  the 
frequency  of  the  normal  distribution. 


j — es.arx  j 

{  -gfi .  4S» . — j 

I  - - 90 . 37% - 1 

Figure  5.8.  The  Standard  Normal  Distribution. 
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b.  Standard  Computer  Random  Functions 

Some  computer  systems  provide  a  random  number  generating 
function  which  observes  the  normal  distribution.  If  this  is  provided,  then  it  can 
be  used  directly  (after  scaling)  as  a  parameter  to  displace  the  Y  coordinate  in  the 
geometry  of  the  normal  to  the  midpoint  displacement  as  described  above. 

Many  computer  systems  only  provide  a  random  number  generating 
function  which  is  uniformly  distributed  over  an  interval  of  integers.  This  is  a 
pseudo- random  number.  Such  a  function,  when  given  a  seed,  will  produce  a 
sequence  of  numbers  distributed  over  the  fixed  interval  defined  by  that  system. 
The  interval  is  typically  proportional  to  the  maximum  integer  defined  in  the 
compilers  of  the  system.  A  normal  distribution  routine  must  then  be  defined  that 
transforms  the  uniform  random  numbers  into  random  numbers  which  behave 
according  to  the  standard  normal  distribution  function. 

There  exist  transformation  functions  that  take  a  uniform  random 
variable  distributed  over  the  interval  (0,1]  into  an  approximate  normal  random 
variable  over  —  oo  <  x  <  oo^^.  This  requires  the  uniform  random  variable  to  be 
mapped  into  the  interval  [0,1]  and  then  transformed  by  the  normal 
approximating  function. 

To  transform  a  uniform  random  variable  distributed  ox’er  an  interval 
[0,maz_int]  while  maintaining  the  distribution  density,  requires  the  following  step: 


UNF,o.„  = 


^NlP|0,woxjn(  I 
max  int 


One  commonly  used  function  that  transforms  uniform  random 
variables  into  normal  random  variables  is  found  in  [Ref.  9:pp  49).  This  function 
uses  two  uniform  variables  from  [0,l|,  denoted  UNFj  and  UNF2,  and  computes 
two  normal  random  variables,  denoted  NORMj  and  NORM2. 


NORM,  =  >/r-21og^UNFi  cos(2xUNF2) 

NORM2  =  sin(2xUNF2) 


S4 


This  is  how  most  standsj-d  system  provided  computer  subroutines  perform  the  operation. 
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Appendix  B  contains  a  C  UNIX  routine  that  implements  an  algorithm  to 
compute  the  uniform[0,l]  -*  norma/ (-00, +oo)  transformation. 

A  programmer  must  be  very  careful  when  dealing  with  random 
number  generators  from  standard  system  subroutines.  These  routines  vary  widely 
and  can  provide  good  to  barely  adequate  results.  When  the  normal 
transformation  routine  is  written,  the  programmer  must  verify  experimentally 
that  his  function  adequately  models  the  normal  distribution.  This  process  is 
illustrated  by  Figure  5.9.  Appendix  B  also  contains  experimental  results  which 
verify  the  transformation. 

The  purist  may  not  accept  the  results  displayed  in  Figure  5.9  as  an 
accurate  transformation  (there  appears  to  be  a  skew  to  the  negative  direction). 
We  must  remind  ourselves  that  we  are  trying  to  approximate  the  roughness  of 
nature  and  minor  random  skewness  will  not  deter  us.  If  the  programmer  demands 
a  better  approximation,  it  is  a  simple  process  to  expand  the  sample  space  of  the 
test  and  build  a  tabic  with  exact  proportions  by  selective  deletion  of  skew 
density. 

c.  Random  Functions  versus  Table  Driven  Methods 

The  application  of  a  random  modifier  in  the  midpoint  displacement 
technique  can  be  achieved  via  two  methods. 

-  By  invoking  the  above  function  iteratively  as  a  variable. 

-  Or  by  a  variable  returned  from  a  table  lookup  operation  from  a  random 
table. 

The  choice  of  which  method  to  use  depends  on  the  programmer’s  application  but 
each  has  its  ramifications. 

In  general,  the  table  lookup  operation  is  considerably  faster  than  the 
functional  method  but  must  by  its  definition  limit  the  amount  of  randoixmess  it 
contains.  The  major  issue  however  is  the  need  to  reproduce  a  figure  under  some 
requirement  for  fixed  terrain.  This  issue  was  the  driving  force  for  Loren  Carpenter 
from  Lucas  Film  in  determining  that  he  needed  to  use  a  table  driven  method  to 
produce  the  planet  images  for  the  film  Star  Trek:  The  Search  for  Spock  [Ref.  7). 
He  had  to  be  able  to  fix  a  space  where  the  images  of  the  actors  could  be  imposed 
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onto  the  fractal  images  and  could  not  allow  the  fixed  space  to  change  with  each 
frame  computed.  This  is  the  major  advantage  of  the  table  method.  By  retaining 
a  seed  to  a  table  of  random  numbers,  you  can  reproduce  the  sequence  of 
dbplacements  along  the  normal  during  the  fractal  recursive  descent. 

When  you  consider  the  existential  qualities  of  randomness  you  are 
confronted  with  basic  questions  about  determinism  and  order  in  the  universe.  It 
ia  not  at  all  clear  which  rules  chance.  Ih  pase,  we  can  deceive  perception 
with  a  relatively  small  table  of  ranJofft  numbers. 

The  question  of  ho\y  much  rando’:u?tesj  is  enough  to  provide  for  a 
visually  appealing  texture  is  not  completely  clears  In  [Ref.  7)  Sniith  demonstrates 
a  variety  of  shapes  computed  with  the  same  algorithm  of  Figure  S.6  using 
random  number  tables  of  different  siaes.  He  demonstrated  that  as  few  as  five 
numbers  can  suffice  to  provide  enough  local  disorder  to  give  the  viewer  the 
accetJtable  texture  of  a  mountain.  If  the  mountain  segments  are  viewed  at  the 
correct  perspective  scale,  this  perception  is  clearly  felt.  A  trained 

matheiiiatician  would  find  the  five  element  mountain  statistically  unappealing 
however.  A  true  stochastic  construction  requires  a  continuous  random  function 
rather  than  a  discrete  table  method.  As  long  as  the  goal  of  our  computations  is 
merely  to  deceive  the  graphics  viewer,  it  sutiice2  to  use  the  random  number  table. 
The  table  must  be  large  enougii  to  provide  for  an  appealing  textural  perception. 
A  complete  C  program  that  computes  a  triangular  mountain  segment  using  the 
random  displacement  midpoint  technique  is  contained  in  appendix  C. 
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VI.  SHORT  CUTS  TO  MOUNTAIN  SHAPES 


Since  the  fractal  mountain  computation  (the  full  approach  with  hidden 
surfaces  etc.)  is  so  costly  in  terms  of  resources,  it  is  important  for  us  to  consider 
shortcuts  that  can  lessen  this  burden.  This  is  best  realized  by  utilizing  the  hidden 
surface  and  curve  fitting  capabilities  that  are  provided  on  some  advanced 
graphics  systems. 

Our  goal  is  to  match  the  well  known  bicubic  surface  procedures  with  the 
structure  computed  by  the  simple  fractal  algorithms.  This  is  best  accomplished 
by  modifying  the  triangular  midpoint  displacement  technique  and  using  a 
rectangle^^  as  the  basic  geometric  building  block.  Most  of  the  cubic  surface 
algorithms  use  the  rectangular  structure  as  their  basis,  so  it  is  easier  to  adapt 
them  to  our  fractal  structure. 

When  the  fractal  algorithm  of  Figure  6.2  has  its  computations  terminated 
before  reaching  the  level  of  pixel  size,  it  yields  a  connected  rectangle  structure 
like  the  one  shown  in  Figufe  6.3.  This  structure  is  a  connected  Euclidean 
structure  that  can  be  used  m  a  tvase  on  which  other  algorithms  can  be  applied. 
Cubic  equations  cmi  fill  the  poi'ygons  to  an  arbitrary  precision  and  standard 
hidden  stjrface  algorithms  can  eliminate  the  hidden  sides  of  the  computed 
surfaces  Simple  lighting  algorithms  can  be  applied  to  the  computed  surface  to 
achieve  a  realist  ic  lijtiiting  effect^.  This  is  how  Vosa  and  Carpenter  created  their 
fFact,al  ifurfacfis  m  (ileiT.  7]. 

A.  RECTANGULAR  MIDPOINT  TBCHNiQUis 

Mociifyiug  the  triangular  midpomt  algorithm  of  chapter  5  is  a  straightforward 
process  that  introduces  no  new  matbematic’S  or  difficulties.  It  consists  of  a 
procedure  to  spUi  llie  midpoints  of  eucii  side  of  the  rectangle  and  »  procedure  to 
find  the  center  of  the  rectangle’.  From  these  fwc  points,  we  construct  four  scaled 

**  Wc  actually  use  a  non-planar  fomr  siUVd  polygon.  Wp  refer  to  ths  basic  structure  as  a  rec- 
l angle  to  simplify  the  lerminolwgy, 

**  Gourauc)  shading  fw  cxaiuiplc. 
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rectangles,  as  demonstrated  in  Figure  6.1.  The  five  shared  midpoints  of  the 
generated  rectangles  are  then  displaced  along  the  normal  to  the  X-Z  plane 
according  to  a  random  Gaussian  value.  This  process  is  exactly  the  same  as  for  the 
triangular  algorithm  of  chapter  5.  The  gapping  problem  still  exists  and  this 
requires  an  algorithm  to  rotate  the  rectangle  relative  to  the  random  number  table 
and  the  starting  seed  to  insure  that  adjacent  midpoints  are  displaced  relative  to 
the  same  random  number.  The  basic  methodology  is  displayed  in  Figure  6.1,  the 
algorithm  is  contained  in  Figure  6.2  with  sample  results  in  Figure  6.3. 
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Pomitive  Y  dispJmcvmeut  Originml  i 

COyPLETED  FRACTAL 

Figure  6.1.  The  Rectangular  Mid 
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MaiiiO 

{ 

LOAD  THE  INITIATING  RECTANGLE 
Se«d  =  1; 

fracj;ectangle(P  |  ,P24*3»P4.Seed) 

) 

fr«c_rectniigle(Pi  ,P2,P3,P4,Seed) 

{ 

DETERMINE  DISTANCE  BETWEEN  ENDPOINTS  OF  AN  INITIATOR 

If  (DIST  <  Pixel.lcngth) 

{ 

Plot_j)oint(); 

returu; 

} 

COMPUTE  THE  MIDPOINTS  (M,,M2,M3,M4,Mc) 

ADJUST  THE  Y  COORDINATE  FOR  M,  Using  Randtablo(Seed) 
ADJUST  THE  Y  COORDINATE  FOR  Mj  Using  Randtable(Seed+l) 
ADJUST  THE  Y  COORDINATE  FOR  M3  Using  Randtable(Seed+2) 
ADJUST  THE  Y  COORDINATE  FOR  M<  Using  Randtable(Seed+3) 
ADJUST  THE  Y  COORDINATE  FOR  Me  Using  Randtablc(Sced+4) 

S4H>d  =  Setni  +  5; 

/•  Rectangle  R,  */ 
frac  j[ectangle{P ,  ,M|  ,Me.M4.Secd) 

/*  Rectangle  R2  */ 
frac  jrectangle(M2,Me.M|  ,P2.Seed) 
j*  Ri*ctangle  R3  */ 
fracjrectangle(M2,Mc,M3,P3,Sc€d) 

/•  Rectangle  R4  */ 
fracjreclangIe(P<,M3,Mc,M4,Seed) 

} 


Figure  6.2.  An  Algorithm  for  the  (Rcct.)  Midpoint  Displacement  Technique. 
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B,  PARAMETRIC  CUBIC  SURFACES 


A  complete  description  of  parametric  cubic  surfaces  is  too  involved  to  be 
described  in  study.  The  theoretical  basis  of  cubic  curves  is  not  directly 
applicable  to  tV;4Ctal  geometry.  For  a  complete  description  refer  to  [Ref.  8:pp. 
514-536].  If  the  reader  is  already  familiar  with  cubic  curves  and  their  derivations, 
he  can  skip  by  the  section  on  cubic  curves  to  the  section  that  details  the 
application  of  ct  bic  surfaces.  For  any  reader  who  has  not  been  exposed  to  the 
derivations  of  parametric  equations  which  yields  cubic  curve  computational 
engines,  it  is  recommended  that  he  read  the  following  section  so  that  he  may  gain 
insight  into  the  mathematics  of  cubic  surfaces.  Detailed  knowledge  of  cubic 
curves  is  not  a  prerequisite  to  the  successful  use  of  cubic  surface  fitting  engines 
with  respect  to  fractal  surfaces.  It  is  helpful,  however,  to  understand  the 
underlying  mathematics  whenever  canned  equations  are  used. 

1.  Cubic  Curves 

The  general  method  of  cubic  curves  has  as  its  basis  that  any  continuous 
curve  in  can  be  expressed  in  parametric  form.  This  form  relates  the  points 

x,y,z  with  a  parameter  t  such  that  as  t  varies  within  some  range  of  values^’  the 
equations  solve  for  unique  points  on  the  curve.  Specifying  two  endpoints  and  two 
control  points  of  a  segment  of  the  curve  allows  us  to  define  certain  constraints  to 
be  applied  to  the  parametric  equations.  These  constraints  allow  us  to  manipulate 
the  parametric  form  of  the  equations  to  yield  a  simple  vector  product  definition 
of  that  segment.  Once  this  vector  product  is  established,  we  can  solve  for  points 
on  the  curve  by  picking  discrete  values  of  t  and  solving  for  in  turn.  This 
yields  a  discrete  approximation  of  the  curve  that  can  be  as  precise  as  needed, 
a.  Parametric  Cubic  Equations  of  a  Curve 

A  parametric  cubic  curve  is  one  for  which  the  jioints  in  {r  ,v  .s ) 
are  each  represented  os  a  third-order  (cubic)  {Kilynomial  of  some  parameter  i. 
Because  we  deal  with  a  finite  s<>gmeiit  of  a  curve,  we  limit  the  range  of  the 


(  may  vary  bclwrrn  0  and  I  for  rKaRiplr 


parameter  t  to  the  range,  0  ^  t  ^1.  This  yields  the  equations: 


x(t)  =  +  bjjt*  +  Cjjt  -f  djt 

y(t)  =  «yt^  +  byt*  +  Cyt  +  dy 

a(t)  =  a,t®  -f  b,t''*  +  c,t  +  d. 

Each  equation  can  be  expressed  as  a  vector  product  as  x(t)  b  below: 


x(t)  =  (t^  t*  t  1) 


This  vector  product  separates  the  dbtinct  parameters  of  the  parametric  equation 
into  the  unknown  coefficients  of  x(t);  parajneter  I  that  wc 

wish  to  inanipubtc.  Through  this  separation,  we  are  able  to  manipulate  them  as 
algebraic  entttiea.  If  you  multiply  the  vector  product  out,  you  find  that  the 
vector  product  is  equivalent  to  the  parametric  equation  that  precedes  it.  Denote 
thb  product  as  x(t)  =  TC\  where 

T  =  (t^  t®  t  Ij 

and 


The  vector  T  is  the  same  for  x(t),  y(t)  and  all). 

We  now  establish  constraints  (as  it  set  of  control  points)  for  the 
equation  x{t)  evaluated  at  the  bounds  of  the  range  of  the  i»arameter  f,  (i.e.  t«0 
and  t~l).  We  consider  four  equations  of  x(t)  and  its  first  derivative  x'(t)  where 
these  boundary  conditions  yield  four  known  ^loints. 


x(t)“  TC’^;  when  evaluaUrtl  ai*0.  -*  x(0)  =  {0  0  0  ijC* 
x(t)« TC^;  when  cvaluaieii  at^-l.  -*  x(l)  =  ji  i  I  ijC^ 
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and  since  the  first  derivative  of  x(t)  is: 

x'{t)  =  |3t2  2t  10]C*  =  T'C^ 

x'(t)  =  T'C„;  when  evaluated  at=0,  -♦  x‘’{0)  =(001  OjCjj 
x'(t)=T'Cjj;  when  evaluated  at=l,  -♦  x'(l)  =(321  OjCjj 

We  now  have  four  equations  that  can  be  grouped  into  a  vector  product: 


x(0) 

■ 

0  0  0  1 

x(l) 

1111 

x'{0) 

0  0  10 

x'(l) 

3  2  10 

We  recognize  that  x(0)  and  x(l)  are  the  endpoints  of  the  curve 
segment  and  x'(0)  and  x'(l)  are  coiiii>onents  of  the  tangent  vector  at  the 
endpoints  (y'(t)  and  2'(t)  are  the  other  components).  With  this  knowledge  we 
are  able  to  solve  the  left  hand  side  of  the  equation  above.  These  points  (that  we 
call  Pj  through  P^)  are  the  control  points  that  we  establish  for  curve  fitting^. 
For  a  given  curve  segment  the  control  {mints  are  ftxed.  We  rewrite  the  equations 
above  with  res{>ect  to  these  kuowm  control  points: 


0  0  0  1 

P, 

1111 

IS 

0  0  10 

X 

3  2  10 

Denote  this  equatiuci  as: 

G,  «  MC„ 

The  matrix  b  often  referred  to  as  the  geometry  of  the  cubic  curve  and  M  as 
the  basis. 

This  equation  has  the  4  by  1  row  vector  C,  as  the  only  unknown. 
The  elements  of  the  0,  vector  are  the  parameters  ^he 

**  tf  <rrr  u  MT>crt  itt«B  ihne  four  kiv  (he  user'*  tnpul  lo  our  rou- 
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parametric  equations.  We  can  solve  this  equation  for  these  parameters  and 
establish  the  parametric  equations  with  the  only  unknown  being  the  parameter  t. 
The  parameter  I  can  be  discretely  varied  over  its  range  of  0  $  t  ^1,  providing  a 
set  of  points  on  the  curve.  It  is  through  these  constraints  that  the  control  points 
eontroi  the  parametric  equations  and  produce  an  equation  thet  can  produce  a 
discretely  sampleable  curve  segment  in  three  space.  Solving  the  equation  for  is 
straightforward: 

C,  =  M  ‘G, 

Subiitituting  Cjf  into  the  equation  for  x(t}  yields^®; 

x(t)  = 

Similar  arguments  yield  the  equatioits  for  y(t)  and  z(t): 

y(tj  =  TM% 

t(t)  =  TM'*G, 

'  The  itmtrix  M  *  is  constaitt  for  ail  three  equations  and  is  usually 

denou><i  by  the  typ^  of  surf^e  that  it  relates  to  Bezier  -•  Mi,,  Hermitc  Mj, 
etc.  It  b  thiough  t.he  control  (loints  aitd  their  interaction  with  the  constraints 
that  the  models  Bezier,  B-Sp!i!fe,  Cardinal  ispline,  Ferguson  {Kermite  or  Coon’s) 
surface  etc.  modify  the  iraraitietric  equatiojis  and  provide  different  curve  fitting 
engines. 

For  each  model,  the  matrix  is  constant  throughout  sfi 

computations.  To  use  the  model  requires  the  determination  of  the  cirntrol  points 
(in  conjunction  with  how  th**y  relate  to  the  curve)  and  a  vector  tnuitiplkation 
engine.  i*ince  vector  pijieiine  computations  are  ideally  suited  to  compute.rs,  this 
nipllioti  becomes  a  fast  technology  for  curve  fitting  with  an  inXiutioe  appeal  for  a 

N 

pfogranuner. 


**  Vi>  bsvv  juM  4H«'rmiBfq  tlir  Krrtttue  mo4el  t^usltoa  for  «(t). 


b.  An  Example:  Bezier  Cubic  Curves 

We  consider  the  model  called  Bezier  [Ref.  8:pp.  514-536].  The  Bezier 
model  defines  the  position  of  the  curve’s  endpoints  and  uses  two  other  points  (not 
on  the  curve)  which  define  tangents  at  the  curve’s  endpoints  (by  the  line  segment 
joining  the  tangent  points  to  the  endpoints). 

The  matrix  M  is  derived  by  setting  the  following  constraints  (see 
Figure  6.4).  One  endpoint  of  the  segment  is  located  at  Pj: 

x(0)=  P, 

The  other  endpoint  is  located  at  P4: 

x(l)=  P4 

The  line  segment  from  P|  to  P^  defines  a  tangent  at  Pj  such  that  x'(0)  relates 
to  the  points  P^P;  as  below; 

x'(0)  =  3(P2-P,) 

And  similarly  for  the  tangent  at  P4  defined  by  PatP^: 

x'(1)  =  3(P4-P3) 

Solving  for  C„  in  terms  of  yields  the  cubic  Bazier  matrix  as: 


-1 

3 

-3 

1 

3 

-6 

3 

0 

-3 

3 

0 

0 

1 

0 

0 

0 

Hence  the  equation  for  x(t}  is: 


x(t)  »  |t’ t  l|  I 


-i  3-3  1  *  « 

3  -6  3  0 

-3  3  0  0  P3 

1  0  0  0 
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Bezier  Curve 


I  Figv;re  6.4.  An  Example  of  a  Bezier  Curve.  J 

The  process  of  creating  a  Bezier  curve  given  above  parametric 
cubic  engine  is  a  simple  process  of  computing  discrete  points  on  the  curve  by 
substituting  values  along  the  range  of  I  and  (ittmg  the  curve  by  connecting  each 
poiiit  with  a  line  seginen;.  Ihis  pco^'ides  an  approximation  to  the  curve  that  can 
be  proci-ssed  at  an  arbitrary  precision  by  incrententing  6t  with  smaller  and 
smaller  lengths. 

l‘he  process  of  shaping  a  curve  is  accomplished  by  increasing  or 
decreasing  the  two  endpoint  tangents  formed  by  the  four  control  points.  It  can  be 
viewed  intuitively  by  thinking  about  each  tangent  as  a  force  which  pulls  the 
curve  in  the  direction  of  the  tangent  until  the  force  from  the  other  endpoint 
overcomes  the  original  at  the  midpoint.  The  two  emijmint  tangents  work  against 
one  another  proijortional  to  the  distance  of  6t  from  each  er'^tipoint. 

" •  Pe/tier  Surfac<*s 

Extending  the  above  method  to  cubic  surface  sections  is  accomplished  by 
adding  a  new  parameter  s  that  we  vary  from  0  ^  s  ^  1  as  we  did  with  the 
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parameter  t  in  cubic  curves.  The  connection  between  cubic  curves  and  surfaces 
can  be  made  by  fixing  one  parameter  and  varying  the  other  over  its  range.  This 
yields  a  cubic  curve.  The  equation  b  of  the  form  x(s,t)  and  is  written  as: 

x(8,t)  =  aijS^t^  +  +  a|38^t  +  a^s® 

+  a2i8*t^  +  a228^t^  +  a23S*t  +  a248* 

+  a3j8t^  +  a328t^  +  a338t  +  a348 

+  a^|t^  a^2^^  ®44 

Written  in  the  algebraic  form: 

x(8,t)  =  SCj,T‘ 

where  S  =  [s^,s^s,l)i  T  =  [t^,t*,t,lj  and  is  the  transpose  of  the  matrix  T. 

The  complete  alge  braic  manipulation  of  the  equation  to  arrive  at  the 
equation  below  is  similar  to  the  curve  process  as  described  in  the  previous 
section.  Its  details  are  covered  in  (Ref.  8:pp.  524-536].  The  equation  for  a  Bezier 
surface  patch  b: 

x(6,t)  =  SM,,Q;,Mb‘T‘ 

where  is  the  same  matrix  as  in  the  curve  equation,  Mt,^  is  its  transpe^e  and 
Qx  the  X  component  of  sixteen  control  points  of  a  surface  patch.  Bezier 
surfaces  are  intuitive  in  their  appeal  and  serve  the  fractal  rectangular  mountain 
well.  To  apply  the  technique  to  the  mountain  of  Figure  6.3  requires  the 
application  of  a  routine  that  takes  the  non-planer  four  sided  shape  of  a  computed 
initiator  and  develops  a  connected  sixteen  point  figure  as  illustrated  in  Figure  6.5. 
The  inclusion  of  a  Bezier  subroutine  at  the  recursive  termination  event  after  thb 
figure  is  developed  matches  the  sixteen  point  figure  with  a  smooth  curve.  To 
achieve  edge  continuity  requires  that  adjacent  sides  ha%’c  the  same  four  points  in 
proi>er  juxta|>osition  in  the  sixteen  point  matrix.  Thb  b  also  demonstrated  in 
Figure  6.5.  Bezier  surfaces  guarantee  such  continuity. 
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VII.  CONCLUSIONS 


A.  DIRECTIONS  FOR  FURTHER  STUDY 

Fractal  geometry  as  an  area  of  research  is  very  new.  Because  of  thb,  there  is 
a  great  need  for  refinement  and  exploration.  What  is  known  needs  to  be  refined 
into  a  set  of  workable  techniques  with  reasonable,  simple  terminology  as  its  root. 
The  areas  that  are  unknown  need  to  be  explored  intrepidly.  With  this  goal  in 
mind,  the  following  paragraphs  quickly  review  some  areas  of  prospective  research. 
The  reader  is  invited  to  explore  their  potential. 

1.  Development  of  New  Fractal  Functional  Methods 

The  current  tools  of  fractal  functions  are  tentative  and  limited  in  their 
ability  to  yield  insight.  New  applications  of  the  recursive  initiator- generator 
paradigm  are  waiting  to  be  discovered.  This  area  of  research  Is  especially  good  for 
the  graphics  programmer  since  the  graphics  medium  b  currently  the  best  method 
for  fractal  experimentation.  As  these  new  functions  are  developed,  they  can  be 
shared,  yielding  a  glossary  of  modeling  functions  that  can  be  molded  into  a 
cohesive  theory^^.  Related  to  this  is  the  need  to  develop  a  functional  language 
(within  the  language  of  mathematics)  of  fractal  geometry  to  aid  in  the 
communication  of  ideas  and  in  the  eventual  coalescence  of  the  theory. 

2.  Fractal  Lighting  Model 

The  current  state  of  ihc  art  in  computer  graphics  lighting  models  lacks  a 
complete  model  for  the  pixel  set  paradigm  that  was  intivoduccd  in  Chapter  5. 
There  arc  a  great  tuaity  practical  applications^’  wtiich  demonstrate  succe^ul 
Ughting  techniques  but  no  pubitslied  model  exists.  This  indicates  a  piecemeal 
undisciplined  adaptation  of  the  Euclidean  based  lighting  models.  Rag  tracing 
techniques  Itxik  promising,  as  does  an  adaptation  of  the  Torrance-Siiarrdw 
lighting  model  that  was  discussed  in  chapter  5.  A  good  pixel  set  lighting  model 
would  oiien  the  avenue  of  complex  ^rrain  modeling  to  a  much  wider  audience. 

Nature'*  fractal  map? 

As  evideticed  Py  the  fractal  ptcture*  that  have  beta  pubttshed. 
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3.  Fractal  Music 


In  [Ref.  7]  Voss  demonstrates  the  application  of  fractal  recursive 
1 

techniques  to  -y  noise  and  has  produced  interesting  if  not  pleasing  tonal  results. 

It  is  safe  to  surmise  that  sound  is  a  roughly  textured  physical  phenomenon  and 
that  it  may  be  possible  to  create  or  decipher  sound  using  a  fractal  model.  Such  a 

discovery  would  aid  science  in  the  area  of  (rapid )  speech  recognition. 

/ 

4.  Fractal  Computer  Graphics  Architectures 

It  is  clear  from  our  discussion  that  new  computer  architectures  need  to 
be  developed  to  support  the  pixel  set  paradigm  and  the  computational  aspects  of 
fractal  functions.  Such  special  architectures  require  parallel  processing  capabilities 
coupled  with  vast  memory  resources.  A  real-time  fractal  terrain  image  generator 
is  one  such  architectural  possibility. 

5.  A  Better  Fractal  Definition 

Fractal  geometry  is  currently  attaining  a  wide  audience.  Because  of  that, 
it  is  time  that  trained  mathematicians  tackle  the  problems  associated  with  the 
imprecise  and  unworkable  current  definition  of  fractal  sets^*.  That  definition  us^ 
competing  definitions  of  dimension,  each  of  which  is  somewhat  difficult.  A  new 
definition  could  be  based  on  a  fractal  set’s  functional  or  statistical  qualities.  Such 
a  definition  scheme  must  provide  tools  to  further  its  workability. 


Sftdiy’,  titers  hiLii  little  attention  from  the  mathrtnktkni  rommunity,  nlihoUfh  thnt  is 
«tlsangia|.  It  is  wilb  great  timidity  tbat  ones  aceepts  fractal  geometry  wiiliout  such  scrutiny. 
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B.  CONCLUSIONS 


Fractai  geometry  is  an  old  idea  that  has  found  a  new  application  with  the 
advent  of  computer  imaging  techniques.  Its  acceptance,  has  spawned  a  great  deal 
of  research  and  has  provided  a  new  tool  to  observe  nature  through  a  different 
perspective.  We  must  be  careful  to  insure  that  our  findings  are  in  fact  valid.  We 
also  must  begin  the  coalescence  of  the  many  techniques  that  have  been  developed 
in  order  to  control  the  growth  of  this  concept  and  to  attain  true  scientific 
acceptance.  Without  this  acceptance  the  theory  will  be  criticized  (validly)  as  an 
imprecise  and  unproven  idea^^.  This  would  be  an  unfortunate  occurrence  because 
of  the  potential  that  fractal  geometry  possesses. 

It  is  the  hope  of  the  author  that  this  work  has  illuminated  the  subject  of 
fractal  geometry  and  that  it  will  aid  others  in  their  research.  The  purpose  and 
essence  of  fractal  geometry  is  based  on  simple  concepts.  The  reader  must  not  be 
overawed  by  the  current  literature  and  should  retain  his  perspective  with  a  mild 
dose  of  skepticism.  He  must  not  be  blinded  by  skepticism  though  as  the  potential 
of  fractal  geometry  has  not  yet  been  realized.  In  the  final  analysis,  we  expect 
that  even  the  skeptical  reader  will  discover  the  mathematical  beauty  and 
applicative  power  that  fractal  geometry  possesses. 


**  Thtt  of  count  U  tbf  cuncftl  sUie  of  «niun  witfa  fractal  gcocnctry. 
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APPENDIX  A.  FRACTAL  COMPUTATION  IN  R* 


The  first  routine  is  the  main  routine  which  initializes  the  data  for  the  Koch  curve 
generator  and  initiates  the  recursive  process  on  each  side  of  the  initiator  triangle. 
The  second  routine  is  the  recursive  subroutine  which  performs  the  generator 
replacement  until  the  recursive  termination  event  is  reached.  The  termination 
event  is  defmed  by  the  precision  of  the  desired  output  medium. 


KOCH.C 


/* 

This  is  the  main  program  which  controls  the  initialization  of 
the  koch  generator  parameters  and  initiates  recursive  operations 
on  each  side  of  the  initiator  triangle. 

V 


/*  Global  generator  and  initiator  data  */ 
ini  Generator  jpoints; 

/*  The  number  of  points  in  the  GENERATOR  */ 
double  Gen _angle(lU); 

/•  The  angle  formed  between  init_pointl 
and  genjioint  •/ 

double  GenjratiujlOl: 

/*  The  between  initjiointl  to  genjjoint  and 
genjioint  to  initjpointZ  */ 

double  Tan_theda{lU{; 

/*  The  tangent  of  the  angle  formed  between 
initjmiutl  and  gen^point  */ 

double  Curj[^»oint|20)(2j; 

/*  Vcrlic^  of  initiator  stfucture  */ 
iiu  Object  points^inb; 

/*  The  number  of  vertices  of  the  initiating  structure  */ 

linclude  <math.h>  /*  Standard  UNIX  tacluuc  file  for  math  library  */ 

I  define  x  (• 

#derme  y  I 
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/•  BEGIN  MAIN  PROGRAM  ♦/ 


mainO 

< 

/•  Local  variables  •/ 
int  I; 

I*  Initialize  global  variables  */ 

/*  Initial  points  of  the  INITIATORS  for  demo  */ 

Curjpoint[Oj[xj  =  4.0; 

Cur_pointjojjyj  =  3.0  +  sqrt(3.0); 

Curjoint(l)(x]  =  S.O; 

Cur_point(l)|yj  =  3.0; 

Curjointj2j(xl  =  3.0; 

Curjpoint|2||yj  =  3.0; 

I*  Remember  to  close  the  side  of  the  triangle  */ 

Cur_point|3|(x)  =  4.0; 

Cur_pointj3jjyj  =  3.0  +  $qrt(3.0); 

Object_pntsjnimb  »  3; 

Generator  j)oints  “  3; 

/*  Angle  (in  radians  formed  between  init^pomtl  and  gen^^oiut 
for  demo)  */ 

Gen  ^iglefll  «  Q.O; 

Gen _angle(2i  «  0.4712388; 

Geu^tglejs]  sa  0.0; 

/*  Ratio  of  distance  bet«k‘e<m  initjiointl  and  gcn^jK>iat{i)  and 
distance  between  geii^oint(i}  and  init ^points  */ 

Gen  j;atio(l{  »  0.5; 

Genjratioj2j  “  1.0; 
tkn_ralia’3j  «  2.0; 


JOf 


/*  Tangent  of  angle  between  initjpintl  and  gen_point(i)  */ 
for  {1=1;  I  <=  Generatorjpoints;  I++) 

Tan_^heda(I]  =  tan(Gen  jwglejlj); 

} 

/♦  BEGIN  RECURSIVE  BUILD  OF  ALL  LNITIATORS  INTO  KOCH  CURVES  7 

/•  The  Koch  curve  is  defined  in  the  infinite  but  our  recursion 
will  terminate  after  the  distance  between  points  becomes  less 
then  the  length  of  the  precisbn.  */ 

for  (I»0;  I  <  ObJect_pnts_nmb;  14-f ) 

{ 

generate(Car  point(Ij{x{,Ctir_p€>int{I|jyj, 

Cur_point[l4  i||5c],Curjpint|l4 1  jlyj); 

} 


r  END  MAIN 


GENERATE.C 


r 

This  subroutine  computes  the  generator  from  a  given  set 
of  points  in  that  dcHne  a  line  segment  which  is  the 
initiator.  The  routine  is  recursive  and  terminates  at  a  predefined 
precision  that  is  input  to  the  subroutine. 

V 


/*  External  global  generator  data;  defmed  in  main  subroutine  */ 
extern  int  Generator _polnts; 

/•  The  number  of  points  in  the  GENERATOR  */ 
extern  double  Gen _angle{lO|; 

/*  The  angle  formed  between  initjpointl 
and  gea_poiat  •/ 

extern  double  G€n_ratio|lO|: 

/*  The  between  initjpointl  to  gen^point  and 
gen_point  to  init_pomt2  */ 

extern  double  TanjthedajlO]; 

/*  The  tangent  of  the  angle  formed  between 
init^pointi  and  gen^point  •/ 

linclude  <math.h>  /*  Standard  math  include  file  for  UNIX  lib  */ 

/*  SEGIN  RECTRSiVE  PROCESS  V 

generatrix  1  I  ,X2,Y”.pr?^bioa) 

I*  Pafamrter  variablts  •/ 
double  Xl.Yi.X2.Y2.precbion; 

/*  toc&l  variabliNi  */ 
long  N: 

double  PasTay|5|l2j; 
double  GjK»isu[J0|f2|.DfST; 
double  Slope  inis  ,Slope^rp.SU>pc^en; 
double  X  JW^^.Y_perp.UJ^s*rp.b_gpn,TEMP; 
double  ten  Jhousand,one.2efi>,minus_one; 
int  Kd; 

/*  «s?sign  con^tarits  */ 

ten  iltousand  **  SOtHXbO;  me  **  1.0;  sero  *  0.0,  minus  one  * 


I*  The  Koch  curve  is  defined  in  the  infinite  but  our  recursion 
will  terminate  after  the  distance  between  points  becomes  less 
then  the  iength  of  a  pixel,  */ 

/*  Determine  distance  between  point  1  and  point  2  */ 

TEMP  -  (X2  -  Xi)*(X2  -  XI)  -f  {Y2  -  Yl)*(y2  -  Y;); 

DfST  =  sqrt(  TEMP  ); 

/•IF  DIST  less  than  the  precision  then  terminate  this 
recursion  and  begin  backtracking  */ 

if  (DIST  <  precision) 

/*  Put  your  Point  plotting  routine  here  / 
printf( "polyline  2"); 
pr  intf(  "%f  0.000000"  ,X  1 ,  Y 1 ) ; 

printf{"%f  0.000000"  »X2.Y2); 

return; 

} 

/•  Put  INITIATOR  points  one  and  two  into  the  fsrst  and  last 
{Ktints  of  the  GENERATOR  points  array  as  they  are  always 
part  of  the  generated  structure  •/ 

Gjpomt(lj{i{  =  Xl; 

G_}Kjintilif2j  ^  Yl; 

G^|ioi»t(Generatorjioints  +  2j|lJ  ==«  X2; 

Gjp«mt(GeneratorjH)mis  +  2j|2|  =  Y2; 


/*  Determine  the  slope  of  the  line  formed  by  the  init_pointl 
and  init_point2.  This  is  the  slope  of  the  INITIATOR  */ 

if  {X2  !=  XI) 

{ 

if  (Y2  !=  Yl) 

{ 

Slope Jnit  =  (Y2  -  Y1)/(X2  -  Xl); 

} 

else 

{ 

Slope_init  =  0.0; 

} 

} 

else 

{ 

/*  We  can’t  have  infinity  in  a  register 
so  settle  with  10k  */ 

Slope_init  =  tenjthousand; 

} 

/*  For  each  GENERATOR  point  (except  end  points  as  they  are  equal 
to  the  INITIATOR  end  points)  find  the  X,Y  values.  This  is 
accomplished  by  using  the  data  from  the  global  external  variables. 
The  constant  data  about  the  ratios  and  angles  between  the 
INITIATOR  and  GENERATOR  remain  the  same  regardless  of  the 
INITIATORS  length  or  position  in  EUCLIDIAN  space  */ 

for  (1=1;  I  <=  Generator_points;  I++) 

{ 

/*  Using  the  ratios  of  the  generator  perpendicular  intercept 
points  on  the  INITIATOR  determine  the  X,Y  values  of  the 
point  of  intersection  of  the  perpendicular  from  the 
GENERATOR  point  to  the  INITIATOR  line.  */ 

X_perp  =  (Xl  +  Gen_ratio(I)  *  X2)/(1.0  +  Gen_ratio(I|); 
Yj)erp  =  (Yl  +  Gen^atiojlj  *  Y2)/(1.0  +  Genj;atio|I)); 
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/*  If  the  ungie  of  the  INITIATOR  point  1  and  the  GENERATOR 
point  in  question  is  zero  then  the  GENERATOR  point  b 
coincident  with  the  INITIATOR  line  and  no  further 
calculations  are  necessary  *  j 

if  (  Gen_angie[I|  ==  zero  ) 

G_point{I  ‘-lj[r]  —  X_perp', 

G_point[I-+lj(2j  ~  Yjperp; 

} 

else 

{ 

/*  There  are  three  STATES  possible  at  this  time.  STATE  1 
where  the  slope  of  the  initiator  line  is  parallel 
to  the  X  or  Y  axis  (which  causes  havoc  with  the  line 
equations).  STATE  2  where  the  slope  of  the  line  formed 
by  the  initiator  point  1  and  the  unknown  ge  lerator  point 
is  parallel  to  the  X  or  Y  axis.  Or  STATE  3  where  no  lines 
are  parallel  to  any  axis.  */ 

/*  Determine  the  slope  of  the  line  through  the  INITIATOR 
point  1  and  the  unknown  GENERATOR  point  using  the 
tangent  of  the  Gen _Mgle  in  Init.h  */ 

Slope_gen  =  (TanJtheda(I)  +  Slope_init)/ 

(one  -  Tan_^heda(I]  *  Slopejnit); 

if  ((Slope^gen  !—  zero)  &A: 

(Slop€_^en  <  ten^housand  )) 

{ 

/*  Condition  one  of  STATE  3  */ 

/*  Determine  Y-intercepi  for  the  generator  line  */ 
b_gen  «  Yl  -  (Slope jgen  *  Xl); 

if  ((Slopejinit  =«=  zero)  j| 

(Slope  iiiit  ten^thousand)) 

{ 

/*  STATE  1  V 
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if  (Slopejnit  ==  ten_thousand  ) 

{ 

/*  STATE  1  condition  1;  INITIATOR  is  parallel 
to  the  Y  axis  */ 

G_point[I+l](2)  =  Y_perp; 

G_||oint{I+lj(l)  =  (G_point(I+ll|2]  -  b  gen)/ 

Slope  _gen; 

} 

else 

{ 

/•  STATE  1  condition  2;  INITIATOR  is  parallel 
to  the  X  axis  */ 

G_point|I+l)(l)  =  X_perp; 

G_point(I+lj{2j  ~  Slopc_gen  * 

G_point[I+l]|l)  +  b_gen; 

} 

}  /*  END  STATE  1  */ 


/•  STATE  3  •/ 

/*  Determine  slope  of  perpendicular  line  through  the 
INITIATOR  perpendicular  intercept.  */ 


Slope_perp  =  (minus_one)/SlopeJnit; 

/*  Determine  Y^intercept  for  perpendicular  line  */ 

b  perp  =  y ■  (Slope_perp  *  Xjperp); 

/*  Determine  the  X,Y  values  of  the  unknowm  GENERATOR 
point.  */ 

G_^oint(I+lj(ll  (bjperp  -  bjgen  )/ 

(81opej;en  -  Slope jperp); 

Gjpoint(I+lj(2|  ®  Slope_gcn  * 

G  j)oint(I+l)|l|  +  b_jen; 


>  /•  END  STATE  3  cond.  1  if  */ 


else 

{ 

/♦  STATE  2  */ 

Slope_perp  =  {mmus_one)/Slope_init; 
b_perp  =  Y_perp  -  (Slope_perp  *  X_perp); 

if  (Slope _gen  ==  one) 

{ 

G_point(I+l][l]  =  Xl; 

Gj>oint[I-fl]{2|  =  Slope_perp  *  G_point[l4-l][ll 
+  b_perp; 

} 

else 

{ 

G_point(I+l]j2]  =  Yl; 

G_point[I+ljjlj  =  (G_point(I+l)(2]  -  b_perp)/ 
Slope_perp; 

} 

} 

}  /♦  END  IF  */ 

}  /*  END  FOR  ♦/ 

/•  Start  recursion  on  each  line  formed  by  the  generator  from 
right  to  left  */ 

for  (J=5l;  J  <=  Generator_points  -f  J++) 

generate(G_point(J)  (l),G_point(J)  (2], 

Gjpoint  (J + 1]  [  i)  ,G_point  (J  -f  l]  (2), precision) ; 

} 

/*  END  gene-ate  */ 

) 
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APPENDIX  B:  RANDOM  NUMBER  GENERATORS 

The  routine  below  is  a  C  UNIX  UCB  implementation  of  the 
unijorm  diatribution  (0,1]  standard  normal  [-00, oc]  transformation.  It 
generates  a  500  entry  table  of  random  numbers  that  observes  the  period  of  the 
standard  normal  distribution.  Following  this  routine  are  statistics  that  verify  the 
transformation. 

RANDOM.TABLE.GENERATOR.C 

/* 

This  subroutine  will  build  a  table  in  memory  that  contains  500  random 
numbers  that  observe  the  period  of  a  standard  normal  variable 
*/ 

iliinclude  <math.h>  /*  Standard  UNIX  include  file  for  math  library  */ 

/*  External  global  variables  */ 
extern  double  RANO(500|: 

/•  BEGIN  MAIN  PROGRAM  */ 
randjable  _gen() 

{ 

/*  Local  Variables  */ 
int  I,J; 

double  UNFl,  UNF2; 
double  range, pi; 
int  factor; 


pi  =»  3.1415026535; 


/*  Determine  the  range  for  the  random  numbers  of  UNIX  UCB  */ 
range  =  2; 

for  (J=l;  J<=30;  J++) 

{ 

range  =  range  *  2; 

} 

range  =  range  -  1; 

/*  Set  the  random  number  generator  seed  */ 
srandom(475836) ; 

/*  Create  a  Table  for  500  entries  */ 

for  (1=0;  I<  500;  I  =  I  +  2) 

{ 

/*  Get  a  uniform  random  number  through  the  Unix  C  subroutine  */ 

UNFl  =  random(); 

UNF2  =  random(); 

/*  Normalize  the  uniform  random  number  to  the  interval  (0,1)  */ 

UNFl  =  UNFl  /  range; 

UNF2  =  UNF2  /  range; 

/•  Mold  the  uniform  random  variable  into  the  approximate  normal 
distribution  */ 

factor  =  1.0; 

if  (log(UNFl)  <  0.0)  factor  =  -1.0; 

RANDjl)  =  sqrt{factor  *  (2.0  *  log(UNFl)))  * 
cos  ((2*pi*UNF2)); 
nAND(I)  =  RAND(Ij  *  factor; 
factor  =  1.0; 

if  (log(UNF2)  <  0.0)  factor  =  -1.0; 

RAND{I-flj  =  sqrt(factor  *  (2.0  *  log(UNFl)))  * 
sin  ({2*pi*UNF2)); 

RAND(I+lj  -  RAND(I+lj  •  factor; 

} 

return; 


VERIFYING  STATISTICS 


The  UNIX  UCfi  operating  system's  uniform  distribution  random  number 
generating  function  spans  the  interval  defined  by  its  integer  range.  For  a  VAX 
11/780  implementation  this  is  equivalent  to  -  1  or  [0,2147483647]. 

The  random  number  seed  was  assigned  the  value  of  475836.  The  UNIX  UCB 
random  number  generator  with  a  fixed  seed  yields  a  fixed  sequence  of  numbers 
returned  from  the  function,  uniformly  distributed  over  the  range.  This  yields  a 
valuable  function  if  the  table  needs  to  be  reproduced  with  the  same  sequence 
after  transformation. 

The  table  below  shows  the  results  of  the  uniform  distribution  sequence  after 
it  was  Breezed  into  the  interval  (0,1).  These  results  show  that  the  uniform 
distribution  has  an  acceptable  distribution  over  its  range.  The  transformation 
into  [0,1]  preserves  the  distribution  from  the  original  range  ([0,2^*  -  1]). 

Analysis  of  the  normalized  uniform  random  numbers 

0.0  -»  0.1  =  52 
0.1  -♦  0.2  =  47 
0.2  0.3  =  44 
0.3  0.4  :=  49 
0.4  -*  0.5  =  49 
0.5  -  0.6  =  47 
0.6  -*  0.7  =  51 
0.7  -  0.8  =  57 
0.8  -  0,9  ^  48 
0.9  -♦  1.0  =>  56 

The  table  below  shows  the  distribution  after  the 
uniform  distribiUion  [0,1]  -*  standard  normal  [-oo,oo|  transformation  given  the 
numbers  as  described  in  the  above  table.  This  i.s  the  data  which  was  used  to 
build  Figure  S.9.  The  transformation  is  acceptable  for  the  purpose  intended,  that 
is,  to  simulate  nature’s  perceived  disorder  in  a  fractal  function. 
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Analysis  of  the  normal  (Gaussian)  random  numbers 


X  <  = 
-2.75  <  X  <= 
-2.25  <  X  <= 
-1.75  <  X  <= 
-1.25  <  X  <= 
-0.75  <  X  <= 
-0.25  <  X  <= 
0.25  <  X  <= 
0.75  <  X  <= 
1.75  <  X  <= 

1.75  <  X  <= 

2.25  <  X  <= 

2.75  <  X 


-2.75  =  5 
-2.25  =  8 
-1.75  =  16 
-1.25  =  29 
-0.75  =  56 
-0.25  =  88 
0.25  =  115 
0.75  =  87 

1.25  =  59 

1.75  =  21 

2.25  =  12 

2.75  =  4 
=  0 
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APPENDIX  C:  THE  TRIANGULAR  MOUNTAIN 


The  first  routine  is  the  main  routine  which  initializes  the  generator  data  for 
the  initiating  triangle  and  initiates  the  recursive  process.  The  second  routine  is 
the  recursive  subroutine  which  performs  the  generator  replacement  until  the 
recursive  termination  event  is  reached,  which  is  defined  by  the  precision 
parameter. 


MOUNT  AIN.C 


This  is  the  main  program  that  controls  the  initialization  of  the  triangle 
initiating  structure  and  initiates  the  recursion  on  that  triangle.  The 
recursion  will  proceed  until  the  recursive  termination  event  (defined 
by  the  precision  global  parameter) 

V 


#include  <math.h>  /*  Standard  UNIX  include  file  for  math  library  */ 

/•  Global  Defines  *( 

f  define  x  0 
f  define  y  \ 

Idcfine  z  2 

/*  Global  Tables  */ 

double  RAND{500]; 
double  Precision: 
double  Scale; 


/♦  BEGIN  MAIN  PROGRAM  */ 
main() 

{ 

/*  Local  Variables  */ 
int  I,J,K; 

double  Pl(3j,P2(3],P3(3); 
int  Seed; 

/*  Create  the  initiating  triangle  structure  *  f 

Pl(x)  =  4.5; 

Pl(y)  =  3.25; 

Pl(z]  =  0.0; 

P2[x)  =  7.0; 

P2iyj  =  3.25; 

P2[2]  0.0; 

P3(x)  =  5.75; 

P3iyj  =  3.25  +  sqrt({{2.5  *  2.5)  -  (1.25  *  1.25))); 

P3|zj  «  0.0; 

/*  Build  the  random  number  table  (appendix  B)  *{ 
rand  jt  able  ^en(); 

/•  Fractalize  until  desired  precision  */ 

Seed  =  0;  f*  Entry  seed  to  the  random  number  table  */ 
Precision  =  0.3;  /*  Recursive  termination  distance  */ 

Scale  =  0.2;  /*  Scaling  factor  for  vertical  Y  displacement 

in  the  mountain  generate  subroutine  */ 

mountain  jgencrate(P  KP2,P3,Secd) ; 


r  END  MAIN  V 

) 


GENERATE.MOUNTAIN.C 


/* 

This  is  the  subroutine  that  computes  the  four  generated  triangles  from 
initiating  triangle.  The  routine  is  recursive  and  terminates  at  a 
predefmed  precision  defined  in  the  global  parameter  Prteiaion 

r 

#include  <xnath.h>  /*  Standard  math  include  file  for  UNIX  lib  */ 

^define  x  0 
^define  y  I 
f  define  z  ? 

/•  Global  Structures  */ 

extern  double  RAND(500]; 
extern  double  Precision; 
extern  double  Scale; 

/*  BEGIN  RECURSIVE  PROCESS  */ 

mountain_generatc(Pl,P2,P3,Seed) 

I*  Parameter  variables  */ 

double  Pl(3).P2|3j.P3(4j; 
int  Seed; 


{ 


/*  Local  variables  */ 
int  l.Jf; 

double  Pgen  1 13|  ,Pgco2|3j  ,Pgen3{3j ; 
double  Pmidl(3{,Fmid2{3|,Pmid3{3|; 
double  TEMP,DISr.TVVO; 

TWO  »  2.0; 
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j*  Determine  distance  between  point  1  and  point  2  */ 


TEMP  -  (P2|xl  -  Pllx))*(P2(x]  -  Pl|x))  + 

(P2(yj  -  Pi(yl)*(P21y)  -  Pi|y))  + 

(P2|zl  -  P1N)*(P2K.  •  Pill]); 

DIST  =  8qrt{  TEMP  ); 

j*  VP  DIST  less  than  one  then  terminate  this  recursion  and 
begin  backtracking  *f 

if  (DIST  <  Precision} 

{ 

/*  Put  your  polygon  plotting  routine  here  *  f 
printf("polygon3”'i ; 
printf(''%f  %f  %r,Pl(xl,Pl[yj,Pl(z)I; 
printfC'^f  %f  %r*,P2(xi,P2jyj,P2(zl); 
printf(»%f  %f  %r\P3(xj,P3iyj,P3|zj); 
return; 

} 

/*  Manage  the  Seed  number  for  a  SIX)  entry  table  */ 

if  (Seed  >  496)  Seed  =  0; 

/*  Find  the  midpoints  of  each  triangle  leg  */ 

for  (1=30;  I<“2;  I++)  /*  0  thru  2  «>  x,y,z  */ 

{ 

PmidI{Ij  =»  (Pi(Ii  +  P2|Ij)  /  TWO; 

) 

for  (I=“0;  i<=2;  I-f +)  /*  0  thru  2  ==>  x^r^z  */ 

{ 

Piiiid2(lj  =  (P2{I1  +  P3(Ij)  /  TWO; 

} 

for  (I®'©;  I<*“2;  1+4-)  /*  0  thru  2  “>  x,y,z  */ 

{ 

Pmid3|I)  «  (P3(!l  +  P1(I))  /  TWO; 

) 
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/*  A(ijust  the  Y  coordinate  =»>  normal  from  IrX  plane  */ 

Pmidl{y]  «=  (Scale  *  RAND(Seedj)  +  Pmidl|y]; 
Pinid2(yj  =  (Scale  *  RAND(Seed-f  Ij)  +  Pmid2{yl; 
Pmidalyj  =  (Scale  *  RANDlSeed+aj)  +  PraidSly); 

Seed  =  S«d  +  1; 

j*  Recurse  on  the  triangles  according  to  the  reverse  order  rule 
for  the  interior  triangle  to  preserve  seed  order  */ 

mountain^enerate(Pmidl ,P2,  Pinid2,Seed); 
mountain^enerate(Pmid3,Pmid2,P3,  Seed); 
inountain_generate(P  1 ,  Pmid  1  ,Pinid3,Seed) ; 
mountain_generate(Pmid2,Pmid3,Pmid  1  .Seed) ; 

/*  END  generate  */ 

} 
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APPENDIX  D:  THE  RECTANGULAR  MOUNTAIN 


The  first  routine  is  the  main  routine  which  initializes  the  generator  data  for 
the  initiating  rectangular  shape  and  initiates  the  recursive  process.  The  second 
routine  is  the  recursive  subroutine  which  performs  the  generator  replacement 
until  the  recursive  termination  event  is  reached,  which  is  defined  by  the  precision 
parameter. 


RECTANGULAR.MOUNTAIN.C 


/• 

This  is  the  main  program  that  controls  the  initialization  of  the  rectangular 
initiating  structure  and  initiates  the  recursion  on  that  rectangle.  The 
recursion  will  proceed  until  the  recursive  termination  event  (denned 
by  the  precision  global  parameter) 

V 

linclude  <math.h>  /*  Standard  UNIX  include  Cte  for  math  library  */ 

/•  Global  Defines  */ 

Idefine  x  0 
^define  y  1 
#defme  s  t 

/•  Global  Tables  */ 

double  RAND|500|; 
double  Precision; 
double  Scale; 


igrnrniamjxmuamM 


/*  BEGIN  MAIN  PROGRAM  */ 

mdin() 

{ 

/*  Local  Variables  */ 
int  I,J,K; 

double  Plf3],P2[3),P3[3],P4(3); 
int  Seed; 

/*  Create  the  four  sided  polygon  initiating  structure  */ 
Pl[x]  =  2.0; 

Pljy]  =  5.0; 

Pl[z]  =  1.0; 

P2[x]  =  6.0; 

P2[yj  =  5.0; 

P2[z]  =  1.0; 

P3[x]  =  7.5; 

P3[y]  =  6.5; 

P3[z]  =  3.0; 

P4(x]  =  3.5; 

P4(y]  =  6.5; 

P4(z]  =  3.0; 

/*  Build  the  random  number  table  (appendix  B)  */ 
rand_table  _gen(); 

/*  Fractalize  until  desired  precision  */ 

Seed  “  100; 

Precision  =  0.5; 

Scale  =  0.07; 


mountain_generate(  P 1  ,P2,P3,P4  .Seed) ; 

/*  END  MAIN  */ 

} 
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RECTANGULAR.GENERATE.MOUNTAIN.C 


/• 

This  is  the  subroutine  that  computes  the  four  generated  rectangles  from 
an  initiating  rectangle.  The  routine  is  recursive  and  terminates  at  a 
predefined  precision  defined  in  the  global  parameter  Precision 

/* 

finclude  <math.h>  /*  Standard  math  include  file  for  UNIX  lib  */ 

^define  x  0 
^define  y  1 
^define  z  2 

/*  Global  Structures  */ 

extern  double  RAND[500); 
extern  double  Precision; 
extern  double  Scale; 

/♦  BEGIN  RECURSIVE  PROCESS  */ 

mountain_genefate(Pl,P2,P3,P4,Seed) 

/*  Parameter  variables  */ 
double  Pl(3l,P2{3l.P3(3l,P4[3j; 
int  Seed; 

{ 

/*  Local  variables  *  I 
int  I,J; 

double  Pmid  1  (3) ,Proid2(3j ,Pmid3(3]  ,Pnud4(3),Centor(3j; 
double  TEMP, DIST, TWO, FOUR; 

TWO  =  2.0;  FOUR  =  4.0; 

/*  Determine  distance  between  point  1  and  point  2  */ 

TEMP  =  (P2|x|  •  Pl(x))*(P2(xj  -  Pl(xj)  4- 
(P2|y]-Pl[y|)*{P2lyl-Plly))  + 
(P2(zI-Pl!zir(P2{2)-Pllzl); 

DIST  «=  sqrt.(  TEMP  ); 


I*  If  DIST  less  than  one  then  terminate  thb  recursion  and 
begin  backtracking  */ 

if  (DIST  <  Precision) 

{ 

/*  Put  your  Polygon  output  routine  here  */ 
printf(  ”polygon4”) ; 
printf("%f  %f  %P»,Pl(x],Pl(y],Pl{z)); 
printfi-'^f  %i  •,P2{xj,P2(yj,P2(zj); 
printf(”%f  %f  %f',P3|xi,P3iy],P3(zj); 
printf(''%f  %f  %P',P4(xj,P4ly],P4iz]); 
return; 


/*  Manage  the  Seed  number  for  a  500  entry  table  */ 

if  (Seed  >  496)  Seed  =  0; 

\ 

/*  Find  the  midpoints  of  each  rectangle  leg  */ 
for  (1=0;  I<=2;  I+-f )  /*  0  thru  2  =>  x,y,z  */ 

{ 

Pmidl(I|  =  (P1(I)  +  P2(I])  /  TWO; 

} 

for  (1=0;  I<=2;  I++)  /*  0  thru  2  =>  x^^  */ 

{ 

Pmid2(Ij  =  (P2{I]  +  P3{Ij)  /  TWO; 

} 

for  (1=0;  I<=2;  1++)  /*  0  thru  2  =>  x,y,z  */ 

{ 

Pmid3|Ij  =  (P3{I]  +  P4(IJ)  /  TWO; 

) 

for  (1=0;  I<=2;  I++)  /*  0  thru  2  =>  x^,z  */ 

{ 

Pmid4(Il  =  (Pl(IJ  +  P4ll|)  /  TWO; 

) 


/*  The  four  sided  polygon  is  non-planar  eo  average  the  xyz-displacement 
for  a  best  fit  approach  */ 
for  (1=0;  I<=2;  I++)  /*  0  thru  2  =>  x,y,z  */ 

{ 

Center(I]  =  (P3(I]  +  Pl[I]  +  P4[I]  +  P2(I])  /  FOUR; 

} 

j*  Adjust  the  Y  coordinate  =>  normal  from  Z-X  plane  */ 

Pmidlfy]  =  (Scale  *  RAND  (Seed))  +  Pmidl(y|; 

Pmid2[yj  =  (Scale  *  RAND{Seed+l])  +  Pmid2[y); 

Pmid3[yj  =  (Scale  *  RAND[Seed+2])  +  Pmid3{y]; 

Pniid4[y)  =  (Scale  *  RAND[Seed+3j)  +  Pmid4(yj; 

Center(y]  =  (Scale  *  RAND(Seed+4])  +  Center(yj; 

Seed  =  Seed  +  4; 

/*  Recurse  on  the  rectangles  according  to  the  reverse  order  rule 
for  the  interior  rectangles  to  preserve  seed  order  */ 

mountainjgenerate(Pl,  Pmidl,  Gen  ter  ,Pmid4, Seed); 
mountain_generate(Pmid2, Center, Pmidl,  P2,  Seed); 
mountain  generate(Pmid2, Center, Pmid3,  P3,  Seed); 
mountain_generate(P4,  Pmid3,  Center ,Pmid4, Seed); 

/*  END  generate  */ 

} 
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APPENDIX  E:  GEOMETRIC  SUPPORT 


Many  fractal  applications  and  computer  graphics  models  use  the  normal  to  a 
plane  as  a  computational  reference  point.  For  this  reason,  this  appendix  is 
devoted  to  two  tools  for  determining  the  plane  equation  of  a  polygon  and  the 
equation  of  the  normal  to  the  computed  plane. 

Determinant  Approach  to  the  Planar  Equation 

One  of  the  most  common  forms  of  a  planar  equation  is  the  general  form.  This 
form  uniquely  describes  a  plane  through  four  coefficients  A,B,C  and  D: 

Ax  +  By  +  Cz  =  D 


With  three  points  on  a  plane,  you  can  determine  the  planar  equation  by 
computing  the  coefficients.  This  approach  utilizes  the  determinant  form  of  the 
planar  equation.  Given  the  points  Pi  =  (^i»yii2i)»  P2  =  (*2»V2»^2) 

P3  “  such  that  Pi  ^  P2  ^  P3,  these  points  determine  a  unique  plane 

in  space  through  the  determinant  equation: 


X  -  X,  y  -  y,  z  -  t, 
X2  -  y2  ”  yi  *2  -  «i 
X2  -  Xj  y2  “  yi  *2  - 


0 


To  simplify  the  equation,  we  replace  the  constant  differences  by  the  expressions: 


Clx  =  *2  -  *1 

C2,  =  X3  -  X, 

Cly  =■  X2  -  y, 

C2y  «  X3  -  y, 

Cl,  =  xj  -  e, 
C2,  =  X3  —  *1 
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Evaluating  the  determinant  using  the  diagonal  approach  yields: 

((x  -  x,)ClyC2,  -  (x  -  x,)Cl,C2yl  + 

((y  -  yi)Cl,C2x  -  (y  -  y,)Cl3jC2,l  + 

((e  -  ei)Cl3jC2y  -  (e  -  *i)ClyC2j  =  0 

Solving  the  equations  for  x,y  and  z  in  terms  of  the  constant  expressions: 

A  =  ClyC2,  -  Cl,C2y 
B  =  Cl,C2x  -  C1,C2, 

C  =  Cl,C2y  -  ClyC2„ 

D  =  -  [Ax,  +  By,  +  Cz,) 

The  Normal  to  the  Plane 

Once  the  parameters  A,B  and  C  have  been  determined,  the  solution  of  the 
linear  equation  for  any  normal  to  the  plane  is  straightforward.  Using  the  plane 
parameters  in  the  parametric  equation  for  the  normal  line  to  the  plane  and  using 
any  known  point  on  the  plane  (xk^,yk,i,^,Zk^)  (the  midpoint  of  the  fractal 
triangle  for  example)  determines  a  normal  line  as: 

X  =  +  c  A 

Y  =  Vkwn  +  c  B 
Z  =  +  c  C 

where  c  is  a  parameter  such  that  c  is  an  element  of  R.  By  varying  the  parameter 
c  we  can  solve  for  unique  points  on  the  normal  line  to  the  plane. 
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